On Sat, May 26, 2012 at 2:22 PM, w_ang_temp <span dir="ltr"><<a href="mailto:w_ang_temp@163.com" target="_blank">w_ang_temp@163.com</a>></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">
<div style="line-height:1.7;font-size:14px;font-family:arial"><div>Hello,Matt<br>    Thanks again.<br>    What you said("Use -ksp_type preonly -pc_type lu") means that I should use a direct solver.And I am puzzled.<br>
</div></div></blockquote><div><br></div><div>Okay, you want to remove variables from the investigation. You say there is a difference between</div><div>the parallel and serial solve. We do this all the time, so it is not a bug. It is an issue of understanding.</div>
<div>First, take away the tolerance issues:</div><div><br></div>  -ksp_rtol 1.0e-15 -ksp_converged_reason -ksp_monitor_true_residual</div><div class="gmail_quote"><br><div>and send output of the serial and parallel run.</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 style="line-height:1.7;font-size:14px;font-family:arial"><div>    (1)I think this problem (single processor is ok and more than one are problematic)is typical,but I cannot <br>
find useful infomation both in the manual and the petsc-users.So can you give me some explanation for this <br>phenomenon.<br>    (2)When I use -pc_type lu,the error message is "Matrix format mpiaij does not have a built-in PETSc LU!".<br>
    (3)The manual says that "This approach prevents the user from the unexpected surprise of having a <br>corrupted matrix after a linear solve"(P72).So do you mean the problem of mine is here?<br>    (4)Generally I want to use the iterative approach other than the direct method.<br>
    As a beginne
 r,I have to learn more about the manual.And also I think an experienced senior will be a more<br> useful help for me because I have gotten stuck on this problem for more than a week.<br>     Thanks.                                Jim<br>
<br></div>
<div></div>
<div></div>
<div><br></div>>At 2012-05-26 20:17:42,"Matthew Knepley" <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> wrote:<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex">>Use -ksp_type preonly -pc_type lu. 
<div><br></div>
<div> >  Matt</div></blockquote></div><br><br><span title="neteasefooter"><span></span></span></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>