<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div><blockquote type="cite" class=""><div class=""><div dir="ltr" class=""><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word" class=""><div class="">Other suggestions on how to best integrate staggered finite differences within the current PETSc framework are ofcourse also highly welcome. </div><div class="">Our current thinking was to pack it into a DMSHELL (which has the problem of not having a restriction interface). </div><div class=""><br class=""></div></div></blockquote><div class=""><br class=""></div><div class="">Using DMShell is the cleanest approach.<br class=""><br class=""></div><div class="">An alternative is to have you user code simply take control of all of the configuration of the PCMG object. E.g. you call your user code which creates the restriction operator, you pull out the PC and call PCMGSetRestriction() on etc. This can be done easily performed in the context of linear problems. For non-linear problems, you could jam this setup code inside your ComputeJacobian function.<br class=""><br class=""></div><div class="">This is all possible, albeit clunky and kinda ugly. It works though if you need something before Barry adds the required support in PCMG.<br class=""></div><div class=""><br class=""></div></div></div></div></div></blockquote><div>that is indeed what we currently do, which does work. </div><div>Yet, as you say, it is clunky and does not allow setting up things like FAS in an easy manner, or add new physics to the code in a straightforward manner.</div><div>Having a cleaner interface would be really nice.</div><div><br class=""></div><div>Boris</div><div><br class=""></div><blockquote type="cite" class=""><div class=""><div dir="ltr" class=""><div class="gmail_extra"><div class="gmail_quote"><div class="">Cheers<br class=""></div><div class=""> Dave<br class=""></div><div class=""><br class=""><br class=""></div></div></div></div>
</div></blockquote></div><br class=""></body></html>