<div dir="ltr"><div dir="ltr">On Wed, Jun 30, 2021 at 6:13 AM Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">On Tue, Jun 29, 2021 at 9:04 PM Adrian Croucher <<a href="mailto:a.croucher@auckland.ac.nz" target="_blank">a.croucher@auckland.ac.nz</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">hi<br>
<br>
I've been testing building my code with PETSc 3.15 and have hit a weird <br>
problem.<br>
<br>
If I create a DMPlex and add boundary ghost cells using <br>
DMPlexConstructGhostCells(), I expect the boundary ghost cells to have <br>
cell type DM_POLYTOPE_FV_GHOST, if I check the cell types using <br>
DMPlexGetCellTypeLabel() and DMLabelGetValue().<br>
<br>
That is true if I build using PETSc 3.14. But with 3.15 it appears the <br>
boundary ghost cells are ending up with cell type <br>
DM_POLYTOPE_INTERIOR_GHOST.<br>
<br>
These types are defined in include/petscdmtypes.h and currently <br>
DM_POLYTOPE_FV_GHOST = 11 and DM_POLYTOPE_INTERIOR_GHOST = 12.<br>
<br>
I've attached a minimal example program in Fortran. It creates a box <br>
mesh and adds boundary ghost cells, then counts how many have type <br>
DM_POLYTOPE_FV_GHOST.<br>
<br>
With PETSc 3.14 it is 24, the same as the number of ghost cells. With <br>
PETSc 3.15 it is zero, and the boundary ghost cells all have type 12.<br>
<br>
Any clues?<br></blockquote><div><br></div><div>Yes, src/dm/f90-mod/petscdm.h is missing PYRAMID. I will make a fix right now.</div></div></div></blockquote><div><br></div><div>Here it is: <a href="https://gitlab.com/petsc/petsc/-/merge_requests/4140">https://gitlab.com/petsc/petsc/-/merge_requests/4140</a></div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_quote"><div>  Thanks,</div><div><br></div><div>      Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Regards, Adrian<br>
<br>
-- <br>
Dr Adrian Croucher<br>
Senior Research Fellow<br>
Department of Engineering Science<br>
University of Auckland, New Zealand<br>
email: <a href="mailto:a.croucher@auckland.ac.nz" target="_blank">a.croucher@auckland.ac.nz</a><br>
tel: +64 (0)9 923 4611<br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>