Searched refs:yTHy (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
403 if (has_objective) PetscCall(VecDotRealPart(Y, W, &yTHy)); in SNESNewtonTRQuadraticDelta()
404 else PetscCall(VecDotRealPart(W, W, &yTHy)); /* Gauss-Newton approximation J^t * J */ in SNESNewtonTRQuadraticDelta()
406 …*deltaqm = -(-(gTy) + 0.5 * (yTHy)); /* difference in quadratic model, -gTy because SNES solves it… in SNESNewtonTRQuadraticDelta()
407 if (yTHy_) *yTHy_ = yTHy; 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 …ouble)delta, (double)fk, (double)fkp1, (double)deltaqm, (double)gTy, (double)yTHy, (double)ynorm)); 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
454 PetscCall(VecDotRealPart(W, W, &yTHy)); /* completes GradY^T J^T J GradY */ in SNESSolve_NEWTONTRDC()
456 mp = f0 - gTy + 0.5 * yTHy; /* quadratic model to satisfy, -gTy because our update is X-Y*/ in SNESSolve_NEWTONTRDC()
|