Nicola,
I would simply do this by
std::mapline(0)), unsigned int> line_count;
for(const auto &cell: tria.active_cell_iterators()) {
for(unsigned int i=0; i< GeometryInfo::lines_per_cell; ++i) {
const it = line_count.find(cell->line(i));
if(it != line_cou
Dear all,
In my code I need to reconstruct the number of cells sharing an edge.
Is there a way to reconstruct this number locally? In principle I would
need something like to cell->line(1)->n_sharing_cells()?
Bests,
Nicola
--
The deal.II project is located at http://www.dealii.org/
For maili