Thanks for your suggestions.
Generally, in FEM you number elements like in the below picture: <https://lh3.googleusercontent.com/-0n-HlF1EsRs/WLBEmN34aaI/AAAAAAAAAGU/GlK9R9OqDJ8FSTLKGXeaI85eU6pdV3MkwCLcB/s1600/System.png> I understand of course, that the cells (elements) in deal.II are parallely distributed, but it should be possible to identify the geometrical position and ID of the elements and their corresponding nodes. For instance the following works: typename parallel::distributed::Triangulation<dim>::active_cell_iterator cell = triangulation.begin_active(), endc = triangulation.end(); for (; cell != endc; ++cell) { std::cout << "CELL NUMBER: " << cell->id() << std::endl; } Of course here the locally owned term is missing, but is it correct to access the cells like this? I am a bit confused since it is a rather easy task which should not be such a big "deal" in deal.II. Best regards, S. A. Mohseni -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
