Searched refs:gTy (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/snes/impls/tr/ |
| H A D | tr.c | 399 PetscReal yTHy, gTy; in SNESNewtonTRQuadraticDelta() local
405 PetscCall(VecDotRealPart(GradF, Y, &gTy)); in SNESNewtonTRQuadraticDelta()
406 …*deltaqm = -(-(gTy) + 0.5 * (yTHy)); /* difference in quadratic model, -gTy because SNES solves it… in SNESNewtonTRQuadraticDelta()
408 if (gTy_) *gTy_ = gTy; 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()
748 …ble)rho, (double)delta, (double)fk, (double)fkp1, (double)deltaqm, (double)gTy, (double)yTHy, (dou… in SNESSolve_NEWTONTR()
|
| /petsc/src/snes/impls/ntrdc/ |
| H A D | ntrdc.c | 304 …PetscReal deltaM, ynnorm, f0, mp, gTy, g, yTHy; /*… in SNESSolve_NEWTONTRDC() local
455 PetscCall(VecDotRealPart(GradF, Y, &gTy)); in SNESSolve_NEWTONTRDC()
456 mp = f0 - gTy + 0.5 * yTHy; /* quadratic model to satisfy, -gTy because our update is X-Y*/ in SNESSolve_NEWTONTRDC()
|