On 2024-08-01 09:24, Matthew Knepley wrote:
On Thu, Aug 1, 2024 at 8:23 AM Eric Chamberland
<eric.chamberl...@giref.ulaval.ca> wrote:
Hi Matthew,
we have our own format that uses MPI I/O for the initial read,
then we would like to do almost exactly what we do in ex47.c
(https://urldefense.us/v3/__https://petsc.org/main/src/dm/impls/plex/tests/ex47.c.html__;!!G_uCfscf7eWS!Yl2BQr5WaJV41Sq7-i2xoMTi_ZGsBeThe3GPDdLjQmRtNXOdQJKpIg1Ec8-av5NcnywNIyr2D9ew6B-O8jC5ICPpWzcZ0mNNE3n3bYIy$
)
excepted the very beginning of the program that will read (MPI
I/O) from the disk. Then, always in parallel:
1- Populate a DMPlex with multiple element types (with a variant
of DMPlexBuildFromCellListParallel ? do you have an example of this?)
...
We can do that. We only need to change step 1. I will put it on my
TODO list. My thinking is the same as Vaclav, namely to replace
numCorners with a PetscSection describing the cells[] array. Will that
work for you?
Hi Matthew,
That sounds fine for me! I can create a mixed mesh partition
description so we add it to ex47.c...
I'll ping @you in a MR for that...
thanks a lot!
Eric
