Dear PETSc developers and users, I am currently exploring the integration of MMG3D with PETSc. Since MMG3D supports only serial execution, I am planning to combine parallel and serial computing in my workflow. Specifically, after solving the linear systems in parallel using PETSc:
1. I intend to use DMPlexGetGatherDM to collect the entire mesh on the root process for input to MMG3D. 2. Additionally, I plan to gather the error field onto the root process using VecScatter. However, I am concerned that the nth value in the gathered error vector (step 2) may not correspond to the nth element in the gathered mesh (step 1). Is this a valid concern? Do you have any suggestions or recommended practices for ensuring correct correspondence between the solution fields and the mesh when switching from parallel to serial mode? Thanks, Xiaodong
