<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Dec 4, 2017 at 1:06 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class="gmail-">On Fri, Dec 1, 2017 at 11:10 AM, Vaclav Hapla <span dir="ltr"><<a href="mailto:vaclav.hapla@erdw.ethz.ch" target="_blank">vaclav.hapla@erdw.ethz.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello<br>
<br>
I noticed DMPlexCreateFromFile ignores the interpolate flag for HDF5. It's the only format for which PetscViewer API is used for loading, and this API has no means to specify interpolation. I guess there should be at least<br>
<br>
    if (interpolate) {<br>
      DM idm = NULL;<br>
<br>
      ierr = DMPlexInterpolate(*dm, &idm);CHKERRQ(ierr);<br>
      ierr = DMPlexCopyCoordinates(*dm, idm);CHKERRQ(ierr);<br>
      ierr = DMCopyLabels(*dm, idm);CHKERRQ(ierr);<br>
      ierr = DMDestroy(dm);CHKERRQ(ierr);<br>
      *dm = idm;<br>
    }<br>
<br>
in DMPlexCreateFromFile after the DMLoad call.<br></blockquote><div><br></div></span><div>Yes.</div><span class="gmail-"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
This brings me to more general thing to discuss. The behavior of different mesh constructors for interpolate=PETSC_TRUE is not unified. The most elaborate behavior is in DMPlexCreateDoublet which after DMPlexInterpolate also copies the DM name, coordinates and labels. So I suggest all this stuff is done in DMPlexInterpolate.<br></blockquote><div><br></div></span><div>Yes.</div><span class="gmail-"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
BTW I think it's kind of weird the output DM passed to DMPlexInterpolate must me nullified - cf. MatReuse.</blockquote><div><br></div></span><div>Yes. My original motivation turned out to be not important. I will get rid of it.</div><div><br></div><div>I will make all changes in a branch and push to next.</div></div></div></div></blockquote><div><br></div><div>It is here.</div><div><br></div><div>  <a href="https://bitbucket.org/petsc/petsc/branch/knepley/fix-plex-interpolation-consistency">https://bitbucket.org/petsc/petsc/branch/knepley/fix-plex-interpolation-consistency</a></div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span class="gmail-m_-1607746285838136088HOEnZb"><font color="#888888"><br>
Vaclav</font></span></blockquote></div><span class="gmail-HOEnZb"><font color="#888888">-- <br><div class="gmail-m_-1607746285838136088gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~<wbr>knepley/</a><br></div></div></div></div></div>
</font></span></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>