Searched refs:GradF (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/snes/impls/tr/ |
| H A D | tr.c | 397 …tonTRQuadraticDelta(SNES snes, Mat J, PetscBool has_objective, Vec Y, Vec GradF, Vec W, PetscReal … in SNESNewtonTRQuadraticDelta() argument
405 PetscCall(VecDotRealPart(GradF, Y, &gTy)); in SNESNewtonTRQuadraticDelta()
488 Vec X, F, Y, G, W, GradF, YU, Yc; in SNESSolve_NEWTONTR() local
513 GradF = !has_objective ? snes->work[2] : snes->vec_func; /* grad f = J^T F */ in SNESSolve_NEWTONTR()
620 PetscCall(MatMultTranspose(J, F, GradF)); /* grad f = J^T F */ in SNESSolve_NEWTONTR()
621 PetscCall(VecNorm(GradF, NORM_2, &gfnorm)); in SNESSolve_NEWTONTR()
623 PetscCall(VecNorm(GradF, neP->norm, &gfnorm_k)); in SNESSolve_NEWTONTR()
630 PetscCall(MatMult(J, GradF, W)); in SNESSolve_NEWTONTR()
631 if (has_objective) PetscCall(VecDotRealPart(GradF, W, &gTBg)); in SNESSolve_NEWTONTR()
639 PetscCall(VecAXPBY(Yc, auk, 0.0, GradF)); in SNESSolve_NEWTONTR()
[all …]
|
| /petsc/src/snes/impls/ntrdc/ |
| H A D | ntrdc.c | 298 Vec X, F, Y, G, W, GradF, YNtmp; in SNESSolve_NEWTONTRDC() local
321 GradF = snes->work[3]; /* grad f = J^T F */ in SNESSolve_NEWTONTRDC()
404 PetscCall(MatMultTranspose(jac, F, GradF)); /* grad f = J^T F */ in SNESSolve_NEWTONTRDC()
413 PetscCall(MatMult(jac, GradF, W)); in SNESSolve_NEWTONTRDC()
415 PetscCall(VecNorm(GradF, NORM_2, &gfnorm)); /* grad f norm <- || grad f || */ in SNESSolve_NEWTONTRDC()
422 PetscCall(VecCopy(GradF, YCtmp)); /* this could be improved */ in SNESSolve_NEWTONTRDC()
455 PetscCall(VecDotRealPart(GradF, Y, &gTy)); in SNESSolve_NEWTONTRDC()
|