From: Tobias Leicht Date: Wed, 13 Jan 2010 13:52:46 +0000 (+0000) Subject: Avoid using level information explicitly as this might go wrong if anisotropic refine... X-Git-Tag: v8.0.0~6636 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=502f48a216f956c275ea0ff035103ccaa14ab8e7;p=dealii.git Avoid using level information explicitly as this might go wrong if anisotropic refinement has been used. git-svn-id: https://svn.dealii.org/trunk@20358 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/fe/fe_poly_tensor.cc b/deal.II/deal.II/source/fe/fe_poly_tensor.cc index fcb832e152..db97ccb8b3 100644 --- a/deal.II/deal.II/source/fe/fe_poly_tensor.cc +++ b/deal.II/deal.II/source/fe/fe_poly_tensor.cc @@ -64,20 +64,7 @@ namespace Triangulation::face_iterator face = cell->face (f); if (!face->at_boundary ()) { - const unsigned int neighbor_level = cell->neighbor (f)->level (); - const unsigned int cell_level = cell->level (); - - // check of neighbor is more - // refined (this case should - // never be encountered) - Assert (neighbor_level <= cell_level, - ExcMessage ("neighbor_level is larger than cell level!")); - - const unsigned int nn = (neighbor_level == cell_level - ? - cell->neighbor_of_neighbor (f) - : - cell->neighbor_of_coarser_neighbor (f).first); + const unsigned int nn = cell->neighbor_face_no(f); if (nn < GeometryInfo::faces_per_cell / 2) for (unsigned int j = 0; j < dofs_per_face; ++j)