Searched refs:deltaqm (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/snes/impls/tr/ |
| H A D | tr.c | 397 …Bool has_objective, Vec Y, Vec GradF, Vec W, PetscReal *yTHy_, PetscReal *gTy_, PetscReal *deltaqm) in SNESNewtonTRQuadraticDelta() argument
406 …*deltaqm = -(-(gTy) + 0.5 * (yTHy)); /* difference in quadratic model, -gTy because SNES solves it… in SNESNewtonTRQuadraticDelta()
491 PetscReal fk, fkp1, deltaqm = 0.0, gTy = 0.0, yTHy = 0.0; in SNESSolve_NEWTONTR() local
737 … PetscCall(SNESNewtonTRQuadraticDelta(snes, J, has_objective, Y, GradF, W, &yTHy, &gTy, &deltaqm)); in SNESSolve_NEWTONTR()
743 …if (deltaqm > 0.0) rho = (fk - fkp1) / deltaqm; /* actual improvement over predicted improvement */ in SNESSolve_NEWTONTR()
748 …rmk=%g\n", (double)rho, (double)delta, (double)fk, (double)fkp1, (double)deltaqm, (double)gTy, (do… in SNESSolve_NEWTONTR()
|