| /petsc/src/snes/impls/ngmres/ |
| H A D | ngmresfunc.c | 38 …vate(SNES snes, PetscInt ivec, PetscInt l, Vec XM, Vec FM, PetscReal fMnorm, Vec X, Vec XA, Vec FA) in SNESNGMRESFormCombinedSolution_Private() argument
97 PetscCall(VecAXPBY(XA, 1.0 - alph_total, 0.0, XM)); in SNESNGMRESFormCombinedSolution_Private()
98 PetscCall(VecMAXPY(XA, l, beta, Xdot)); in SNESNGMRESFormCombinedSolution_Private()
100 PetscCall(VecWAXPY(Y, -1.0, X, XA)); in SNESNGMRESFormCombinedSolution_Private()
101 PetscCall(SNESLineSearchPostCheck(snes->linesearch, X, Y, XA, &changed_y, &changed_w)); in SNESNGMRESFormCombinedSolution_Private()
104 PetscCall(SNESApplyNPC(snes, XA, NULL, FA)); in SNESNGMRESFormCombinedSolution_Private()
106 PetscCall(SNESComputeFunction(snes, XA, FA)); 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
129 if (xAnorm) PetscCall(VecNormBegin(XA, NORM_2, xAnorm)); in SNESNGMRESNorms_Private()
132 PetscCall(VecWAXPY(D, -1.0, XA, X)); in SNESNGMRESNorms_Private()
[all …]
|
| H A D | anderson.c | 27 Vec XA, FA, XM, FM; in SNESSolve_Anderson() local
47 XA = snes->work[2]; 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()
152 PetscCall(VecCopy(XA, X)); in SNESSolve_Anderson()
|
| H A D | snesngmres.c | 134 Vec XA, FA, XM, FM; /* candidate linear combination answers */ in SNESSolve_NGMRES() local
159 XA = snes->work[2]; 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()
254 PetscCall(SNESComputeObjective(snes, XA, &objA)); in SNESSolve_NGMRES()
259 …PetscCall(SNESNGMRESSelect_Private(snes, k_restart, XM, FM, xMnorm, fMnorm, yMnorm, objM, XA, FA, … in SNESSolve_NGMRES()
|