<div dir="ltr"><div>Hi everyone,</div><div>this time I am using PETSc to do something that is more complicated than my usual and I want to do it at the highest possible abstraction level.</div><div><br></div><div>To put it in a nutshell, my intent is to build a parallel multigrid to solve a linear system via DM, KSP and PCMG (I would like to use DMMG but probably I should have the same problems or more).</div><div><br></div><div>I created the distributed object, <b>da</b>, with DMDACreate3d, even if it is distributed (as yet) only in the x-dimension and has 3 dof.</div><div>Then I create the KSP (type KSPRICHARDSON) and set the nonzero initial guess and PCMG as preconditioner. Here I start to tune the MG. </div><div><br></div><div>The point is that I need to define all the operators as matrix-free, since they will do several operations on <b>x</b> to obtain <b>y</b>, and I am not familiar with the way to access all the elements or informations in the two levels involved and/or among the processors with a so-high level interface.</div><div> </div><div>So please (please please please please) tell me if I correctly understand the mechanism or I am on the wrong way and clear my doubts.</div><div><br></div><div>That is, let's say that my operation for the shell are:</div><div><ul><li><b>A_mult(Mat mat,Vec x, Vec y) </b>//coefficients matrix</li></ul>in this case the level is only one but should I write it taking into account only the local data (I think so) and accessing them via the informations in <b>da</b>? </div><div><br></div><div>For example, if I use DMDAVecGetArray, DMDAVecGetCorners (or DMDAGetGhostCorners) and DMDAVecRestoreArray, will they retrieve informations from the right level each time (I am pretty sure that in some official examples it is done in this way)?</div><div><br></div><div>Or should I handle just <b>Vec</b>s as local structures with their usual indices (through VecGetArray and VecRestoreArray)?<ul><li><b>P_mult(Mat mat,Vec x, Vec y) </b>//interpolation matrix that is NOT conceptually the traspose of Restriction</li></ul>in this case <b>x</b> and <b>y</b> will be from two different levels (respectively L and L+1), so, if I retrieve informations from the <b>da</b>... how can I access the two at different levels?<br></div><div><br></div><div>I am sorry if it seems that they are trivial questions, and I will be grateful to anyone will help me. </div><div><br></div><div>Thanks a lot,</div><div>Valeria</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div style="font-size:12.8px"><font color="#666666">---------------------------------------------------------------------------------------------<span style="background-color:rgb(255,255,255)"><br></span></font></div><div style="font-size:12.8px"><font color="#666666"><span style="background-color:rgb(255,255,255)">PhD Valeria Mele</span><br><br>University of Naples Federico II</font></div><div style="font-size:12.8px"><font color="#666666">Department of Mathematics and Applications "R. Caccioppoli"<br></font></div><div style="font-size:12.8px"><font color="#666666">Complesso Universitario M.S. Angelo, Via Cinthia<br></font></div><div style="font-size:12.8px"><font color="#666666">80126 Naples </font></div><div style="font-size:12.8px"><font color="#666666">---------------------------------------------------------------------------------------------</font><font color="#000000"><br></font></div></div></div></div>
</div>