Thanks, Dr. Smith,<div>   I guess it is wrong since I validate with my own code to solve a Poisson equation: Delta p = 1.0. The result from PETSc is exactly the negative value of what I get. The file is in /src/ksp/ksp/example/tutorial/ex29.c</div>
<div><br></div><div>best,</div><div>Alan</div><div><br><div class="gmail_quote">On Wed, Sep 21, 2011 at 12:03 PM, Barry Smith <span dir="ltr">&lt;<a href="mailto:bsmith@mcs.anl.gov">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;"><br>
  Alan,<br>
<br>
   It is very possible that the comment in the example code is wrong and has an incorrect sign. If you tell us what example this is we&#39;ll check it and fix the comment it if is wrong.<br>
<font color="#888888"><br>
   Barry<br>
</font><div><div></div><div class="h5"><br>
On Sep 21, 2011, at 11:28 AM, Alan Wei wrote:<br>
<br>
&gt; Oh, Thanks, Matt,<br>
&gt;    I got a little bit confused, since in the code, it described:<br>
&gt;          div \rho grad u = f,  0 &lt; x,y &lt; 1,<br>
&gt;    But you said, the solver solves -\Delta u = f  (Eq.1), which means:<br>
&gt;    for example, to solve a equation like Delta p = 1, I should put rhs = -1 = f in (Eq.1) in the code, therefore -\Delta u = -1, which, then, will give me a good result for Delta p = 1, is that true?<br>
&gt;<br>
&gt; thanks in advance,<br>
&gt; Alan<br>
&gt;<br>
&gt; On Wed, Sep 21, 2011 at 11:19 AM, Matthew Knepley &lt;<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>&gt; wrote:<br>
&gt; On Wed, Sep 21, 2011 at 4:16 PM, Alan Wei &lt;<a href="mailto:zhenglun.wei@gmail.com">zhenglun.wei@gmail.com</a>&gt; wrote:<br>
&gt; However, why signs for v[] in ComputeMatrix, which 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:<br>
&gt; 1)  u[j][i] = (u[j+1][i] + u[j-1][i] + u[j][i+1] + u[j][i-1] - rhs * dx*dy)/4<br>
&gt; or<br>
&gt; 2)  4u[j][i] - u[j+1][i] - u[j-1][i] - u[j][i+1] - u[j][i-1] + rhs * dx*dy = 0<br>
&gt;<br>
&gt; The Laplacian is a negative definite operator, so we usually solver -\Delta u = f since that<br>
&gt; is a positive definite problem.<br>
&gt;<br>
&gt;   Thanks,<br>
&gt;<br>
&gt;      Matt<br>
&gt;<br>
&gt; thanks,<br>
&gt; Alan<br>
&gt;<br>
&gt;<br>
&gt; On Wed, Sep 21, 2011 at 8:22 AM, Barry Smith &lt;<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>&gt; wrote:<br>
&gt;<br>
&gt; On Sep 20, 2011, at 10:47 PM, Alan Wei wrote:<br>
&gt;<br>
&gt; &gt; Dear Dr. Smith,<br>
&gt; &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; &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>
&gt;<br>
&gt;   We do not change the sign of the right hand side.<br>
&gt;<br>
&gt;   Barry<br>
&gt;<br>
&gt; &gt;<br>
&gt; &gt; thanks in advance,<br>
&gt; &gt; Alan<br>
&gt; &gt;<br>
&gt; &gt; On Mon, Sep 19, 2011 at 8:43 PM, Barry Smith &lt;<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>&gt; wrote:<br>
&gt; &gt;<br>
&gt; &gt; On Sep 19, 2011, at 6:25 PM, Alan Wei wrote:<br>
&gt; &gt;<br>
&gt; &gt; &gt; Dear folks,<br>
&gt; &gt; &gt;     I hope you guys are having a nice day.<br>
&gt; &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; &gt;<br>
&gt; &gt;   -ksp_rtol 1.e-10 for example<br>
&gt; &gt;<br>
&gt; &gt;   Run with -ksp_monitor_true_residual -ksp_converged_reason<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &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; &gt;<br>
&gt; &gt;    Hmm, multigrid should likely converge well. Are you sure you&#39;ve set the BC&#39;s correctly?<br>
&gt; &gt;<br>
&gt; &gt;   Barry<br>
&gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; thanks in advance,<br>
&gt; &gt; &gt; Alan<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
&gt; -- Norbert Wiener<br>
&gt;<br>
<br>
</div></div></blockquote></div><br></div>