]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Make Triangulation::get_boundary_ids() work for non-hypercube cells.
authorWolfgang Bangerth <bangerth@colostate.edu>
Tue, 27 Apr 2021 21:47:51 +0000 (15:47 -0600)
committerWolfgang Bangerth <bangerth@colostate.edu>
Wed, 28 Apr 2021 03:14:24 +0000 (21:14 -0600)
source/grid/tria.cc

index e2dd694c26f557dc6b3524d404758b5377f6d2c1..a933df9a37c73f19a1b3e57d373d562043ce65d6 100644 (file)
@@ -10370,7 +10370,7 @@ Triangulation<dim, spacedim>::get_boundary_ids() const
       std::set<types::boundary_id> b_ids;
       for (auto cell : active_cell_iterators())
         if (cell->is_locally_owned())
-          for (const unsigned int face : GeometryInfo<dim>::face_indices())
+          for (const unsigned int face : cell->face_indices())
             if (cell->at_boundary(face))
               b_ids.insert(cell->face(face)->boundary_id());
       std::vector<types::boundary_id> boundary_ids(b_ids.begin(), b_ids.end());

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.