<html><head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body>
<p>hi Matt,<br>
</p>
<div class="moz-cite-prefix">On 28/03/23 00:14, Matthew Knepley
wrote:<br>
</div>
<blockquote type="cite" cite="mid:CAMYG4GnBgP4_+ypTxJVxvhWngf9geWuJ54MdRyeDcrFuXW80RQ@mail.gmail.com">
<div dir="ltr">
<div class="gmail_quote">Yes, it is a new check. I do not expect
cells that are shared until we create an overlap. We
<div>can just put in a flag to shut it off, since otherwise it
is a great check for finding bugs. Can</div>
<div>you show me how you add cells, so we can make a small
example that tests this to go in</div>
<div>PETSc?</div>
</div>
</div>
</blockquote>
<p>Essentially I'm adding (one or more) cells nested "inside" some
of the original DM cells to simulate dual porosity. As well as new
cells I also create a new face between the original cell and the
new one, and a new edge and point for the face (so the DAG is the
same depth everywhere - otherwise it causes trouble).<br>
</p>
<p>In the code, I create a new DMPlex, figure out the chart
including the new cells etc., set the cone sizes and the cones. I
symmetrize it, copy labels over from the original DM, set cell
types (using DM_POLYTOPE_INTERIOR_GHOST for the new points) and
set up the point SF for the new DM. If something is wrong with
that point SF, would it cause the error I'm seeing when I
redistribute?</p>
<p>The code for setting up the new point SF is here (hope this link
works):<br>
</p>
<p><a class="moz-txt-link-freetext" href="https://github.com/waiwera/waiwera/blob/55d26483b3392b70985fb25f998fbda3297c6154/src/mesh.F90#L3408">https://github.com/waiwera/waiwera/blob/55d26483b3392b70985fb25f998fbda3297c6154/src/mesh.F90#L3408</a></p>
<p>- Adrian<br>
</p>
<pre class="moz-signature" cols="72">--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: <a class="moz-txt-link-abbreviated" href="mailto:a.croucher@auckland.ac.nz">a.croucher@auckland.ac.nz</a>
tel: +64 (0)9 923 4611
</pre>
</body>
</html>