[petsc-users] Mixing PETSc Parallelism with Serial MMG3D Workflow

[neil liu]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20250418/f46553d3/attachment.html>