<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"><DIV>Hello<BR>     As I become familiar with PETSc, I make a conclusion about the question and hope that it is right.<BR>     In order to get satisfactory results when solving a linear system with KSP, the Krylov space methods</DIV>
<DIV>are typically used in conjunction with a preconditioner. So the user should  select the effective Krylov</DIV>
<DIV>method(such as cg) and preconditioner(such as jacobi) for this linear system Ax=b. Or it maybe not</DIV>
<DIV>convergence.<BR>     Thanks.<BR>                              Jim<BR><BR></DIV>
<DIV></DIV>
<DIV id="divNeteaseMailCard"></DIV>
<DIV><BR></DIV>>At 2012-05-26 23:05:26,w_ang_temp <w_ang_temp@163.com> wrote:<BR>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" id="isReplyContent">
<DIV style="LINE-HEIGHT: 1.7; FONT-FAMILY: arial; COLOR: #000000; FONT-SIZE: 14px">
<DIV>>Hello, Matt</DIV>
<DIV>>     The results of the parallel and the serial  are <SPAN id="src_1_3" data-aligning="#src_1_3,#tran_1_4" jQuery17103119459660716153="97">identical by the "-ksp_rtol 1.0e-15 -ksp_converged_reason -ksp_monitor_true_residual".</SPAN></DIV>
<DIV><SPAN data-aligning="#src_1_3,#tran_1_4" jQuery17103119459660716153="97">>     I will make a <SPAN id="tran_2_8" data-aligning="#src_2_8,#tran_2_8" jQuery17107117125654621121="29">careful</SPAN> <SPAN id="tran_2_9" data-aligning="#src_2_9,#tran_2_9" jQuery17107117125654621121="30">analysis</SPAN> of these  <SPAN id="tran_1_6" data-aligning="#src_1_6,#tran_1_6" jQuery17107117125654621121="64">parameters next day because it is late night here now.</SPAN></SPAN></DIV>
<DIV><SPAN data-aligning="#src_1_3,#tran_1_4" jQuery17103119459660716153="97"><SPAN data-aligning="#src_1_6,#tran_1_6" jQuery17107117125654621121="64">>     Thanks again.                                                       Jim</SPAN></SPAN><BR></DIV>>ÔÚ 2012-05-26 22:26:20£¬"Matthew Knepley" <<A href="mailto:knepley@gmail.com">knepley@gmail.com</A>> Ð´µÀ£º<BR>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" id="isReplyContent">>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 style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<DIV style="LINE-HEIGHT: 1.7; FONT-FAMILY: arial; FONT-SIZE: 14px">
<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></BLOCKQUOTE></DIV></BLOCKQUOTE></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>