<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Aug 6, 2015 at 10:09 AM, Gautam Bisht <span dir="ltr"><<a href="mailto:gbisht@lbl.gov" target="_blank">gbisht@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Matt,<div><br></div><div>During the last week discussion with Mark, it became obvious that to add PETSc's AMR capability in PFLOTRAN would require incorporation of DMPlex. </div><div><br></div><div>Currently, PFLOTRAN has multiple `modes` which include:</div><div>- RICHARDS: Single phase mass conservation,</div><div>-TH: Single phase mass + energy conservation,</div><div>-GENERAL: Two phase mass + energy conservation, etc.</div><div><br></div><div>All the different modes have portions of duplicate code for Residual/Jacobian evaluation. I believe incorporation of DMComposite would possibly shrink the current PFLOTRAN code.</div><div><br></div><div>Thus, I was curious to check if an example for DMPlex+DMComposite exists.</div></div></blockquote><div><br></div><div>We have only used DMComposite to deal with things that have really regular discretizations and topologies, namely</div><div>that work for DMDA. DMComposite will not calculate the Jacobian layout (you can tell it), and its not trivial to calculate</div><div>the blocks yourself (in my opinion).</div><div><br></div><div>I tend to use Plex by itself for these kind of situations. I guess it depends on what you want to be in charge of. If you</div><div>want to run the whole show (cell iteration, local numbering, etc.) then DMComposite is a good choice. If you want to</div><div>give up some control, and only manage element-wise residuals/Jacobians then Plex may be a good choice.</div><div><br></div><div>Maybe I am missing something in how you want to set this up.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><span class="HOEnZb"><font color="#888888"><div>-Gautam.</div><div><br></div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 6, 2015 at 8:08 AM, Gautam Bisht <span dir="ltr"><<a href="mailto:gbisht@lbl.gov" target="_blank">gbisht@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra">I'm going to move this thread over to dev mailing list.</div><span><font color="#888888"><div class="gmail_extra"><br></div><div class="gmail_extra">-Gautam.</div></font></span><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 6, 2015 at 4:44 AM, 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span>On Wed, Aug 5, 2015 at 10:21 PM, Gautam Bisht <span dir="ltr"><<a href="mailto:gbisht@lbl.gov" target="_blank">gbisht@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Hi Matt,</div><div><br></div><div>Instead of using gcc4.9, I reinstalled PETSc using clang on mac os x 10.10 and the example runs fine.</div><div><br></div><div>Btw, are there any examples that use DMPlex+DMComposite?</div></div></blockquote><div><br></div></span><div>I don't think so. What would you anticipate using it for?</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Thanks,</div><div>-Gautam.</div></div>
</blockquote></div><span><br><br clear="all"><div><br></div>-- <br><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>
</span></div></div>
</blockquote></div><br></div></div></div></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>