Searched refs:YNtmp (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/snes/impls/ntrdc/ |
| H A D | ntrdc.c | 298 Vec X, F, Y, G, W, GradF, YNtmp; in SNESSolve_NEWTONTRDC() local
322 YNtmp = snes->work[4]; /* Newton solution */ in SNESSolve_NEWTONTRDC()
327 PetscCall(VecGetBlockSize(YNtmp, &bs)); in SNESSolve_NEWTONTRDC()
375 PetscCall(KSPSolve(snes->ksp, F, YNtmp)); /* Quasi Newton Solution */ in SNESSolve_NEWTONTRDC()
384 PetscCall(VecStrideNormAll(YNtmp, NORM_INFINITY, inorms)); in SNESSolve_NEWTONTRDC()
390 PetscCall(VecStrideScale(YNtmp, j, 1.0 / inorms[j])); in SNESSolve_NEWTONTRDC()
405 PetscCall(VecNorm(YNtmp, NORM_2, &ynnorm)); /* ynnorm <- || Y_newton || */ in SNESSolve_NEWTONTRDC()
411 PetscCall(VecCopy(YNtmp, Y)); in SNESSolve_NEWTONTRDC()
429 PetscCall(VecAXPY(YNtmp, -1.0, YCtmp)); /* YCtmp = A, YNtmp = B */ in SNESSolve_NEWTONTRDC()
430 PetscCall(VecNorm(YNtmp, NORM_2, &c0)); /* this could be improved */ in SNESSolve_NEWTONTRDC()
[all …]
|