From: Wolfgang Bangerth Date: Tue, 27 Apr 2021 21:47:51 +0000 (-0600) Subject: Make Triangulation::get_boundary_ids() work for non-hypercube cells. X-Git-Tag: v9.3.0-rc1~166^2~1 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=5003efdf5420efb895ddcdd2986afaa117132ea1;p=dealii.git Make Triangulation::get_boundary_ids() work for non-hypercube cells. --- 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());