* appropriately.
*
* It is necessary that the finite element underlying the MeshType used has
- * degrees of freedom that are logically associated to a vertex, line, quad,
+ * degrees of freedom that are logically associated to a vertex, line, quad,
* or hex. This includes both nodal degrees of freedom and other modal types
- * of dofs that are associated to an edge, etc. The result is the patch of
- * cells representing the support of the basis element associated to the
+ * of dofs that are associated to an edge, etc. The result is the patch of
+ * cells representing the support of the basis element associated to the
* degree of freedom. For instance using an FE_Q finite element, we obtain
- * the standard patch of cells touching the degree of freedom and then others
+ * the standard patch of cells touching the degree of freedom and then others
* that take care of possible hanging node constraints. Using a FE_DGQ finite
* element, the degrees of freedom are logically considered to be "interior" to
- * the cells so the patch would be the cell on which the degree of freedom is
+ * the cells so the patch would be the cell on which the degree of freedom is
* located.
*
* @tparam MeshType The MeshType should be a DoFHandler<dim> or hp::DoFHandler<dim>.
// be added on account of adaptivity hanging node
// constraints.
typename MeshType::active_cell_iterator cell = dof_handler.begin_active(),
- endc = dof_handler.end();
+ endc = dof_handler.end();
for (; cell!=endc; ++cell)
{
// Need to loop through all cells that could
if (MeshType::dimension == 3)
- {
- // finally, restore user flags that were changed above
- // to when we constructed the pointers to parent of lines
- // Since dof_handler is const, we must leave it unchanged.
- const_cast<dealii::Triangulation<MeshType::dimension,MeshType::space_dimension> &>(dof_handler.get_triangulation()).load_user_flags (user_flags);
- }
+ {
+ // finally, restore user flags that were changed above
+ // to when we constructed the pointers to parent of lines
+ // Since dof_handler is const, we must leave it unchanged.
+ const_cast<dealii::Triangulation<MeshType::dimension,MeshType::space_dimension> &>(dof_handler.get_triangulation()).load_user_flags (user_flags);
+ }
// Finally, we copy map of sets to
// map of vectors using assign()