[petsc-users] How to combine different element types into a single DMPlex?

[Eric Chamberland]

On 2024-08-01 09:24, Matthew Knepley wrote:
> On Thu, Aug 1, 2024 at 8:23 AM Eric Chamberland 
> <Eric.Chamberland at 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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20240801/053fb210/attachment.html>