From cb2cbe073ae3a30234e0f85563e0c5c98fbc55ef Mon Sep 17 00:00:00 2001 From: hartmann Date: Thu, 28 Oct 2004 14:19:28 +0000 Subject: [PATCH] Fix check for neighbor_child_on_subface git-svn-id: https://svn.dealii.org/trunk@9730 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/bits/mesh_3d_15.cc | 30 ++++++++---------------------- 1 file changed, 8 insertions(+), 22 deletions(-) diff --git a/tests/bits/mesh_3d_15.cc b/tests/bits/mesh_3d_15.cc index 30eb28a3af..b77239e3d1 100644 --- a/tests/bits/mesh_3d_15.cc +++ b/tests/bits/mesh_3d_15.cc @@ -62,28 +62,14 @@ void check_this (Triangulation<3> &tria) const unsigned int neighbor_neighbor = cell->neighbor_of_neighbor (face_no); - // see whether face and - // the neighbor's - // counterface share the - // same indexing of - // children. if not so, - // translate child - // indices - const bool face_orientations_match - = (neighbor->face_orientation(neighbor_neighbor) == - cell->face_orientation(face_no)); - static const unsigned int subface_translation[4] - = { 0, 3, 2, 1 }; - const unsigned int neighbor_child_index - = (GeometryInfo<3>:: - child_cell_on_face(neighbor_neighbor, - (face_orientations_match ? - subface_no : - subface_translation[subface_no]))); - const DoFHandler<3>::active_cell_iterator neighbor_child - = neighbor->child(neighbor_child_index); - - Assert (neighbor_child == + const DoFHandler<3>::face_iterator subface= + cell->face(face_no)->child(subface_no); + unsigned int i=0; + for (; i::children_per_cell; ++i) + if (neighbor->child(i)->face(neighbor_neighbor)==subface) + break; + + Assert (neighbor->child(i)== cell->neighbor_child_on_subface (face_no, subface_no), ExcInternalError()); -- 2.39.5