| /petsc/src/mat/tests/ |
| H A D | ex88.c | 52 Vec W1, W2, W3, diff; in TestMatrix() local
71 PetscCall(VecDuplicate(X, &W2)); in TestMatrix()
83 PetscCall(MatMultTranspose(A, W1, W2)); in TestMatrix()
84 PetscCall(VecView(W2, viewer)); in TestMatrix()
87 PetscCall(MatMultHermitianTranspose(A, W1, W2)); in TestMatrix()
88 PetscCall(VecConjugate(W2)); in TestMatrix()
89 PetscCall(VecView(W2, viewer)); in TestMatrix()
94 PetscCall(MatMultAdd(A, W1, W1, W2)); in TestMatrix()
95 PetscCall(VecView(W2, viewer)); in TestMatrix()
96 PetscCall(VecAXPY(W2, -1.0, diff)); in TestMatrix()
[all …]
|
| /petsc/src/tao/linesearch/impls/gpcglinesearch/ |
| H A D | gpcglinesearch.c | 10 PetscCall(VecDestroy(&ctx->W2)); in TaoLineSearchDestroy_GPCG()
48 if (!neP->W2) { in TaoLineSearchApply_GPCG()
49 PetscCall(VecDuplicate(x, &neP->W2)); in TaoLineSearchApply_GPCG()
57 PetscCall(VecDestroy(&neP->W2)); in TaoLineSearchApply_GPCG()
60 PetscCall(VecDuplicate(x, &neP->W2)); in TaoLineSearchApply_GPCG()
73 PetscCall(VecCopy(x, neP->W2)); in TaoLineSearchApply_GPCG()
96 PetscCall(VecWAXPY(neP->W2, ls->step, s, x)); in TaoLineSearchApply_GPCG()
99 PetscCall(VecMedian(neP->W2, ls->lower, ls->upper, neP->W2)); in TaoLineSearchApply_GPCG()
107 PetscCall(TaoLineSearchComputeObjective(ls, neP->W2, f)); in TaoLineSearchApply_GPCG()
110 PetscCall(TaoLineSearchComputeObjectiveAndGTS(ls, neP->W2, f, &gdx)); in TaoLineSearchApply_GPCG()
[all …]
|
| H A D | gpcglinesearch.h | 13 Vec W2; member
|
| /petsc/src/ksp/ksp/impls/bcgs/pipebcgs/ |
| H A D | pipebcgs.c | 19 Vec X, B, S, R, RP, Y, Q, P2, Q2, R2, S2, W, Z, W2, Z2, T, V; in KSPSolve_PIPEBCGS() local
38 W2 = ksp->work[11]; in KSPSolve_PIPEBCGS()
83 PetscCall(KSP_PCApply(ksp, W, W2)); /* w2 <- K w */ in KSPSolve_PIPEBCGS()
84 PetscCall(KSP_MatMult(ksp, pc->mat, W2, T)); /* t <- A w2 */ in KSPSolve_PIPEBCGS()
96 PetscCall(VecCopy(W2, S2)); /* s2 <- w2 */ in KSPSolve_PIPEBCGS()
101 …PetscCall(VecAXPBYPCZ(S2, 1.0, -beta * omega, beta, W2, Z2)); /* s2 <- beta * s2 + w2 - beta * ome… in KSPSolve_PIPEBCGS()
138 PetscCall(VecWAXPY(R2, -alpha, Z2, W2)); /* r2 <- w2 - alpha z2 */ in KSPSolve_PIPEBCGS()
151 PetscCall(KSP_PCApply(ksp, W, W2)); /* w2 <- K w */ in KSPSolve_PIPEBCGS()
152 PetscCall(KSP_MatMult(ksp, pc->mat, W2, T)); /* t <- A w2 */ in KSPSolve_PIPEBCGS()
174 PetscCall(KSP_PCApply(ksp, W, W2)); /* w2 <- K w */ in KSPSolve_PIPEBCGS()
[all …]
|
| /petsc/src/ts/tests/ |
| H A D | ex13.c | 21 Vec W, W2, Wdot; in main() local
35 PetscCall(VecDuplicate(W, &W2)); in main()
37 PetscCall(TSSetSolution(ts, W2)); in main()
123 PetscCall(TSSetSolution(ts, W2)); in main()
167 PetscCall(VecDestroy(&W2)); in main()
|
| /petsc/src/snes/impls/ls/ |
| H A D | ls.c | 98 Vec W1, W2; in SNESNEWTONLSCheckResidual_Private() local
101 PetscCall(VecDuplicate(F, &W2)); in SNESNEWTONLSCheckResidual_Private()
106 PetscCall(MatMultTranspose(A, W1, W2)); in SNESNEWTONLSCheckResidual_Private()
108 PetscCall(VecNorm(W2, NORM_2, &a2)); in SNESNEWTONLSCheckResidual_Private()
111 PetscCall(VecDestroy(&W2)); in SNESNEWTONLSCheckResidual_Private()
|
| /petsc/src/tao/matrix/ |
| H A D | adamat.c | 10 Vec W2; member
27 PetscCall(VecPointwiseMult(ctx->W2, ctx->D2, a)); in MatMult_ADA()
28 PetscCall(VecAXPY(y, one, ctx->W2)); in MatMult_ADA()
63 PetscCall(VecDestroy(&ctx->W2)); in MatDestroy_ADA()
370 PetscCall(VecDuplicate(d2, &ctx->W2)); in MatCreateADA()
374 ctx->W2 = NULL; in MatCreateADA()
|
| /petsc/src/ksp/ksp/impls/lsqr/ |
| H A D | lsqr.c | 73 Vec X, B, V, V1, U, U1, TMP, W, W2, Z = NULL; in KSPSolve_LSQR() local
95 W2 = lsqr->vwork_n[3]; in KSPSolve_LSQR()
203 PetscCall(VecCopy(W, W2)); in KSPSolve_LSQR()
204 PetscCall(VecSquare(W2)); in KSPSolve_LSQR()
205 PetscCall(VecScale(W2, 1.0 / (rho * rho))); in KSPSolve_LSQR()
206 PetscCall(VecAXPY(lsqr->se, 1.0, W2)); /* lsqr->se <- lsqr->se + (w^2/rho^2) */ in KSPSolve_LSQR()
|