On Wed, Sep 21, 2011 at 4:16 PM, Alan Wei <span dir="ltr">&lt;<a href="mailto:zhenglun.wei@gmail.com">zhenglun.wei@gmail.com</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;">
However, why signs for v[] in ComputeMatrix, which <span style="font-family:arial, sans-serif;font-size:13px;background-color:rgb(255, 255, 255)">contains the values of a row of the matrix. They all have a negative signs. Therefore, I got confused which equation does this program solve for:</span><div>

<font face="arial, sans-serif">1)  u[j][i] = (u[j+1][i] + u[j-1][i] + u[j][i+1] + u[j][i-1] - rhs * dx*dy)/4</font></div><div><font face="arial, sans-serif">or</font></div>
<div><font face="arial, sans-serif">2)  4u[j][i] - u[j+1][i] - u[j-1][i] - u[j][i+1] - u[j][i-1] + rhs * dx*dy = 0</font></div></blockquote><div><br></div><div>The Laplacian is a negative definite operator, so we usually solver -\Delta u = f since that</div>
<div>is a positive definite problem.</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><font face="arial, sans-serif">thanks,</font></div><div><font face="arial, sans-serif">Alan</font></div><div><font face="arial, sans-serif"><br>
</font><br><div class="gmail_quote">On Wed, Sep 21, 2011 at 8:22 AM, Barry Smith <span dir="ltr">&lt;<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><br>
On Sep 20, 2011, at 10:47 PM, Alan Wei wrote:<br>
<br>
&gt; Dear Dr. Smith,<br>
&gt;     I figure out this problem. Actually, I made my own RHS, but I did not multiply them by the volume (dx*dy).<br>
&gt;     However, I met another problem. All values calculated from this example are exactly opposite to values from my own code. I wonder if the RHS I input by ComputeRHS are the really RHS or -1.*RHS?<br>
<br>
</div>   We do not change the sign of the right hand side.<br>
<div><div></div><div><br>
   Barry<br>
<br>
&gt;<br>
&gt; thanks in advance,<br>
&gt; Alan<br>
&gt;<br>
&gt; On Mon, Sep 19, 2011 at 8:43 PM, Barry Smith &lt;<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>&gt; wrote:<br>
&gt;<br>
&gt; On Sep 19, 2011, at 6:25 PM, Alan Wei wrote:<br>
&gt;<br>
&gt; &gt; Dear folks,<br>
&gt; &gt;     I hope you guys are having a nice day.<br>
&gt; &gt;     I&#39;m reading the /src/ksp/ksp/examples/tutorials/ex29.c.html and wonder how to set up the convergence criteria?<br>
&gt;<br>
&gt;   -ksp_rtol 1.e-10 for example<br>
&gt;<br>
&gt;   Run with -ksp_monitor_true_residual -ksp_converged_reason<br>
&gt;<br>
&gt;<br>
&gt; &gt; Currently I use it as a poisson solver to solve a Poisson Equation with three direction Neumann BC&#39;s and one direction Diriechlet BC&#39;s. It seems very bad on it.<br>
&gt;<br>
&gt;    Hmm, multigrid should likely converge well. Are you sure you&#39;ve set the BC&#39;s correctly?<br>
&gt;<br>
&gt;   Barry<br>
&gt;<br>
&gt; &gt;<br>
&gt; &gt; thanks in advance,<br>
&gt; &gt; Alan<br>
&gt;<br>
&gt;<br>
<br>
</div></div></blockquote></div><br>
</div>
</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>