if ((cell->neighbor(f)->level() == cell->level())
&&
(cell->neighbor(f)->has_children() == false)
- &&
+ &&
cell_is_in_fluid_domain (cell->neighbor(f)))
{
// same size
&&
(cell->neighbor(f)->has_children() == true))
{
- // neighbor has children. loop over
+ // neighbor has children. loop over
// the cells adjacent to the commone
// interface and see which subdomain
// they belong to
for (unsigned int subface=0; subface<cell->face(f)->n_children(); ++subface)
if (cell_is_in_fluid_domain (cell->neighbor_child_on_subface (f, subface)))
{
- elasticity_fe_subface_values.reinit (cell,
+ elasticity_fe_subface_values.reinit (cell,
f,
subface);
stokes_fe_face_values.reinit (cell->neighbor_child_on_subface (f, subface),
local_dof_indices,
neighbor_dof_indices,
system_matrix);
-
+
}
}
}
(((cell->neighbor(f)->level() == cell->level())
&&
(cell->neighbor(f)->has_children() == false)
- &&
+ &&
cell_is_in_fluid_domain (cell->neighbor(f)))
||
((cell->neighbor(f)->level() == cell->level())
(((cell->neighbor(f)->level() == cell->level())
&&
(cell->neighbor(f)->has_children() == false)
- &&
+ &&
cell_is_in_solid_domain (cell->neighbor(f)))
||
((cell->neighbor(f)->level() == cell->level())
))
estimated_error_per_cell(cell_index) = 0;
}
- }
-
+ }
+
GridRefinement::refine_and_coarsen_fixed_number (triangulation,
estimated_error_per_cell,
0.3, 0.0);