<!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, I've made a bit of progress with this, using a debugger. I wasn't setting the overlap for my dual-porosity DM, so it was defaulting to zero, and when it got to the DMPlexCheckPointSF() call it was doing the check for cells in the overlap
</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="pfptBanner1wj7rdo" 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="pfptBanner1wj7rdo" style="all: unset !important; float:left !important; display:block !important; margin: 0px 0px 1px 0px !important; max-width: 600px !important;">
<div id="pfptBanner1wj7rdo" 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="pfptBanner1wj7rdo" 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>
#pfptBanner1wj7rdo { all: revert !important; display: block !important;
visibility: visible !important; opacity: 1 !important;
background-color: #D0D8DC !important;
max-width: none !important; max-height: none !important }
.pfptPrimaryButton1wj7rdo:hover, .pfptPrimaryButton1wj7rdo:focus {
background-color: #b4c1c7 !important; }
.pfptPrimaryButton1wj7rdo:active {
background-color: #90a4ae !important; }
</style>
<!-- BaNnErBlUrFlE-HeAdEr-end -->
<meta charset="UTF-8"></head><body><pre style="font-family: sans-serif; font-size: 100%; white-space: pre-wrap; word-wrap: break-word">hi Matt,
I've made a bit of progress with this, using a debugger.
I wasn't setting the overlap for my dual-porosity DM, so it was
defaulting to zero, and when it got to the DMPlexCheckPointSF() call it
was doing the check for cells in the overlap when it shouldn't have been
(this is only triggered "if (!overlap)").
Now I do a DMPlexSetOverlap() to make sure the dual-porosity DM has the
same overlap as the original DM. That gets rid of the error raised in
DMPlexCheckPointSF(), as the overlap is now 1.
However it is still falling over further on, when it tries to stratify
the redistributed DM:
[0]PETSC ERROR: Petsc has generated inconsistent data
[0]PETSC ERROR: New depth 2 range [0,141) overlaps with depth 1 range
[0,141)
[0]PETSC ERROR: #1 DMPlexCreateDepthStratum() at
/home/acro018/software/PETSc/code/src/dm/impls/plex/plex.c:4087
[0]PETSC ERROR: #2 DMPlexStratify() at
/home/acro018/software/PETSc/code/src/dm/impls/plex/plex.c:4211
[0]PETSC ERROR: #3 DMPlexDistributeCones() at
/home/acro018/software/PETSc/code/src/dm/impls/plex/plexdistribute.c:1164
[0]PETSC ERROR: #4 DMPlexMigrate() at
/home/acro018/software/PETSc/code/src/dm/impls/plex/plexdistribute.c:1659
[0]PETSC ERROR: #5 DMPlexDistribute() at
/home/acro018/software/PETSc/code/src/dm/impls/plex/plexdistribute.c:1778
I have checked the depth label on my dual-porosity DM for a simple
example in which this error is raised, and the label looks correct to me.
However I have a memory of trying DMPlexStratify() on my dual-porosity
DMs in the past and it not liking them. I create the depth label by
copying its values from the original DM and then explicitly assigning
depth values for the new dual-porosity cells, vertices, faces & edges.
When the DM is redistributed I'd have thought it shouldn't be necessary
(at least in my case) to stratify as the depth label is already there
and just needs to be redistributed as well? But I guess that might not
always be the case.
Has anything changed in this part of the PETSc code recently that might
make this stop working?
- Adrian
--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
Waipapa Taumata Rau / University of Auckland, New Zealand
email: a.croucher@auckland.ac.nz
tel: +64 (0)9 923 4611
</pre></body></html>