| /petsc/src/snes/impls/ngmres/ |
| H A D | anderson.c | 27 Vec XA, FA, XM, FM; in SNESSolve_Anderson() local
53 FM = snes->work[4]; in SNESSolve_Anderson()
111 PetscCall(SNESGetNPCFunction(snes, FM, &fMnorm)); in SNESSolve_Anderson()
114 PetscCall(VecCopy(F, FM)); in SNESSolve_Anderson()
115 PetscCall(VecWAXPY(XM, -ngmres->andersonBeta, FM, X)); in SNESSolve_Anderson()
119 PetscCall(SNESNGMRESFormCombinedSolution_Private(snes, ivec, l, XM, FM, fMnorm, X, XA, FA)); in SNESSolve_Anderson()
122 …PetscCall(SNESNGMRESNorms_Private(snes, l, X, F, XM, FM, XA, FA, D, &dnorm, &dminnorm, NULL, NULL,… in SNESSolve_Anderson()
128 …PetscCall(SNESNGMRESNorms_Private(snes, l, X, F, XM, FM, XA, FA, D, NULL, NULL, NULL, NULL, NULL, … in SNESSolve_Anderson()
134 …PetscCall(SNESNGMRESNorms_Private(snes, l, X, F, XM, FM, XA, FA, D, NULL, NULL, NULL, NULL, NULL, … in SNESSolve_Anderson()
146 PetscCall(SNESNGMRESUpdateSubspace_Private(snes, ivec, l, FM, fnorm, XM)); in SNESSolve_Anderson()
|
| H A D | ngmresfunc.c | 38 …CombinedSolution_Private(SNES snes, PetscInt ivec, PetscInt l, Vec XM, Vec FM, PetscReal fMnorm, V… in SNESNGMRESFormCombinedSolution_Private() argument
56 PetscCall(VecMDotBegin(FM, l, Fdot, xi)); in SNESNGMRESFormCombinedSolution_Private()
58 PetscCall(VecMDotEnd(FM, l, Fdot, xi)); in SNESNGMRESFormCombinedSolution_Private()
109 PetscCall(VecAXPBY(FA, 1.0 - alph_total, 0.0, FM)); in SNESNGMRESFormCombinedSolution_Private()
115 PetscErrorCode SNESNGMRESNorms_Private(SNES snes, PetscInt l, Vec X, Vec F, Vec XM, Vec FM, Vec XA,… in SNESNGMRESNorms_Private() argument
124 if (fMnorm) PetscCall(VecNormBegin(FM, NORM_2, fMnorm)); in SNESNGMRESNorms_Private()
146 if (fMnorm) PetscCall(VecNormEnd(FM, NORM_2, fMnorm)); in SNESNGMRESNorms_Private()
163 PetscErrorCode SNESNGMRESSelect_Private(SNES snes, PetscInt k_restart, Vec XM, Vec FM, PetscReal xM… in SNESNGMRESSelect_Private() argument
175 PetscCall(VecCopy(FM, F)); in SNESNGMRESSelect_Private()
178 …PetscCall(VecDotRealPart(FM, Y, &ngmres->descent_ls_test)); /* this is actually < F(XM), XM - XA >… in SNESNGMRESSelect_Private()
[all …]
|
| H A D | snesngmres.c | 134 Vec XA, FA, XM, FM; /* candidate linear combination answers */ in SNESSolve_NGMRES() local
166 FM = snes->work[4]; in SNESSolve_NGMRES()
226 PetscCall(SNESGetNPCFunction(snes, FM, &fMnorm)); in SNESSolve_NGMRES()
230 PetscCall(VecCopy(F, FM)); in SNESSolve_NGMRES()
235 PetscCall(SNESLineSearchApply(snes->linesearch, XM, FM, &fMnorm, Y)); in SNESSolve_NGMRES()
244 PetscCall(SNESNGMRESFormCombinedSolution_Private(snes, ivec, l, XM, FM, fMnorm, X, XA, FA)); in SNESSolve_NGMRES()
248 …PetscCall(SNESNGMRESNorms_Private(snes, l, X, F, XM, FM, XA, FA, D, &dnorm, &dminnorm, &xMnorm, NU… in SNESSolve_NGMRES()
250 …PetscCall(SNESNGMRESNorms_Private(snes, l, X, F, XM, FM, XA, FA, D, NULL, NULL, &xMnorm, NULL, &yM… in SNESSolve_NGMRES()
259 …PetscCall(SNESNGMRESSelect_Private(snes, k_restart, XM, FM, xMnorm, fMnorm, yMnorm, objM, XA, FA, … in SNESSolve_NGMRES()
290 PetscCall(SNESNGMRESUpdateSubspace_Private(snes, 0, 0, FM, fMnorm, XM)); in SNESSolve_NGMRES()
[all …]
|