46 static PetscErrorCode SNESNEWTONLSCheckLocalMin_Private(SNES snes, Mat A, Vec F, PetscReal fnorm, P…  in SNESNEWTONLSCheckLocalMin_Private()  argument
63 …cInfo(snes, "|| J^T F|| %14.12e near zero implies found a local minimum\n", (double)(a1 / fnorm)));  in SNESNEWTONLSCheckLocalMin_Private()
64       if (a1 / fnorm < 1.e-4) *ismin = PETSC_TRUE;  in SNESNEWTONLSCheckLocalMin_Private()
76       a1 = PetscAbsScalar(result) / (fnorm * wnorm);  in SNESNEWTONLSCheckLocalMin_Private()
129   PetscReal            fnorm, xnorm, ynorm;  in SNESSolve_NEWTONLS()  local
165     PetscCall(VecNormBegin(F, NORM_2, &fnorm));  in SNESSolve_NEWTONLS()
166     PetscCall(VecNormEnd(F, NORM_2, &fnorm));  in SNESSolve_NEWTONLS()
172   PetscCall(VecNorm(F, NORM_2, &fnorm)); /* fnorm <- ||F||  */  in SNESSolve_NEWTONLS()
173   SNESCheckFunctionDomainError(snes, fnorm);  in SNESSolve_NEWTONLS()
175   snes->norm = fnorm;  in SNESSolve_NEWTONLS()
177   PetscCall(SNESLogConvergenceHistory(snes, fnorm, 0));  in SNESSolve_NEWTONLS()
180   PetscCall(SNESConverged(snes, 0, 0.0, 0.0, fnorm));  in SNESSolve_NEWTONLS()
181   PetscCall(SNESMonitor(snes, 0, fnorm));  in SNESSolve_NEWTONLS()
191 …PetscCall(VecNorm(snes->vec_func, NORM_2, &fnorm)); /* no-op unless update() function changed f() …  in SNESSolve_NEWTONLS()
205         PetscCall(SNESGetNPCFunction(snes, F, &fnorm));  in SNESSolve_NEWTONLS()
228     gnorm = fnorm;  in SNESSolve_NEWTONLS()
234     PetscCall(SNESLineSearchApply(linesearch, X, F, &fnorm, Y));  in SNESSolve_NEWTONLS()
236     PetscCall(SNESLineSearchGetNorms(linesearch, &xnorm, &fnorm, &ynorm));  in SNESSolve_NEWTONLS()
237 …16e, gnorm=%18.16e, ynorm=%18.16e, lsreason=%d\n", (double)gnorm, (double)fnorm, (double)ynorm, (i…  in SNESSolve_NEWTONLS()
238     SNESCheckFunctionDomainError(snes, fnorm);  in SNESSolve_NEWTONLS()
260         PetscCall(SNESNEWTONLSCheckLocalMin_Private(snes, snes->jacobian, F, fnorm, &ismin));  in SNESSolve_NEWTONLS()
268     snes->norm  = fnorm;  in SNESSolve_NEWTONLS()
274     PetscCall(SNESConverged(snes, snes->iter, xnorm, ynorm, fnorm));  in SNESSolve_NEWTONLS()