<div><font face="&#39;courier new&#39;, monospace">[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</font></div>
<div><br></div><div>The relative tolerance can&#39;t possibly be what you want, sorry I didn&#39;t notice this earlier. </div><div><br></div><div><font face="&#39;courier new&#39;, monospace">Linear solve converged due to CONVERGED_RTOL iterations 0</font></div>
<div><font face="&#39;courier new&#39;, monospace">[0] SNESSolve_LS(): iter=0, linear solve iterations=0</font></div><div><font face="&#39;courier new&#39;, monospace">[0] SNESLSCheckResidual_Private(): ||J^T(F-Ax)||/||F-AX|| 1 near zero implies inconsistent rhs</font></div>
<div><font face="&#39;courier new&#39;, monospace">[0] SNESLineSearchCubic(): Search direction and size is 0</font></div><div><font face="&#39;courier new&#39;, monospace">[0] SNESSolve_LS(): fnorm=2.3699692311947174e-02, gnorm=2.3699692311947174e-02, ynorm=0.0000000000000000e+00, lssucceed=0</font></div>
<div><font face="&#39;courier new&#39;, monospace">[0] SNESLSCheckLocalMin_Private(): || J^T F|| 1 near zero implies found a local minimum</font></div><div><font face="&#39;courier new&#39;, monospace">SNES Object:</font></div>
<div><font face="&#39;courier new&#39;, monospace">  type: ls</font></div><div><font face="&#39;courier new&#39;, monospace">    line search variant: SNESLineSearchCubic</font></div><div><font face="&#39;courier new&#39;, monospace">    alpha=0.0001, maxstep=1e+08, minlambda=1e-12</font></div>
<div><font face="&#39;courier new&#39;, monospace">  maximum iterations=50, maximum function evaluations=10000</font></div><div><font face="&#39;courier new&#39;, monospace">  tolerances: relative=1e-08, absolute=1e-50, solution=1e-08</font></div>
<div><font face="&#39;courier new&#39;, monospace">  total number of linear solver iterations=0</font></div><div><font face="&#39;courier new&#39;, monospace">  total number of function evaluations=1</font></div><div><font face="&#39;courier new&#39;, monospace">  KSP Object:</font></div>
<div><font face="&#39;courier new&#39;, monospace">    type: gmres</font></div><div><font face="&#39;courier new&#39;, monospace">      GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement</font></div>
<div><font face="&#39;courier new&#39;, monospace">      GMRES: happy breakdown tolerance 1e-30</font></div><div><font face="&#39;courier new&#39;, monospace">    maximum iterations=20, initial guess is zero</font></div><div>
<font face="&#39;courier new&#39;, monospace">    tolerances:  relative=800695, absolute=5.33098e-35, divergence=5.33098e-35</font></div><div><br></div><div><font face="arial, helvetica, sans-serif">How are you setting the linear solve tolerance? Perhaps you are passing arguments of the incorrect type to KSPSetTolerances()? (Fortran doesn&#39;t warn about this, sadly.)</font></div>