<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><br class=""><div><br class=""><blockquote type="cite" class=""><div class="">5. 12. 2017 v 12:29, Matthew Knepley <<a href="mailto:knepley@gmail.com" class="">knepley@gmail.com</a>>:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" style="font-family: Menlo-Regular; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><div class="gmail_extra"><div class="gmail_quote">On Mon, Dec 4, 2017 at 1:06 PM, Matthew Knepley<span class="Apple-converted-space"> </span><span dir="ltr" class=""><<a href="mailto:knepley@gmail.com" target="_blank" class="">knepley@gmail.com</a>></span><span class="Apple-converted-space"> </span>wrote:<br class=""><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;"><div dir="ltr" class=""><div class="gmail_extra"><div class="gmail_quote"><span class="gmail-">On Fri, Dec 1, 2017 at 11:10 AM, Vaclav Hapla<span class="Apple-converted-space"> </span><span dir="ltr" class=""><<a href="mailto:vaclav.hapla@erdw.ethz.ch" target="_blank" class="">vaclav.hapla@erdw.ethz.ch</a>></span><span class="Apple-converted-space"> </span>wrote:<br class=""><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;">Hello<br class=""><br class="">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 class=""><br class=""> if (interpolate) {<br class=""> DM idm = NULL;<br class=""><br class=""> ierr = DMPlexInterpolate(*dm, &idm);CHKERRQ(ierr);<br class=""> ierr = DMPlexCopyCoordinates(*dm, idm);CHKERRQ(ierr);<br class=""> ierr = DMCopyLabels(*dm, idm);CHKERRQ(ierr);<br class=""> ierr = DMDestroy(dm);CHKERRQ(ierr);<br class=""> *dm = idm;<br class=""> }<br class=""><br class="">in DMPlexCreateFromFile after the DMLoad call.<br class=""></blockquote><div class=""><br class=""></div></span><div class="">Yes.</div><span class="gmail-"><div class=""> </div><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: 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 class=""></blockquote><div class=""><br class=""></div></span><div class="">Yes.</div><span class="gmail-"><div class=""> </div><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: 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 class=""><br class=""></div></span><div class="">Yes. My original motivation turned out to be not important. I will get rid of it.</div><div class=""><br class=""></div><div class="">I will make all changes in a branch and push to next.</div></div></div></div></blockquote><div class=""><br class=""></div><div class="">It is here.</div><div class=""><br class=""></div><div class=""> <a href="https://bitbucket.org/petsc/petsc/branch/knepley/fix-plex-interpolation-consistency" class="">https://bitbucket.org/petsc/petsc/branch/knepley/fix-plex-interpolation-consistency</a></div><div class=""><br class=""></div><div class=""> Matt</div></div></div></div></div></blockquote><div><br class=""></div><div>It looks good, thanks!</div><div><br class=""></div><div>Vaclav</div><br class=""><blockquote type="cite" class=""><div class=""><div dir="ltr" style="font-family: Menlo-Regular; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><div class="gmail_extra"><div class="gmail_quote"><div class=""> </div><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;"><div dir="ltr" class=""><div class="gmail_extra"><div class="gmail_quote"><div class=""> Thanks,</div><div class=""><br class=""></div><div class=""> Matt</div><div class=""> </div><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;"><span class="gmail-m_-1607746285838136088HOEnZb"><font color="#888888" class=""><br class="">Vaclav</font></span></blockquote></div><span class="gmail-HOEnZb"><font color="#888888" class="">--<span class="Apple-converted-space"> </span><br class=""><div class="gmail-m_-1607746285838136088gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class="">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br class="">-- Norbert Wiener</div><div class=""><br class=""></div><div class=""><a href="http://www.caam.rice.edu/~mk51/" target="_blank" class="">https://www.cse.buffalo.edu/~<wbr class="">knepley/</a><br class=""></div></div></div></div></div></font></span></div></div></blockquote></div><br class=""><br clear="all" class=""><div class=""><br class=""></div>--<span class="Apple-converted-space"> </span><br class=""><div class="gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class="">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br class="">-- Norbert Wiener</div><div class=""><br class=""></div><div class=""><a href="http://www.caam.rice.edu/~mk51/" target="_blank" class="">https://www.cse.buffalo.edu/~knepley/</a></div></div></div></div></div></div></div></div></blockquote></div><br class=""></body></html>