<html><head></head><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:16px"><div id="yiv5237025425"><div id="yui_3_16_0_1_1468641671620_32857"><div style="color:#000;background-color:#fff;font-family:times new roman, new york, times, serif;font-size:16px;" id="yui_3_16_0_1_1468641671620_32856"><div id="yiv5237025425yui_3_16_0_1_1468641671620_29819"><span id="yiv5237025425yui_3_16_0_1_1468641671620_29895">Thanks for all the pointers. </span></div><div id="yiv5237025425yui_3_16_0_1_1468641671620_29819"><span><br clear="none"></span></div><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29819"><span id="yiv5237025425yui_3_16_0_1_1468641671620_29990">I am happy to switch to </span><span id="yiv5237025425yui_3_16_0_1_1468641671620_29991" style="font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;">ksp/examples/tutorials/ex25.c in a first instance as you suggest.</span></div><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29819"><span style="font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;"><br clear="none"></span></div><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29819"><span id="yiv5237025425yui_3_16_0_1_1468641671620_30077" style="font-family:HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif;"> I am still stuck with the same issue as before though. I am trying to extract the hierarchy </span></div><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29819"> of coarser grid matrices and the intergrid transfer operators from the DMDA data structure. I would </div><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29819"> like to modify these operators and define a multigrid cycle with the modified operators. </div><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29819"><br clear="none"></div><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29819"> Given A^h (Helmholtz) and M^h (shifted Laplace), I would like to define a multigrid cycle involving </div><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29819"> both A^H and M^H. Can I rely on the multilevel DMDA structure to construct A^H and M^H for me </div><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29819"> in a set-up phase, plug them into a user-defined context, and plug them back out in a solve phase? </div><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29819"><br clear="none"></div><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29819">Thanks, Domenico. </div><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29819"><br clear="none"></div><div class="yiv5237025425qtdSeparateBR" id="yiv5237025425yui_3_16_0_1_1468641671620_29808"><br clear="none"></div><div class="yiv5237025425yqt8939338427" id="yiv5237025425yqt03652"></div></div></div></div><div class=".yiv5237025425yahoo_quoted" id="yui_3_16_0_1_1468641671620_32858"> <div id="yiv5237025425yui_3_16_0_1_1468641671620_29812" style="font-family:times new roman, new york, times, serif;font-size:16px;"> <div id="yiv5237025425yui_3_16_0_1_1468641671620_29811" style="font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:16px;"> <div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29810"> <font id="yiv5237025425yui_3_16_0_1_1468641671620_29809" size="2" face="Arial"> </font><hr id="yiv5237025425yui_3_16_0_1_1468641671620_29827" size="1"> <b id="yiv5237025425yui_3_16_0_1_1468641671620_29899"><span id="yiv5237025425yui_3_16_0_1_1468641671620_29898" style="font-weight:bold;">From:</span></b> Matthew Knepley <knepley@gmail.com><br clear="none"> <b id="yiv5237025425yui_3_16_0_1_1468641671620_29897"><span id="yiv5237025425yui_3_16_0_1_1468641671620_29896" style="font-weight:bold;">To:</span></b> Barry Smith <bsmith@mcs.anl.gov> <br clear="none"><b><span style="font-weight:bold;">Cc:</span></b> domenico lahaye <domenico_lahaye@yahoo.com>; "petsc-users@mcs.anl.gov" <petsc-users@mcs.anl.gov><br clear="none"> <b><span style="font-weight:bold;">Sent:</span></b> Sunday, July 17, 2016 2:29 PM<br clear="none"> <b id="yiv5237025425yui_3_16_0_1_1468641671620_29815"><span id="yiv5237025425yui_3_16_0_1_1468641671620_29814" style="font-weight:bold;">Subject:</span></b> Re: [petsc-users] Regarding ksp ex42 - Citations<br clear="none"> </div> <div class="yiv5237025425y_msg_container" id="yiv5237025425yui_3_16_0_1_1468641671620_29816"><br clear="none"><div id="yiv5237025425"><div id="yiv5237025425yui_3_16_0_1_1468641671620_29824"><div dir="ltr" id="yiv5237025425yui_3_16_0_1_1468641671620_29823"><div class="yiv5237025425gmail_extra" id="yiv5237025425yui_3_16_0_1_1468641671620_29822"><div class="yiv5237025425gmail_quote" id="yiv5237025425yui_3_16_0_1_1468641671620_29821">On Sat, Jul 16, 2016 at 10:11 PM, Barry Smith <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:bsmith@mcs.anl.gov" target="_blank" href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br clear="none"><blockquote class="yiv5237025425gmail_quote" id="yiv5237025425yui_3_16_0_1_1468641671620_29825" 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="yiv5237025425" id="yiv5237025425yui_3_16_0_1_1468641671620_29826"><br clear="none">
> On Jul 14, 2016, at 12:21 PM, domenico lahaye <<a rel="nofollow" shape="rect" ymailto="mailto:domenico_lahaye@yahoo.com" target="_blank" href="mailto:domenico_lahaye@yahoo.com">domenico_lahaye@yahoo.com</a>> wrote:<br clear="none">
><br clear="none">
</span><div id="yiv5237025425yui_3_16_0_1_1468641671620_29944"><div class="yiv5237025425h5" id="yiv5237025425yui_3_16_0_1_1468641671620_29943">> Dear PETSc team,<br clear="none">
><br clear="none">
> 1) I am looking into ks/examples/tutorials/ex42.c I am still new to the DMDA structure<br clear="none">
> and likely not giving it as much time as it deserves. However, I do not see immediately<br clear="none">
> what function is responsible for calling PCMGSetSmoother and PCMGSetResidual.<br clear="none">
><br clear="none">
> I tried to call PCMGGetCoarseSolve(pc, &kcpc) and subsequently<br clear="none">
> KSPGetOperators (kspc, ... ) to check how the coarse grid operator is defined<br clear="none">
> after calling DMCoarsenHierarchy, but that failed.<br clear="none">
><br clear="none">
> I am solving Helmholtz with shifted Laplace, and managed to exploit DMDA to perform<br clear="none">
> a multigrid solve on the preconditioner. In a next stage I want to implement the deflation<br clear="none">
> using DMDA as well.<br clear="none">
><br clear="none">
> 2) On <a rel="nofollow" shape="rect" target="_blank" href="http://www.mcs.anl.gov/petsc/documentation/referencing.html">http://www.mcs.anl.gov/petsc/documentation/referencing.html</a> I see<br clear="none">
><br clear="none">
> @Misc{petsc-web-page,<br clear="none">
> author = {Satish Balay and Shrirang Abhyankar and Mark~F. Adams and Jed Brown and Peter Brune<br clear="none">
> and Kris Buschelman and Lisandro Dalcin and Victor Eijkhout and William~D. Gropp<br clear="none">
> and Dinesh Kaushik and Matthew~G. Knepley<br clear="none">
> and Lois Curfman McInnes and Karl Rupp and Barry~F. Smith<br clear="none">
> and Stefano Zampini and Hong Zhang and Hong Zhang},<br clear="none">
> title = {{PETS}c {W}eb page},<br clear="none">
> url = {<a rel="nofollow" shape="rect" target="_blank" href="http://www.mcs.anl.gov/petsc">http://www.mcs.anl.gov/petsc</a>},<br clear="none">
> howpublished = {\url{<a rel="nofollow" shape="rect" target="_blank" href="http://www.mcs.anl.gov/petsc}">http://www.mcs.anl.gov/petsc}</a>},<br clear="none">
> year = {2016}<br clear="none">
> }<br clear="none">
><br clear="none">
><br clear="none">
><br clear="none">
> Is the last author mentioned twice intentionally?<br clear="none">
><br clear="none">
> 3) On <a rel="nofollow" shape="rect" target="_blank" href="http://www.mcs.anl.gov/petsc/publications/petscapps-bib.html#OpenFOAM%202.2.1">http://www.mcs.anl.gov/petsc/publications/petscapps-bib.html#OpenFOAM%202.2.1</a> I see<br clear="none">
><br clear="none">
> @misc{OpenFOAM<br clear="none">
> ,<br clear="none">
><br clear="none">
><br clear="none">
> title = "OpenFOAM",<br clear="none">
><br clear="none">
> howpublished = "\url{<a rel="nofollow" shape="rect" target="_blank" href="http://www.openfoam.com/">http://www.openfoam.com</a>}",<br clear="none">
><br clear="none">
> url = {<a rel="nofollow" shape="rect" target="_blank" href="http://www.openfoam.com/">http://www.openfoam.com</a>},<br clear="none">
><br clear="none">
> note = "OpenFOAM is a free, open source CFD software package. It allows PETSc linear algebra and solvers to be used underneath.",<br clear="none">
><br clear="none">
> key = "OpenFOAM 2.2.1"<br clear="none">
><br clear="none">
> }<br clear="none">
><br clear="none">
><br clear="none">
> Do you have more information on the use of PETSc within OpenFoam?<br clear="none">
<br clear="none">
</div></div> Very good question. It seems that this citation is wrong or no longer valid; I have removed it from the PETSc repository. I could find no mention of PETSc usage in the OpenFoam and its third party packages. I think we should not have been listing this citation.</blockquote><div><br clear="none"></div><div>This suggests that people are using it with OpenFOAM: <a rel="nofollow" shape="rect" target="_blank" href="http://powerlab.fsb.hr/ped/kturbo/OpenFOAM/slides/PatersonNuTTS2009.pdf">http://powerlab.fsb.hr/ped/kturbo/OpenFOAM/slides/PatersonNuTTS2009.pdf</a></div><div><br clear="none"></div><div>In fact, they use PETSc in the dynamic overset grid implementation for OpenFOAM, which I think is an approved extension:</div><div><br clear="none"></div><div> <a rel="nofollow" shape="rect" target="_blank" href="http://web.student.chalmers.se/groups/ofw5/Abstracts/DavidBogerAbstractOFW5.pdf">http://web.student.chalmers.se/groups/ofw5/Abstracts/DavidBogerAbstractOFW5.pdf</a></div><div class="yiv5237025425yqt3783790807" id="yiv5237025425yqtfd60844"><div><br clear="none"></div><div> Matt</div><div> </div></div><blockquote class="yiv5237025425gmail_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 class="yiv5237025425yqt3783790807" id="yiv5237025425yqtfd65670"><span class="yiv5237025425"><font color="#888888"><br clear="none">
Barry<br clear="none">
</font></span></div><div class="yiv5237025425"><div class="yiv5237025425h5"><div class="yiv5237025425yqt3783790807" id="yiv5237025425yqtfd67569"><br clear="none">
><br clear="none">
> 4) @matt in response to a question he raised in Vienna<br clear="none">
><br clear="none">
> MIPSE is a BEM solver. Details are on:<br clear="none">
> <a rel="nofollow" shape="rect" target="_blank" href="http://www.g2elab.grenoble-inp.fr/plateforms/mipse-modeling-of-interconnected-power-systems-632862.kjsp?RH=G2ELAB_R-MAGE">http://www.g2elab.grenoble-inp.fr/plateforms/mipse-modeling-of-interconnected-power-systems-632862.kjsp?RH=G2ELAB_R-MAGE</a><br clear="none">
><br clear="none">
> Cheers, Domenico Lahaye.</div><br clear="none">
><br clear="none">
<br clear="none">
</div></div></blockquote></div><br clear="none"><br clear="all"><div><br clear="none"></div>-- <br clear="none"><div class="yiv5237025425gmail_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 clear="none">-- Norbert Wiener</div><div class="yiv5237025425yqt3783790807" id="yiv5237025425yqtfd34108">
</div></div></div></div></div><br clear="none"><br clear="none"></div> </div> </div> </div></div></body></html>