On Sun, Sep 4, 2011 at 10:08 AM, Алексей Рязанов <span dir="ltr">&lt;<a href="mailto:ram@ibrae.ac.ru">ram@ibrae.ac.ru</a>&gt;</span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<font size="2" face="arial, helvetica, sans-serif">Hello!</font><div><font size="2" face="arial, helvetica, sans-serif">Help me please to solve my tiny problem with dmmg.</font></div>

<div><font size="2" face="arial, helvetica, sans-serif"><br></font></div><div><font size="2" face="arial, helvetica, sans-serif">In the example src/ksp/ksp/examples/tutorials/ex22.c, which is discussed on the 111th page of users manual, we use function </font></div>


<div><font size="2" face="&#39;courier new&#39;, monospace">(line 36): DACreate3d(PETSC_COMM_WORLD,DA_NONPERIODIC,DA_STENCIL_STAR,-3,-3,-3,PETSC_DECIDE,PETSC_DECIDE,PETSC_DECIDE,1,1,0,0,0,&amp;da);</font></div>

<div><font size="2" face="arial, helvetica, sans-serif">to &quot;<i>create the DA that stores information about the coarsest grid you wish to use</i>&quot;.</font></div><div><font size="2" face="arial, helvetica, sans-serif">We can see that negative dimension values (<span style="background-color:rgb(255, 255, 255)">-M, -N, -P) </span>are used <span style="background-color:rgb(255, 255, 255)"> &quot;<i>to indicate that it may be set to a different value from the command line with -da_grid_x &lt;M&gt; -da_grid_y &lt;N&gt; -da_grid_z &lt;P&gt;&quot;. </i></span></font></div>


<div><font size="2" face="arial, helvetica, sans-serif">The question is: when I use -da_grid_ M,N,P options as the size of the coarsest matrix, what&#39;s the size of the finest (most accurate, original problem) matrix?</font></div>


<div><font size="2" face="arial, helvetica, sans-serif"><br></font></div><div><font size="2" face="arial, helvetica, sans-serif">I&#39;ve tried to clear this up, using VecVew on the residual vector and figured out, that when i&#39;m using M,N,P as input parameters (so define the coarsest grid as field of M*N*P dots), I get the residual vector of si<font style="line-height:23px">ze (4*M-3)*(4*N-3)*(4*P-3).</font></font></div>


<div><font style="line-height:23px" size="2" face="arial, helvetica, sans-serif">Is my result correct?  If yes, could you please explain why it happens in this way.</font></div></blockquote><div><br></div><div>When refinement happens, it adds a vertex between any two vertices, starting with the mesh you provide.</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><font style="line-height:23px" size="2" face="arial, helvetica, sans-serif">Sorry if my question is incorrect. I could mix up everything at all.</font></div>
<div><font style="line-height:23px" size="2" face="arial, helvetica, sans-serif">Thanks in advance!</font></div>

<div><font face="&#39;Times New Roman&#39;" size="4" style="line-height:23px"><br></font></div><font color="#888888"><div>-- <br><span style="border-collapse:collapse;font-family:arial, sans-serif;font-size:13px"><div>

<span style="border-collapse:collapse;font-family:arial, sans-serif;font-size:13px"><div>Best regards,</div><div>Alexey Ryazanov</div><div>______________________________________</div><div><span style="color:rgb(51, 51, 51);font-family:Arial, Helvetica, sans-serif"><font size="1">Nuclear Safety Institute of Russian Academy of Sciences </font></span></div>


</span></div><div><br></div></span><br>
</div>
</font></blockquote></div><br><br clear="all"><div><br></div>-- <br>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<br>