// We loop over
// subfaces
for (unsigned int subface_no=0;
- subface_no<GeometryInfo<dim>::subfaces_per_face;
- ++subface_no)
+ subface_no<face->n_children(); ++subface_no)
{
// and set the
// cell
cell->neighbor_of_neighbor(face_no);
for (unsigned int subface_no=0;
- subface_no<GeometryInfo<dim>::subfaces_per_face;
- ++subface_no)
+ subface_no<face->n_children(); ++subface_no)
{
typename DoFHandler<dim>::cell_iterator neighbor_child
= cell->neighbor_child_on_subface (face_no, subface_no);
const unsigned int
n_q_points = face_data.fe_face_values_cell.n_quadrature_points;
+ const typename DoFHandler<dim>::face_iterator
+ face = cell->face(face_no);
const typename DoFHandler<dim>::cell_iterator
neighbor = cell->neighbor(face_no);
Assert (neighbor.state() == IteratorState::valid,
// for one face for all the
// sub-faces now:
for (unsigned int subface_no=0;
- subface_no<GeometryInfo<dim>::subfaces_per_face;
- ++subface_no)
+ subface_no<face->n_children(); ++subface_no)
{
// Start with some checks
// again: get an iterator
// computed and do not carry an
// invalid value.
double sum = 0;
- typename DoFHandler<dim>::face_iterator face = cell->face(face_no);
for (unsigned int subface_no=0;
- subface_no<GeometryInfo<dim>::subfaces_per_face;
- ++subface_no)
+ subface_no<face->n_children(); ++subface_no)
{
Assert (face_integrals.find(face->child(subface_no)) !=
face_integrals.end(),
// We loop over
// subfaces
for (unsigned int subface_no=0;
- subface_no<GeometryInfo<dim>::subfaces_per_face;
- ++subface_no)
+ subface_no<face->n_children(); ++subface_no)
{
// and set the
// cell
cell->neighbor_of_neighbor(face_no);
for (unsigned int subface_no=0;
- subface_no<GeometryInfo<dim>::subfaces_per_face;
- ++subface_no)
+ subface_no<face->n_children(); ++subface_no)
{
typename hp::DoFHandler<dim>::cell_iterator neighbor_child=
neighbor->child(GeometryInfo<dim>::child_cell_on_face(