Make Triangulation::get_boundary_ids() work for non-hypercube cells.

diff --git a/source/grid/tria.cc b/source/grid/tria.cc
index e2dd694c26f557dc6b3524d404758b5377f6d2c1..a933df9a37c73f19a1b3e57d373d562043ce65d6 100644 (file)
--- a/source/grid/tria.cc
+++ b/source/grid/tria.cc
@@ -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());