cell = dof_handler.begin_active(),
endc = dof_handler.end();
for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
+ for (const unsigned int f : cell->face_indices())
if (cell->at_boundary(f))
{
- const unsigned int dofs_per_face = cell->get_fe().dofs_per_face;
+ const unsigned int dofs_per_face = cell->get_fe().n_dofs_per_face();
dofs_on_face.resize(dofs_per_face);
cell->face(f)->get_dof_indices(dofs_on_face,
cell->active_fe_index());
if (interior_only)
{
// Exclude degrees of freedom on faces opposite to the vertex
- exclude.resize(fe.dofs_per_cell);
+ exclude.resize(fe.n_dofs_per_cell());
std::fill(exclude.begin(), exclude.end(), false);
- const unsigned int dpf = fe.dofs_per_face;
+ const unsigned int dpf = fe.n_dofs_per_face();
- for (const unsigned int face : GeometryInfo<dim>::face_indices())
+ for (const unsigned int face : cell->face_indices())
if (cell->at_boundary(face) ||
cell->neighbor(face)->level() != cell->level())
for (unsigned int i = 0; i < dpf; ++i)
// Identify all vertices active on this level and remember some data
// about them
for (cell = dof_handler.begin(level); cell != endc; ++cell)
- for (const unsigned int v : GeometryInfo<dim>::vertex_indices())
+ for (const unsigned int v : cell->vertex_indices())
{
const unsigned int vg = cell->vertex_index(v);
- vertex_dof_count[vg] += cell->get_fe().dofs_per_cell;
+ vertex_dof_count[vg] += cell->get_fe().n_dofs_per_cell();
++vertex_cell_count[vg];
for (unsigned int d = 0; d < dim; ++d)
{
for (cell = dof_handler.begin(level); cell != endc; ++cell)
{
const FiniteElement<dim> &fe = cell->get_fe();
- indices.resize(fe.dofs_per_cell);
+ indices.resize(fe.n_dofs_per_cell());
cell->get_mg_dof_indices(indices);
- for (const unsigned int v : GeometryInfo<dim>::vertex_indices())
+ for (const unsigned int v : cell->vertex_indices())
{
const unsigned int vg = cell->vertex_index(v);
const unsigned int block = vertex_mapping[vg];
cell = dof.begin_active(),
endc = dof.end();
for (; cell != endc; ++cell)
- for (const unsigned int f : GeometryInfo<dim>::face_indices())
+ for (const unsigned int f : cell->face_indices())
if (cell->at_boundary(f))
{
- const unsigned int dofs_per_face = cell->get_fe().dofs_per_face;
+ const unsigned int dofs_per_face = cell->get_fe().n_dofs_per_face();
dofs_on_this_face.resize(dofs_per_face);
cell->face(f)->get_dof_indices(dofs_on_this_face,
cell->active_fe_index());