[petsc-users] Problem with SNES

Jed Brown jedbrown at mcs.anl.gov
Mon Aug 29 11:50:59 CDT 2011


[0] KSPDefaultConverged(): Linear solver has converged. Residual norm
0.0236997 is less than relative tolerance 800695 times initial right hand
side norm 0.0236997 at iteration 0

The relative tolerance can't possibly be what you want, sorry I didn't
notice this earlier.

Linear solve converged due to CONVERGED_RTOL iterations 0
[0] SNESSolve_LS(): iter=0, linear solve iterations=0
[0] SNESLSCheckResidual_Private(): ||J^T(F-Ax)||/||F-AX|| 1 near zero
implies inconsistent rhs
[0] SNESLineSearchCubic(): Search direction and size is 0
[0] SNESSolve_LS(): fnorm=2.3699692311947174e-02,
gnorm=2.3699692311947174e-02, ynorm=0.0000000000000000e+00, lssucceed=0
[0] SNESLSCheckLocalMin_Private(): || J^T F|| 1 near zero implies found a
local minimum
SNES Object:
  type: ls
    line search variant: SNESLineSearchCubic
    alpha=0.0001, maxstep=1e+08, minlambda=1e-12
  maximum iterations=50, maximum function evaluations=10000
  tolerances: relative=1e-08, absolute=1e-50, solution=1e-08
  total number of linear solver iterations=0
  total number of function evaluations=1
  KSP Object:
    type: gmres
      GMRES: restart=30, using Classical (unmodified) Gram-Schmidt
Orthogonalization with no iterative refinement
      GMRES: happy breakdown tolerance 1e-30
    maximum iterations=20, initial guess is zero
    tolerances:  relative=800695, absolute=5.33098e-35,
divergence=5.33098e-35

How are you setting the linear solve tolerance? Perhaps you are passing
arguments of the incorrect type to KSPSetTolerances()? (Fortran doesn't warn
about this, sadly.)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110829/b57166d6/attachment.htm>


More information about the petsc-users mailing list