On Jan 16, 2013, at 7:55 AM, w_ang_temp  <w_ang_temp at 163.com> wrote:

> At 2013-01-16 03:02:32,"Barry Smith" <bsmith at mcs.anl.gov> wrote:
> >
> >> 1.701448294063e+04  >   1.e4*1.145582415879e+00  hence it declares 
> >> divergence
>   Hello, Barry
>   I made some tests and it is true.  Thanks.
>  
>   But in the mannual, both version 3.2 and 3.3, the rule 
> is:||rk||>dtol*||b||. It is ||b||, not r_0.
> My misunderstanding? Or the error in the manual?

   Error in the manual. Depending on the solver it is either the preconditioned 
norm or the unpreconditioned norm

>   Besides, in the mannual, both version 3.2 and 3.3, the default dtol=1.0E+5. 
> But from the results
> in the example, it is 1.0E+4.

    Another error in the manual

> I donot know the reason. My petsc version is 3.2-p7.

   None of this should really matter. The solver has gone bad at this point 
regardless of specifics of how you measure it.

   Barry

>   Thanks.                            Jim
>  
> 
> >Note that at iteration 171 the preconditioned residual is 9.348832909193e-13 
> >< 1.e-12 * 1.145582415879e+00  very good convergence.
> >
> >You seem to have set an unreasonably tight convergence criteria. In double 
> >precision you can never realistically expect to use a rtol smaller than 
> >e-12. In fact normally it is not reasonable to use more than like 1.e-8. 
> >Those extra digits don't mean anything.
> >
> >   Barry
> >
> >
> 
> 

Reply via email to