<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>