/**
* Sort the degrees of freedom by component. It does the same thing as the
- * above function, only that it does this for one single level of a multi-
- * level discretization. The non-multigrid part of the the DoFHandler is not
- * touched.
+ * above function, only that it does this for one single level of a multilevel
+ * discretization. The non-multigrid part of the the DoFHandler is not touched.
*/
template <class DH>
void
/**
* Sort the degrees of freedom by vector block. It does the same thing as the
- * above function, only that it does this for one single level of a multi-
- * level discretization. The non-multigrid part of the the DoFHandler is not
- * touched.
+ * above function, only that it does this for one single level of a multilevel
+ * discretization. The non-multigrid part of the the DoFHandler is not touched.
*/
template <int dim, int spacedim>
void
- block_wise (DoFHandler<dim,spacedim> &dof_handler, unsigned int level);
+ block_wise (DoFHandler<dim,spacedim> &dof_handler, const unsigned int level);
/**
* Sort the degrees of freedom by block. It does the same thing as the above
if (is_level_operation)
{
//we are dealing with mg dofs, skip foreign level cells:
- if ((start->get_dof_handler().get_tria().locally_owned_subdomain() != numbers::invalid_subdomain_id)
- &&
- (cell->level_subdomain_id()!=start->get_dof_handler().get_tria().locally_owned_subdomain()))
+ if (!cell->is_locally_owned_on_level())
continue;
}
else
compute_block_wise (std::vector<types::global_dof_index> &new_indices,
const ITERATOR &start,
const ENDITERATOR &end,
- bool is_level_operation)
+ const bool is_level_operation)
{
const hp::FECollection<dim,spacedim>
fe_collection (start->get_dof_handler().get_fe ());
if (is_level_operation)
{
//we are dealing with mg dofs, skip foreign level cells:
- if ((start->get_dof_handler().get_tria().locally_owned_subdomain() != numbers::invalid_subdomain_id)
- && (cell->level_subdomain_id()
- != start->get_dof_handler().get_tria().locally_owned_subdomain()))
+ if (!cell->is_locally_owned_on_level())
continue;
}
else