Hello Leo, I think you are right, the global_active_cell_index() depends on the partitioning of the mesh. You need to check but I believe that id() <https://dealii.org/current/doxygen/deal.II/classCellAccessor.html#aee7d5ec9cac8af95a136554501ea4f0a> does what you want.
Best, Bruno On Sunday, September 21, 2025 at 5:09:17 AM UTC-4 [email protected] wrote: Deal all, I would like to confirm that if the cell_index is different while running with different number of MPI cores, that is the reordering of triangulation cell would differ according to the MPI cores? I just found that when I want to output the belonged cell_index and quadrature point position in real space of some boundary, the qudrature points are identical in general if I did the x-y plot of all the poins but the belonged cell_index are different while using serial running or different number of MPI cores. The data were extracted using cell-> global_active_cell_index() to get cell_index and fe_face_values_ref. quadrature_point(f_q_point) to get the quadratuer_point f_q_point in FEFaceValue fe_face_values_ref. The grid was not generated by deal.ii functions but read from ABAQUS .inp file. If so, is there any way to keep the reordering identical in all the cases? Best Regards, Leo -- 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 [email protected]. To view this discussion visit https://groups.google.com/d/msgid/dealii/3c14816b-5dc4-4e3f-9751-57e7f4bc8952n%40googlegroups.com.
