[petsc-dev] DMPlex overlap redistribution

[Adrian Croucher]

hi Matt,

On 7/03/24 2:55 am, Matthew Knepley wrote:
> Let's do this. Can you send me a dual porosity mesh, and I will figure 
> out why the stratification
> algorithm does not work? I think you can just output the Plex as 
> binary or HDF5, but you could
> also give me code that makes one.

OK, I've attached an example dual porosity / MINC DM, output before the 
attempted redistribution.

  The original mesh is a vertical column of 11 hexahedral cells. Six of 
the cells (2 - 7) have two levels of interior MINC cells added (i.e. an 
interior cell, and then another interior cell inside that). So there are 
a total of 11 + 2*6 = 23 cells. The DM was distributed over two processes.

For each extra interior face I also create one interior edge and one 
interior vertex. I give the interior vertices arbitrary (0,0,0) 
locations. I assign all the extra interior DM points (cells, vertices, 
faces, edges) cell type DM_POLYTOPE_INTERIOR_GHOST.

- Adrian

-- 
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / University of Auckland, New Zealand
email:a.croucher at auckland.ac.nz
tel: +64 (0)9 923 4611
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20240307/510f7912/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: minc_dm.h5
Type: application/x-hdf
Size: 132888 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20240307/510f7912/attachment-0001.hdf>