]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Avoid using level information explicitly as this might go wrong if anisotropic refine...
authorTobias Leicht <tobias.leicht@dlr.de>
Wed, 13 Jan 2010 13:52:46 +0000 (13:52 +0000)
committerTobias Leicht <tobias.leicht@dlr.de>
Wed, 13 Jan 2010 13:52:46 +0000 (13:52 +0000)
git-svn-id: https://svn.dealii.org/trunk@20358 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_poly_tensor.cc

index fcb832e152a1420a60c11a692e535e8e8610966c..db97ccb8b34ae52799c4bc75665bd05d94a76a81 100644 (file)
@@ -64,20 +64,7 @@ namespace
        Triangulation<dim,spacedim>::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<dim>::faces_per_cell / 2)
              for (unsigned int j = 0; j < dofs_per_face; ++j)

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.