Dear PETSc team:
Greetings.
I keep working on mesh I/O utilities using DMPlex.
Specifically for the output stage, I need a solid grasp on the global numbers
and ideally how to set them into the DMPlex during an input operation and
carrying the global numbers through API calls to DMPlexDistribute() or
DMPlexMigrate() and hopefully also through some of the mesh adaption APIs. I
was wondering if the GlobalToNatural PetscSF manages these global numbers. The
next most useful object is the PointSF, but to me, it seems to only help
establish DAG point ownership, not DAG point global indices.
Otherwise, I have been working with the IS obtained from
DMPlexGetPointNumbering() and manually determining global stratum sizes,
offsets, and numbers by looking at the signs of the involuted index list that
comes with that IS. It's working for now (I can monolithically write meshes to
CGNS in parallel), but it is resulting in repetitive code that I will need for
another mesh format that I want to support.
Sincerely:
J.A. Ferrand
Embry-Riddle Aeronautical University - Daytona Beach - FL
Ph.D. Candidate, Aerospace Engineering
M.Sc. Aerospace Engineering
B.Sc. Aerospace Engineering
B.Sc. Computational Mathematics
Phone: (386)-843-1829
Email(s): ferra...@my.erau.edu
jesus.ferr...@gmail.com
