From 5003efdf5420efb895ddcdd2986afaa117132ea1 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Tue, 27 Apr 2021 15:47:51 -0600 Subject: [PATCH] Make Triangulation::get_boundary_ids() work for non-hypercube cells. --- source/grid/tria.cc | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/source/grid/tria.cc b/source/grid/tria.cc index e2dd694c26..a933df9a37 100644 --- a/source/grid/tria.cc +++ b/source/grid/tria.cc @@ -10370,7 +10370,7 @@ Triangulation::get_boundary_ids() const std::set b_ids; for (auto cell : active_cell_iterators()) if (cell->is_locally_owned()) - for (const unsigned int face : GeometryInfo::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 boundary_ids(b_ids.begin(), b_ids.end()); -- 2.39.5