<!DOCTYPE html><!-- BaNnErBlUrFlE-BoDy-start -->
<!-- Preheader Text : BEGIN -->
<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">
hi Matt, On 20/01/24 6: 03 am, Matthew Knepley wrote: I put in a test for redistribution with overlap, Plex test ex12, test 7. You can run it using make -f ./gmakefile test search="dm_impls_plex_tests-ex12_7" and see the output in
</div>
<!-- Preheader Text : END -->
<!-- Email Banner : BEGIN -->
<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">ZjQcmQRYFpfptBannerStart</div>
<!--[if ((ie)|(mso))]>
<table border="0" cellspacing="0" cellpadding="0" width="100%" style="padding: 16px 0px 16px 0px; direction: ltr" ><tr><td>
<table border="0" cellspacing="0" cellpadding="0" style="padding: 0px 10px 5px 6px; width: 100%; border-radius:4px; border-top:4px solid #90a4ae;background-color:#D0D8DC;"><tr><td valign="top">
<table align="left" border="0" cellspacing="0" cellpadding="0" style="padding: 4px 8px 4px 8px">
<tr><td style="color:#000000; font-family: 'Arial', sans-serif; font-weight:bold; font-size:14px; direction: ltr">
This Message Is From an External Sender
</td></tr>
<tr><td style="color:#000000; font-weight:normal; font-family: 'Arial', sans-serif; font-size:12px; direction: ltr">
This message came from outside your organization.
</td></tr>
</table>
</td></tr></table>
</td></tr></table>
<![endif]-->
<![if !((ie)|(mso))]>
<div dir="ltr" id="pfptBanneroumnxqr" style="all: revert !important; display:block !important; text-align: left !important; margin:16px 0px 16px 0px !important; padding:8px 16px 8px 16px !important; border-radius: 4px !important; min-width: 200px !important; background-color: #D0D8DC !important; background-color: #D0D8DC; border-top: 4px solid #90a4ae !important; border-top: 4px solid #90a4ae;">
<div id="pfptBanneroumnxqr" style="all: unset !important; float:left !important; display:block !important; margin: 0px 0px 1px 0px !important; max-width: 600px !important;">
<div id="pfptBanneroumnxqr" style="all: unset !important; display:block !important; visibility: visible !important; background-color: #D0D8DC !important; color:#000000 !important; color:#000000; font-family: 'Arial', sans-serif !important; font-family: 'Arial', sans-serif; font-weight:bold !important; font-weight:bold; font-size:14px !important; line-height:18px !important; line-height:18px">
This Message Is From an External Sender
</div>
<div id="pfptBanneroumnxqr" style="all: unset !important; display:block !important; visibility: visible !important; background-color: #D0D8DC !important; color:#000000 !important; color:#000000; font-weight:normal; font-family: 'Arial', sans-serif !important; font-family: 'Arial', sans-serif; font-size:12px !important; line-height:18px !important; line-height:18px; margin-top:2px !important;">
This message came from outside your organization.
</div>
</div>
<div style="clear: both !important; display: block !important; visibility: hidden !important; line-height: 0 !important; font-size: 0.01px !important; height: 0px"> </div>
</div>
<![endif]>
<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">ZjQcmQRYFpfptBannerEnd</div>
<!-- Email Banner : END -->
<!-- BaNnErBlUrFlE-BoDy-end -->
<html><head><!-- BaNnErBlUrFlE-HeAdEr-start -->
<style>
#pfptBanneroumnxqr { all: revert !important; display: block !important;
visibility: visible !important; opacity: 1 !important;
background-color: #D0D8DC !important;
max-width: none !important; max-height: none !important }
.pfptPrimaryButtonoumnxqr:hover, .pfptPrimaryButtonoumnxqr:focus {
background-color: #b4c1c7 !important; }
.pfptPrimaryButtonoumnxqr:active {
background-color: #90a4ae !important; }
</style>
<!-- BaNnErBlUrFlE-HeAdEr-end -->
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body>
<p>hi Matt,<br>
</p>
<div class="moz-cite-prefix">On 20/01/24 6:03 am, Matthew Knepley
wrote:<br>
</div>
<blockquote type="cite" cite="mid:CAMYG4G=GziHyCWzKurY-2NCSnbSEoo57Z-_D9pCOGEh9QND3Xg@mail.gmail.com">
<div dir="ltr">
<div class="gmail_quote"><br>
<div>I put in a test for redistribution with overlap, Plex
test ex12, test 7. You can run it using</div>
<div><br>
</div>
<div> make -f ./gmakefile test
search="dm_impls_plex_tests-ex12_7"</div>
<div><br>
</div>
<div>and see the output in
src/dm/impls/plex/tests/output/ex12_7.out. The second
DMPlexDistribute() call</div>
<div>is made on a mesh with overlap.</div>
<div><br>
</div>
<div>If this is not a solution, can we make a small test that
goes in ex12 that shows what is wrong?</div>
</div>
</div>
</blockquote>
<p>I just got around to testing this with PETSc 3.20.5.</p>
<p>I am still getting problems with redistribution of MINC (dual
porosity) meshes. The function DMPlexCheckPointSF() is now raising
errors like this:</p>
<p>[0]PETSC ERROR: Petsc has generated inconsistent data<br>
[0]PETSC ERROR: Point SF contains 25 which is a cell<br>
</p>
<p>Looking in the PETSc code it appears some extra checks were added
about a year ago (commit 7726db96) including the bit commented
"Check there are no cells in interface", which is what is raising
the error.</p>
<p>My understanding of the point SF is a bit rusty now but I would
have thought that for finite volume meshes with overlap 1 (which
is what I have), there should be cells in the interface?</p>
<p>Cheers, Adrian<br>
</p>
<pre class="moz-signature" cols="72">--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / 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>