typename DH::cell_iterator cell = dof.begin(level),
endc = dof.end(level);
for (; cell!=endc; ++cell)
- if (dof.get_tria().locally_owned_subdomain()==numbers::invalid_subdomain_id || cell->level_subdomain_id()==dof.get_tria().locally_owned_subdomain())
+ if (dof.get_tria().locally_owned_subdomain()==numbers::invalid_subdomain_id
+ || cell->level_subdomain_id()==dof.get_tria().locally_owned_subdomain())
{
cell->get_mg_dof_indices (dofs_on_this_cell);
// make sparsity pattern for this cell
endc = dof.end();
for (; cell!=endc; ++cell)
{
- if (cell->level_subdomain_id()!=dof.get_tria().locally_owned_subdomain())
+ if (cell->level_subdomain_id()!=dof.get_tria().locally_owned_subdomain()
+ && dof.get_tria().locally_owned_subdomain()!=numbers::invalid_subdomain_id)
continue;
const FiniteElement<dim> &fe = cell->get_fe();
const unsigned int level = cell->level();
cell = dof.begin(),
endc = dof.end();
for (; cell!=endc; ++cell)
+ if (cell->level_subdomain_id()==dof.get_tria().locally_owned_subdomain()
+ || dof.get_tria().locally_owned_subdomain()==numbers::invalid_subdomain_id)
for (unsigned int face_no = 0; face_no < GeometryInfo<dim>::faces_per_cell;
++face_no)
{
for (; cell!=endc; ++cell)
{
- if (cell->level_subdomain_id()!=mg_dof_handler.get_tria().locally_owned_subdomain())
+ if (cell->level_subdomain_id()!=mg_dof_handler.get_tria().locally_owned_subdomain()
+ && mg_dof_handler.get_tria().locally_owned_subdomain()!=numbers::invalid_subdomain_id)
continue;
std::fill (cell_dofs.begin(), cell_dofs.end(), false);
for (; cell!=endc; ++cell)
{
- if (cell->level_subdomain_id()!=mg_dof_handler.get_tria().locally_owned_subdomain())
+ if (cell->level_subdomain_id()!=mg_dof_handler.get_tria().locally_owned_subdomain()
+ && mg_dof_handler.get_tria().locally_owned_subdomain()!=numbers::invalid_subdomain_id)
continue;
std::fill (cell_dofs.begin(), cell_dofs.end(), false);
for (; cell!=endc; ++cell)
{
- if (cell->level_subdomain_id()!=mg_dof_handler.get_tria().locally_owned_subdomain())
+ if (cell->level_subdomain_id()!=mg_dof_handler.get_tria().locally_owned_subdomain()
+ && mg_dof_handler.get_tria().locally_owned_subdomain()!=numbers::invalid_subdomain_id)
continue;
bool has_coarser_neighbor = false;