Reorder differently the triangulation on different MPI cores?

11 views
Skip to first unread message

Léonhard YU

unread,
Sep 21, 2025, 5:09:17 AM (2 days ago) Sep 21
to deal.II User Group
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

Bruno Turcksin

unread,
Sep 21, 2025, 9:30:37 PM (2 days ago) Sep 21
to deal.II User Group
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() does what you want.

Best,

Bruno
Reply all
Reply to author
Forward
0 new messages