| /petsc/src/mat/tests/output/ |
| H A D | ex202_1.out | 76 <Ax, y> = 6.96084
78 <Ax, y> = 7.92757
80 <Ax, y> = 6.76323
82 <Ax, y> = 5.58672
84 <Ax, y> = 3.77397
86 <Ax, y> = 5.89481
88 <Ax, y> = 1.9272
90 <Ax, y> = 6.9047
92 <Ax, y> = 5.0697
94 <Ax, y> = 3.59347
|
| /petsc/src/mat/utils/ |
| H A D | multequal.c | 11 Vec Ax = NULL, Bx = NULL, s1 = NULL, s2 = NULL, Ay = NULL, By = NULL; in MatMultEqual_Private() local
38 PetscCall(MatCreateVecs(A, &s1, &Ax)); in MatMultEqual_Private()
41 PetscCall(MatCreateVecs(A, &Ax, &s1)); in MatMultEqual_Private()
53 PetscCall(VecSetRandom(Ax, rctx)); in MatMultEqual_Private()
54 PetscCall(VecCopy(Ax, Bx)); in MatMultEqual_Private()
69 PetscCall(MatMultTransposeAdd(A, Ax, Aadd, s1)); in MatMultEqual_Private()
72 PetscCall(MatMultTranspose(A, Ax, s1)); in MatMultEqual_Private()
77 PetscCall(MatMultHermitianTransposeAdd(A, Ax, Aadd, s1)); in MatMultEqual_Private()
80 PetscCall(MatMultHermitianTranspose(A, Ax, s1)); in MatMultEqual_Private()
85 PetscCall(MatMultAdd(A, Ax, Aadd, s1)); in MatMultEqual_Private()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex202.c | 13 Vec x, y, Ax, ATy; in TestInitialMatrix() local
30 PetscCall(VecDuplicate(y, &Ax)); in TestInitialMatrix()
51 PetscCall(MatMult(A, x, Ax)); in TestInitialMatrix()
52 PetscCall(VecDot(Ax, y, &dot1)); in TestInitialMatrix()
61 PetscCall(VecDestroy(&Ax)); in TestInitialMatrix()
75 PetscCall(MatCreateVecs(A, NULL, &Ax)); in TestInitialMatrix()
76 PetscCall(MatMult(A, x, Ax)); in TestInitialMatrix()
80 PetscCall(VecAXPY(y, -1.0, Ax)); in TestInitialMatrix()
81 PetscCall(VecDestroy(&Ax)); in TestInitialMatrix()
|
| H A D | ex66.c | 33 Vec v, x, y, Ax, Ay, Bx, By; in main() local
159 PetscCall(MatCreateVecs(A, &Ax, &Ay)); in main()
161 PetscCall(VecSetRandom(Ax, NULL)); in main()
162 PetscCall(VecCopy(Ax, Bx)); in main()
163 PetscCall(MatMult(A, Ax, Ay)); in main()
216 PetscCall(MatMultTranspose(A, Ay, Ax)); in main()
218 PetscCall(VecViewFromOptions(Ax, NULL, "-multtrans_vec_view")); in main()
220 PetscCall(VecNorm(Ax, NORM_INFINITY, &nX)); in main()
221 PetscCall(VecAXPY(Ax, -1.0, Bx)); in main()
222 PetscCall(VecViewFromOptions(Ax, NULL, "-multtrans_vec_view")); in main()
[all …]
|
| H A D | ex70.c | 67 Vec Ax, Ay; in proj_mult() local
80 PetscCall(MatCreateVecs(A, &Ax, &Ay)); in proj_mult()
92 PetscCall(VecCopy(Py, Ax)); in proj_mult()
93 PetscCall(MatMult(A, Ax, Ay)); in proj_mult()
103 PetscCall(VecDestroy(&Ax)); in proj_mult()
|
| /petsc/src/tao/constrained/impls/admm/ |
| H A D | admm.c | 49 PetscCall(VecNorm(am->Ax, NORM_2, &Axnorm)); in TaoADMMToleranceUpdate()
77 PetscCall(VecWAXPY(tempJR, -1., am->Axold, am->Ax)); in AdaptiveADMMPenaltyUpdate()
152 static PetscErrorCode ADMMUpdateConstraintResidualVector(Tao tao, Vec x, Vec z, Vec Ax, Vec Bz, Vec… in ADMMUpdateConstraintResidualVector() argument
161 PetscCall(MatMult(mis->jacobian_equality, x, Ax)); in ADMMUpdateConstraintResidualVector()
165 PetscCall(VecWAXPY(residual, 1., Bz, Ax)); in ADMMUpdateConstraintResidualVector()
182 …PetscCall(ADMMUpdateConstraintResidualVector(parent, x, am->subsolverZ->solution, am->Ax, am->Bz, … in SubObjGradUpdate()
211 …PetscCall(ADMMUpdateConstraintResidualVector(parent, am->subsolverX->solution, z, am->Ax, am->Bz, … in RegObjGradUpdate()
414 PetscCall(MatMult(am->subsolverX->jacobian_equality, am->subsolverX->solution, am->Ax)); in TaoSolve_ADMM()
425 PetscCall(VecWAXPY(am->workJacobianRight, 1 / am->mu, am->y, am->Ax)); in TaoSolve_ADMM()
431 …esidualVector(tao, am->subsolverX->solution, am->subsolverZ->solution, am->Ax, am->Bz, am->residua… in TaoSolve_ADMM()
[all …]
|
| H A D | admm.h | 19 Vec z, zold, Ax, Bz, Axold, Bzold, Bz0; member
|
| /petsc/src/ksp/ksp/guess/impls/fischer/ |
| H A D | fischer.c | 12 Vec Ax; /* method 2 */ member
38 PetscCall(VecDestroy(&itg->Ax)); in KSPGuessReset_Fischer()
55 if (!itg->Ax && itg->method == 2) PetscCall(VecDuplicate(itg->xtilde[0], &itg->Ax)); in KSPGuessSetUp_Fischer()
71 PetscCall(VecDestroy(&itg->Ax)); in KSPGuessDestroy_Fischer()
165 PetscCall(KSP_MatMult(guess->ksp, guess->A, x, itg->Ax)); /* norm = sqrt(x'Ax) */ in KSPGuessUpdate_Fischer_2()
166 PetscCall(VecDot(x, itg->Ax, &norm)); in KSPGuessUpdate_Fischer_2()
176 PetscCall(KSP_MatMult(guess->ksp, guess->A, itg->xtilde[curl], itg->Ax)); in KSPGuessUpdate_Fischer_2()
177 PetscCall(VecMDot(itg->Ax, curl, itg->xtilde, itg->alpha)); in KSPGuessUpdate_Fischer_2()
181 …PetscCall(KSP_MatMult(guess->ksp, guess->A, itg->xtilde[curl], itg->Ax)); /* norm = sqrt(xtilde[cu… in KSPGuessUpdate_Fischer_2()
182 PetscCall(VecDot(itg->xtilde[curl], itg->Ax, &norm)); in KSPGuessUpdate_Fischer_2()
[all …]
|
| /petsc/src/ksp/ksp/tutorials/output/ |
| H A D | ex56_seqaijmkl.out | 4 [0]main |b-Ax|/|b|=9.264005e-05, |b|=5.391826e+00, emax=9.987464e-01
|
| H A D | ex56_nns.out | 4 [0]main |b-Ax|/|b|=4.543854e-05, |b|=5.391826e+00, emax=9.949156e-01
|
| H A D | ex56_2.out | 4 [0]main |b-Ax|/|b|=9.023568e-05, |b|=4.969822e+00, emax=1.016793e+00
|
| H A D | ex56_hypre.out | 4 [0]main |b-Ax|/|b|=9.722297e-05, |b|=4.351790e+00, emax=9.991454e-01
|
| H A D | ex56_gamg.out | 4 [0]main |b-Ax|/|b|=4.993291e-05, |b|=4.351790e+00, emax=9.986640e-01
|
| H A D | ex56_1.out | 4 [0]main |b-Ax|/|b|=9.081303e-05, |b|=4.969822e+00, emax=9.982432e-01
|
| H A D | ex56_latebs-2.out | 285 [0]main |b-Ax|/|b|=2.425235e-04, |b|=5.391826e+00, emax=9.946388e-01
|
| H A D | ex56_latebs.out | 496 [0]main |b-Ax|/|b|=1.165490e-04, |b|=5.391826e+00, emax=9.990809e-01
|
| /petsc/src/ksp/ksp/tests/output/ |
| H A D | ex63_1.out | 4 Norm2 of Ax - b = {3.88578058618804789e-16}
|
| H A D | ex63_2_alt.out | 16 Norm2 of Ax - b = {1.11022302462515654e-16}
|
| H A D | ex63_2.out | 16 Norm2 of Ax - b = {0.00000000000000000e+00}
|
| /petsc/src/ksp/ksp/utils/lmvm/ |
| H A D | lmbasis.c | 489 PETSC_INTERN PetscErrorCode LMBasisSetCachedProduct(LMBasis A, Vec x, Vec Ax) in LMBasisSetCachedProduct() argument
499 PetscCall(PetscObjectReference((PetscObject)Ax)); in LMBasisSetCachedProduct()
501 A->cached_product = Ax; in LMBasisSetCachedProduct()
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | KSP.pyx | 239 the preconditioned residual P⁻¹(b - Ax).
241 Uses the l₂ norm of the true b - Ax residual.
273 ∥r∥ <= rtolnorm(b) or rtolnorm(b - Ax₀)
1755 Ax = b where b = bₚ + bₜ where bₜ is not in the range of A
1759 KSP first removes bₜ producing the linear system Ax = bₚ (which
|
| H A D | Mat.pyx | 5485 """Solve Ax=b, given a factored matrix.
|
| /petsc/doc/manual/ |
| H A D | ksp.md | 276 r \equiv b - Ax = b - A M_R^{-1} \, M_R x,
454 …rm of the preconditioned \$ B(b - A x) \$ or unconditioned residual \$ b - Ax\$, depending on the …
566 formula $r = b - Ax$, the routine is slow and should be used only
1908 the usual formula $b - Ax$, the user may wish to provide an
|
| H A D | tao.md | 1966 \text{subject to} & Ax + Bz = c
1978 L_{\rho}(x,z,y) = f(x) + g(z) + y^T(Ax+Bz-c) + (\rho/2)||Ax+Bz-c||_2^2
1992 y^{k+1} := y^k + \rho(Ax^{k+1}+Bz^{k+1}-c)
|
| H A D | getting_started.md | 402 `Ax` $=$ `b`, the user can then use `KSP` to solve the
|