| /petsc/src/snes/impls/ngmres/ |
| H A D | anderson.c | 27 Vec XA, FA, XM, FM; in SNESSolve_Anderson() local
52 XM = snes->work[3]; in SNESSolve_Anderson()
99 PetscCall(VecCopy(X, XM)); in SNESSolve_Anderson()
102 PetscCall(PetscLogEventBegin(SNES_NPCSolve, snes->npc, XM, B, 0)); in SNESSolve_Anderson()
103 PetscCall(SNESSolve(snes->npc, B, XM)); in SNESSolve_Anderson()
104 PetscCall(PetscLogEventEnd(SNES_NPCSolve, snes->npc, XM, B, 0)); in SNESSolve_Anderson()
112 PetscCall(VecAXPBY(XM, 1.0 - ngmres->andersonBeta, ngmres->andersonBeta, X)); 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()
[all …]
|
| H A D | ngmresfunc.c | 38 PetscErrorCode SNESNGMRESFormCombinedSolution_Private(SNES snes, PetscInt ivec, PetscInt l, Vec XM,… in SNESNGMRESFormCombinedSolution_Private() argument
97 PetscCall(VecAXPBY(XA, 1.0 - alph_total, 0.0, XM)); 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
123 if (xMnorm) PetscCall(VecNormBegin(XM, NORM_2, xMnorm)); in SNESNGMRESNorms_Private()
126 PetscCall(VecWAXPY(D, -1.0, XM, X)); in SNESNGMRESNorms_Private()
136 PetscCall(VecWAXPY(D, -1.0, XM, XA)); in SNESNGMRESNorms_Private()
145 if (xMnorm) PetscCall(VecNormEnd(XM, NORM_2, xMnorm)); in SNESNGMRESNorms_Private()
163 PetscErrorCode SNESNGMRESSelect_Private(SNES snes, PetscInt k_restart, Vec XM, Vec FM, PetscReal xM… in SNESNGMRESSelect_Private() argument
176 PetscCall(VecCopy(XM, X)); in SNESNGMRESSelect_Private()
184 PetscCall(VecWAXPY(Y, -1.0, X, XM)); in SNESNGMRESSelect_Private()
[all …]
|
| H A D | snesngmres.c | 134 Vec XA, FA, XM, FM; /* candidate linear combination answers */ in SNESSolve_NGMRES() local
165 XM = snes->work[3]; in SNESSolve_NGMRES()
214 PetscCall(VecCopy(X, XM)); in SNESSolve_NGMRES()
217 PetscCall(PetscLogEventBegin(SNES_NPCSolve, snes->npc, XM, B, 0)); in SNESSolve_NGMRES()
218 PetscCall(SNESSolve(snes->npc, B, XM)); in SNESSolve_NGMRES()
219 PetscCall(PetscLogEventEnd(SNES_NPCSolve, snes->npc, XM, B, 0)); in SNESSolve_NGMRES()
231 PetscCall(VecCopy(X, XM)); in SNESSolve_NGMRES()
235 PetscCall(SNESLineSearchApply(snes->linesearch, XM, FM, &fMnorm, Y)); in SNESSolve_NGMRES()
239 PetscCall(SNESComputeObjective(snes, XM, &objM)); in SNESSolve_NGMRES()
244 PetscCall(SNESNGMRESFormCombinedSolution_Private(snes, ivec, l, XM, FM, fMnorm, X, XA, FA)); in SNESSolve_NGMRES()
[all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex40f90.F90 | 20 PetscScalar f(in%DOF, in%XS + 1:in%XS + in%XM, in%YS + 1:in%YS + in%YM)
23 do i = in%XS + 1, in%XS + in%XM
|
| H A D | ex5f.F90 | 175 xe = xs + info%XM - 1
|