<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"><DIV><BR> </DIV>
<DIV>Hello, Barry</DIV>
<DIV><BR>    Is it true that BJacobi preconditioner is only one type of ASM preconditioner(overlap=0), it is only one part of ASM?<BR>    I find that sometimes the two have the same results while sometimes they are not. Is the reason that PETSc can <BR>internally get the better parameter when using ASM, and sometimes it use overlap=0 then ASM and BJacobi have the same<BR>results.<BR>    If it is true, ASM is always better than BJacobi?</DIV>
<DIV>Thanks.                   Jim<BR><BR><BR></DIV>
<DIV></DIV>
<DIV id="divNeteaseMailCard"></DIV>
<DIV><BR></DIV><PRE><BR>At 2013-02-17 01:57:46,"Barry Smith" <bsmith@mcs.anl.gov> wrote:
>
>  This is completely possible. The convergence of iterative methods is a complicated business. Even something seemingly simple like block Jacobi whose conver> gence depends not only on the number of blocks, but the "shapes" of the blocks and what part of the matrix is "discarded" in the preconditioner depending on the > number and "shape" of the blocks. Decomposing the domain into 4 and 8 pieces results in different "shapes" of the blocks and different parts of the matrix being "> discarded" in the preconditioner.
>
>   Barry
>
>On Feb 16, 2013, at 11:36 AM, w_ang_temp <w_ang_temp@163.com> wrote:
>
>> 
>> 
>> I use the two commands in the same project. The first is divergent while the second is convergent.
>> nohup mpiexec -n 4 ./ex4f -ksp_type bcgs -pc_type bjacobi -ksp_rtol 1.0e-5 -ksp_converged_reason >out.txt &
>> nohup mpiexec -n 8 ./ex4f -ksp_type bcgs -pc_type bjacobi -ksp_rtol 1.0e-5 -ksp_converged_reason >out.txt &
>> So what is the reason?
>> Thanks.                    Jim
>> 
>> 
>> 
>> >> 2013-02-07 13:16:51£¬"Matthew Knepley" <knepley@gmail.com> Ð´µÀ£º
>> >> On Thu, Feb 7, 2013 at 12:11 AM, w_ang_temp <w_ang_temp@163.com> wrote:
>> >> Hello,
>> >>     I use the same project, but I find that when different number of processors is choosed,
>> >> the convergence is different. For example, when the processors are 4, it is divergent; when
>> >> the processors are 8, it is convergent.
>> >>     So what is the reason?
>> 
>> >It is likely that your preconditioner changed.
>> 
>>   > Matt
>>  
>> >>     Thanks.
>> >>                                    Jim.
>> 
>> 
>> 
>> 
>> 
>> -- 
>> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
>> -- Norbert Wiener
>> 
>> 
>
</PRE></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>