Searched refs:gTBg (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/snes/impls/tr/ |
| H A D | tr.c | 492 PetscReal auk, tauk, gfnorm, gfnorm_k, ycnorm, gTBg, objmin = 0.0, beta_k = 1.0; in SNESSolve_NEWTONTR() local
631 if (has_objective) PetscCall(VecDotRealPart(GradF, W, &gTBg)); in SNESSolve_NEWTONTR()
632 else PetscCall(VecDotRealPart(W, W, &gTBg)); /* B = J^t * J */ in SNESSolve_NEWTONTR()
635 if (gTBg < 0.0) tauk = 1.0; in SNESSolve_NEWTONTR()
636 else tauk = PetscMin(gfnorm * gfnorm * gfnorm_k / (delta * gTBg), 1); in SNESSolve_NEWTONTR()
668 …ted on boundary. delta: %g, ycnorm: %g, gTBg: %g\n", (double)delta, (double)ycnorm, (double)gTBg)); in SNESSolve_NEWTONTR()
685 if (fallback == SNES_TR_FALLBACK_CAUCHY || gTBg <= 0.0) { in SNESSolve_NEWTONTR()
687 …, ynorm: %g, ycnorm: %g, gTBg: %g\n", (double)delta, (double)ynorm, (double)ycnorm, (double)gTBg)); in SNESSolve_NEWTONTR()
689 auk = gfnorm * gfnorm / gTBg; in SNESSolve_NEWTONTR()
726 …CP" : "DL", (double)tneg, (double)tpos, (double)tau, (double)ynorm, (double)ycnorm, (double)gTBg)); in SNESSolve_NEWTONTR()
|
| /petsc/src/snes/impls/ntrdc/ |
| H A D | ntrdc.c | 305 …PetscReal auk, gfnorm, ycnorm, c0, c1, c2, tau, tau_pos, tau_neg, gTBg; /*… in SNESSolve_NEWTONTRDC() local
414 PetscCall(VecDotRealPart(W, W, &gTBg)); /* completes GradF^T J^T J GradF */ in SNESSolve_NEWTONTRDC()
416 if (gTBg <= 0.0) { in SNESSolve_NEWTONTRDC()
419 auk = PetscSqr(gfnorm) / gTBg; in SNESSolve_NEWTONTRDC()
|