<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<p>hi Matt,<br>
</p>
<div class="moz-cite-prefix">On 29/06/19 12:07 AM, Matthew Knepley
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:CAMYG4GnPKvv5zb+EHKMpe6U+2RMYgi4Hddp4k4b=n=r-GsUWjQ@mail.gmail.com">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<div dir="ltr">
<div class="gmail_quote"><br>
<div>Okay, I ran through this. I am attaching my C version
which was easier for me to play with in the debugger, but it
should not be</div>
<div>hard to put the few lines into yours. The problem comes
from using overlap in the redistribution. You can see this
by running</div>
<div><br>
</div>
<div> mpiexec -n 2 ./distgeom2 -overlap 1 -overlap2 0</div>
<div><br>
</div>
<div>which works fine. If overlap2 is 1, then the face data is
shifted.</div>
</div>
</div>
</blockquote>
<p><br>
</p>
<p>Yes, that's the same thing I found.<br>
</p>
<blockquote type="cite"
cite="mid:CAMYG4GnPKvv5zb+EHKMpe6U+2RMYgi4Hddp4k4b=n=r-GsUWjQ@mail.gmail.com">
<div dir="ltr">
<div class="gmail_quote">
<div>Something in DistributeField() does not know that we
changed</div>
<div>the topology after we distributed by adding overlap
cells. This is strange since DMPlexMigrate() works on the
coordinates for</div>
<div>the overlapped mesh, however they always start at 0 in
the vector, so it must be that some offset is missing. I am
looking for it now.</div>
</div>
</div>
</blockquote>
<p><br>
</p>
<p>It looked to me like the problem might lie not in
DMPlexDistributeField() but in DMPlexDistribute(), because if you
view the redistribution SF coming out of the second call to
DMPlexDistribute() it looks wrong.</p>
<p>DMPlexDistributeField() appeared to be doing what it was told,
but the SF it was given didn't look right to me.</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>