Searched refs:fkp1 (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/snes/impls/tr/ |
| H A D | tr.c | 415 …(SNES snes, PetscBool has_objective, Vec X, Vec Y, Vec W, Vec G, PetscReal *gnorm, PetscReal *fkp1) in SNESNewtonTRObjective() argument
430 PetscCall(SNESComputeObjective(snes, W, fkp1)); in SNESNewtonTRObjective()
431 SNESCheckObjectiveDomainError(snes, *fkp1); in SNESNewtonTRObjective()
432 } else *fkp1 = 0.5 * PetscSqr(*gnorm); in SNESNewtonTRObjective()
491 PetscReal fk, fkp1, deltaqm = 0.0, gTy = 0.0, yTHy = 0.0; in SNESSolve_NEWTONTR() local
740 PetscCall(SNESNewtonTRObjective(snes, has_objective, X, Y, W, G, &gnorm, &fkp1)); in SNESSolve_NEWTONTR()
741 if (PetscIsInfOrNanReal(fkp1)) rho = neP->eta1; in SNESSolve_NEWTONTR()
743 …if (deltaqm > 0.0) rho = (fk - fkp1) / deltaqm; /* actual improvement over predicted improvement */ in SNESSolve_NEWTONTR()
748 …g, yTHy=%g, ynormk=%g\n", (double)rho, (double)delta, (double)fk, (double)fkp1, (double)deltaqm, (… in SNESSolve_NEWTONTR()
777 fk = fkp1; in SNESSolve_NEWTONTR()
|