Home
last modified time | relevance | path

Searched refs:Amat (Results 1 – 25 of 76) sorted by relevance

1234

/petsc/src/ksp/ksp/tutorials/
H A Dex55.c12 Mat Amat; in main() local
49 PetscCall(MatCreate(comm, &Amat)); in main()
50 PetscCall(MatSetSizes(Amat, m, m, M, M)); in main()
51 PetscCall(MatSetType(Amat, MATAIJ)); in main()
52 PetscCall(MatSetOption(Amat, MAT_SPD, PETSC_TRUE)); in main()
53 PetscCall(MatSetOption(Amat, MAT_SPD_ETERNAL, PETSC_TRUE)); in main()
54 PetscCall(MatSetFromOptions(Amat)); in main()
55 PetscCall(MatSetBlockSize(Amat, 2)); in main()
56 PetscCall(MatSeqAIJSetPreallocation(Amat, 18, NULL)); in main()
57 PetscCall(MatMPIAIJSetPreallocation(Amat, 18, NULL, 18, NULL)); in main()
[all …]
H A Dex54.c9 Mat Amat, Pmat; in main() local
39 PetscCall(MatCreate(comm, &Amat)); in main()
40 PetscCall(MatSetSizes(Amat, PETSC_DECIDE, PETSC_DECIDE, M, M)); in main()
41 PetscCall(MatSetType(Amat, MATAIJ)); in main()
42 PetscCall(MatSetOption(Amat, MAT_SPD, PETSC_TRUE)); in main()
43 PetscCall(MatSetOption(Amat, MAT_SPD_ETERNAL, PETSC_TRUE)); in main()
44 PetscCall(MatSetFromOptions(Amat)); in main()
45 PetscCall(MatSeqAIJSetPreallocation(Amat, 81, NULL)); in main()
46 PetscCall(MatMPIAIJSetPreallocation(Amat, 81, NULL, 57, NULL)); in main()
56 PetscCall(MatCreateVecs(Amat, &bb, &xx)); in main()
[all …]
H A Dex56.c26 Mat Amat; in main() local
119 PetscCall(MatCreate(comm, &Amat)); in main()
120 PetscCall(MatSetSizes(Amat, m, m, M, M)); in main()
121 if (!test_late_bs) PetscCall(MatSetBlockSize(Amat, 3)); in main()
122 PetscCall(MatSetType(Amat, MATAIJ)); in main()
123 PetscCall(MatSetOption(Amat, MAT_SPD, PETSC_TRUE)); in main()
124 …PetscCall(MatSetOption(Amat, MAT_SPD_ETERNAL, PETSC_TRUE)); // this keeps CG after switch to negat… in main()
125 PetscCall(MatSetFromOptions(Amat)); in main()
126 PetscCall(MatSeqAIJSetPreallocation(Amat, 0, d_nnz)); in main()
127 PetscCall(MatMPIAIJSetPreallocation(Amat, 0, d_nnz, 0, o_nnz)); in main()
[all …]
H A Dex54f.F9024 Mat Amat
84 PetscCallA(MatCreate(PETSC_COMM_WORLD, Amat, ierr))
85 PetscCallA(MatSetSizes(Amat, PETSC_DECIDE, PETSC_DECIDE, M, M, ierr))
86 PetscCallA(MatSetType(Amat, MATAIJ, ierr))
87 PetscCallA(MatSetOption(Amat, MAT_SPD, PETSC_TRUE, ierr))
88 PetscCallA(MatSetOption(Amat, MAT_SPD_ETERNAL, PETSC_TRUE, ierr))
90 PetscCallA(MatSetType(Amat, MATAIJ, ierr))
92 PetscCallA(MatSetType(Amat, MATMPIAIJ, ierr))
94 …PetscCallA(MatMPIAIJSetPreallocation(Amat, f9, PETSC_NULL_INTEGER_ARRAY, f6, PETSC_NULL_INTEGER_AR…
95 PetscCallA(MatSetFromOptions(Amat, ierr))
[all …]
/petsc/src/ts/interface/sensitivity/
H A Dtssen.c28 PetscErrorCode TSSetRHSJacobianP(TS ts, Mat Amat, TSRHSJacobianPFn *func, PetscCtx ctx) in TSSetRHSJacobianP() argument
32 PetscValidHeaderSpecific(Amat, MAT_CLASSID, 2); in TSSetRHSJacobianP()
36 if (Amat) { in TSSetRHSJacobianP()
37 PetscCall(PetscObjectReference((PetscObject)Amat)); in TSSetRHSJacobianP()
39 ts->Jacprhs = Amat; in TSSetRHSJacobianP()
64 PetscErrorCode TSGetRHSJacobianP(TS ts, Mat *Amat, TSRHSJacobianPFn **func, PetscCtxRt ctx) in TSGetRHSJacobianP() argument
69 if (Amat) *Amat = ts->Jacprhs; in TSGetRHSJacobianP()
90 PetscErrorCode TSComputeRHSJacobianP(TS ts, PetscReal t, Vec U, Mat Amat) in TSComputeRHSJacobianP() argument
93 if (!Amat) PetscFunctionReturn(PETSC_SUCCESS); in TSComputeRHSJacobianP()
97 …llback JacobianP for sensitivity analysis", (*ts->rhsjacobianp)(ts, t, U, Amat, ts->rhsjacobianpct… in TSComputeRHSJacobianP()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipecgrr/
H A Dpipecgrr.c27 Mat Amat, Pmat; in KSPSolve_PIPECGRR() local
46 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_PIPECGRR()
50 PetscCall(KSP_MatMult(ksp, Amat, X, R)); /* r <- b - Ax */ in KSPSolve_PIPECGRR()
63 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECGRR()
71 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECGRR()
79 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECGRR()
86 PetscCall(KSP_MatMult(ksp, Amat, U, W)); in KSPSolve_PIPECGRR()
98 PetscCall(MatNorm(Amat, NORM_INFINITY, &Anorm)); in KSPSolve_PIPECGRR()
140 PetscCall(KSP_MatMult(ksp, Amat, M, N)); /* n <- Am */ in KSPSolve_PIPECGRR()
222 PetscCall(KSP_MatMult(ksp, Amat, X, R)); /* r <- Ax - b */ in KSPSolve_PIPECGRR()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipecg/
H A Dpipecg.c26 Mat Amat, Pmat; in KSPSolve_PIPECG() local
45 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_PIPECG()
49 PetscCall(KSP_MatMult(ksp, Amat, X, R)); /* r <- b - Ax */ in KSPSolve_PIPECG()
61 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECG()
67 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECG()
73 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECG()
79 PetscCall(KSP_MatMult(ksp, Amat, U, W)); in KSPSolve_PIPECG()
103 PetscCall(KSP_MatMult(ksp, Amat, M, N)); /* n <- Am */ in KSPSolve_PIPECG()
/petsc/src/dm/impls/stag/tests/
H A Dex40.c63 PetscErrorCode FormJacobian1DNoCoupling(SNES snes, Vec x, Mat Amat, Mat Pmat, PetscCtx ctx) in FormJacobian1DNoCoupling() argument
87 …PetscCall(DMStagMatSetValuesStencil(dm, Amat, 1, &row_vertex, 1, &row_vertex, &val, INSERT_VALUES)… in FormJacobian1DNoCoupling()
99 …PetscCall(DMStagMatSetValuesStencil(dm, Amat, 1, &row_element, 1, &row_element, &val, INSERT_VALUE… in FormJacobian1DNoCoupling()
105 PetscCall(MatAssemblyBegin(Amat, MAT_FINAL_ASSEMBLY)); in FormJacobian1DNoCoupling()
106 PetscCall(MatAssemblyEnd(Amat, MAT_FINAL_ASSEMBLY)); in FormJacobian1DNoCoupling()
107 …PetscCheck(Amat == Pmat, PETSC_COMM_WORLD, PETSC_ERR_SUP, "Not implemented for distinct Amat and P… in FormJacobian1DNoCoupling()
229 PetscErrorCode FormJacobian1D(SNES snes, Vec x, Mat Amat, Mat Pmat, PetscCtx ctx) in FormJacobian1D() argument
248 PetscCall(MatZeroEntries(Amat)); in FormJacobian1D()
275 … PetscCall(DMStagMatSetValuesStencil(dm, Amat, 1, &row_vertex, 1, &col, &val, ADD_VALUES)); in FormJacobian1D()
293 … PetscCall(DMStagMatSetValuesStencil(dm, Amat, 1, &row_vertex, 1, &col, &val, ADD_VALUES)); in FormJacobian1D()
[all …]
/petsc/src/ksp/ksp/impls/cr/pipecr/
H A Dpipecr.c26 Mat Amat, Pmat; in KSPSolve_PIPECR() local
43 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_PIPECR()
48 PetscCall(KSP_MatMult(ksp, Amat, X, W)); /* w <- b - Ax */ in KSPSolve_PIPECR()
59 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECR()
63 PetscCall(KSP_MatMult(ksp, Amat, U, W)); in KSPSolve_PIPECR()
84 PetscCall(KSP_MatMult(ksp, Amat, M, N)); /* n <- Am */ in KSPSolve_PIPECR()
/petsc/src/ksp/ksp/impls/cg/pipeprcg/
H A Dpipeprcg.c45 Mat Amat, Pmat; in KSPSolve_PIPEPRCG() local
69 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_PIPEPRCG()
74 PetscCall(KSP_MatMult(ksp, Amat, X, R)); /* r <- b - Ax */ in KSPSolve_PIPEPRCG()
81 PetscCall(KSP_MatMult(ksp, Amat, RT, W)); /* w <- A rt */ in KSPSolve_PIPEPRCG()
88 PetscCall(KSP_MatMult(ksp, Amat, ST, U)); /* u <- Ast */ in KSPSolve_PIPEPRCG()
156 PetscCall(KSP_MatMult(ksp, Amat, ST, U)); /* u <- A st */ in KSPSolve_PIPEPRCG()
162 PetscCall(KSP_MatMult(ksp, Amat, RT, W)); /* w <- A rt */ in KSPSolve_PIPEPRCG()
/petsc/src/ksp/ksp/impls/cg/cgne/
H A Dcgne.c48 Mat Amat, Pmat; in KSPSolve_CGNE() local
73 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_CGNE()
76 PetscCall(KSP_MatMultTranspose(ksp, Amat, B, T)); in KSPSolve_CGNE()
78 PetscCall(KSP_MatMult(ksp, Amat, X, P)); in KSPSolve_CGNE()
79 PetscCall(KSP_MatMultTranspose(ksp, Amat, P, R)); in KSPSolve_CGNE()
134 PetscCall(KSP_MatMult(ksp, Amat, P, T)); in KSPSolve_CGNE()
135 PetscCall(KSP_MatMultTranspose(ksp, Amat, T, Z)); in KSPSolve_CGNE()
/petsc/src/ksp/ksp/impls/bicg/
H A Dbicg.c17 Mat Amat, Pmat; in KSPSolve_BiCG() local
32 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_BiCG()
35 PetscCall(KSP_MatMult(ksp, Amat, X, Rr)); /* r <- b - Ax */ in KSPSolve_BiCG()
77 PetscCall(KSP_MatMult(ksp, Amat, Pr, Zr)); /* z <- Kp */ in KSPSolve_BiCG()
78 PetscCall(KSP_MatMultHermitianTranspose(ksp, Amat, Pl, Zl)); in KSPSolve_BiCG()
/petsc/src/ts/tests/
H A Dex3.c23 Mat Amat; /* left hand side matrix */ member
91 PetscCall(MatCreateSeqAIJ(PETSC_COMM_WORLD, m, m, 3, NULL, &appctx.Amat)); in main()
92 PetscCall(MatSetFromOptions(appctx.Amat)); in main()
93 PetscCall(MatSetUp(appctx.Amat)); in main()
141 PetscCall(TSSetIJacobian(ts, appctx.Amat, appctx.Amat, TSComputeIJacobianConstant, &appctx)); in main()
183 PetscCall(MatDestroy(&appctx.Amat)); in main()
276 PetscCall(KSPSetOperators(ksp, obj->Amat, obj->Amat)); in Petsc_KSPSolve()
480 PetscCall(MatSetValue(obj->Amat, i, j, add_term, ADD_VALUES)); in femA()
487 PetscCall(MatAssemblyBegin(obj->Amat, MAT_FINAL_ASSEMBLY)); in femA()
488 PetscCall(MatAssemblyEnd(obj->Amat, MAT_FINAL_ASSEMBLY)); in femA()
/petsc/src/snes/tutorials/
H A Dex73f90t.F9056 Mat::Amat, AmatLin, Bmat, CMat, Dmat
195 Mat:: Amat
203 …PetscCall(MatCreateSubMatrix(jac_prec, solver%isPhi, solver%isPhi, MAT_INITIAL_MATRIX, Amat, ierr))
205 PetscCall(FormJacobianLocal(Xsub(1), Amat, solver, .true., ierr))
206 PetscCall(MatDestroy(Amat, ierr)) ! discard our reference
424 Mat:: Amat, Bmat, Cmat, Dmat, KKTMat, matArray(4)
489 PetscCallA(MatCreate(PETSC_COMM_WORLD, Amat, ierr))
490 PetscCallA(MatSetSizes(Amat, PETSC_DECIDE, PETSC_DECIDE, N1, N1, ierr))
491 PetscCallA(MatSetUp(Amat, ierr))
502 PetscCallA(DMShellSetMatrix(daphi, Amat, ierr))
[all …]
H A Dex56.c194 Mat Amat; in main() local
349 PetscCall(DMCreateMatrix(dm, &Amat)); in main()
350 …PetscCall(MatSetOption(Amat, MAT_SYMMETRIC, PETSC_TRUE)); /* Some matrix kernels can take a… in main()
351 …PetscCall(MatSetOption(Amat, MAT_SYMMETRY_ETERNAL, PETSC_TRUE)); /* Inform PETSc that Amat is alwa… in main()
352 PetscCall(MatSetBlockSize(Amat, Ncomp)); in main()
353 PetscCall(MatSetOption(Amat, MAT_SPD, PETSC_TRUE)); in main()
354 PetscCall(MatSetOption(Amat, MAT_SPD_ETERNAL, PETSC_TRUE)); in main()
369 if (attach_nearnullspace) PetscCall(MatSetNearNullSpace(Amat, nearNullSpace)); in main()
373 PetscCall(SNESSetJacobian(snes, Amat, Amat, NULL, NULL)); in main()
399 PetscCall(SNESComputeJacobian(snes, xx, Amat, Amat)); in main()
[all …]
/petsc/src/ksp/ksp/impls/cg/groppcg/
H A Dgroppcg.c29 Mat Amat, Pmat; in KSPSolve_GROPPCG() local
45 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_GROPPCG()
49 PetscCall(KSP_MatMult(ksp, Amat, x, r)); /* r <- b - Ax */ in KSPSolve_GROPPCG()
59 PetscCall(KSP_MatMult(ksp, Amat, p, s)); /* s <- Ap */ in KSPSolve_GROPPCG()
112 PetscCall(KSP_MatMult(ksp, Amat, z, Z)); /* Z <- Az */ in KSPSolve_GROPPCG()
/petsc/src/ts/interface/
H A Dtsrhssplit.c216 PetscErrorCode TSRHSSplitSetIJacobian(TS ts, const char splitname[], Mat Amat, Mat Pmat, TSIJacobia… in TSRHSSplitSetIJacobian() argument
223 if (Amat) PetscValidHeaderSpecific(Amat, MAT_CLASSID, 3); in TSRHSSplitSetIJacobian()
225 if (Amat) PetscCheckSameComm(ts, 1, Amat, 3); in TSRHSSplitSetIJacobian()
243 PetscCall(TSSetIJacobian(isplit->ts, Amat, Pmat, ijac, ctx)); in TSRHSSplitSetIJacobian()
/petsc/src/ksp/ksp/interface/
H A Ditres.c37 Mat Amat, Pmat; in KSPInitialResidual() local
45 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPInitialResidual()
48 PetscCall(KSP_MatMult(ksp, Amat, vsoln, vt1)); in KSPInitialResidual()
H A Ditcreate.c539 PetscErrorCode KSPSetOperators(KSP ksp, Mat Amat, Mat Pmat) in KSPSetOperators() argument
543 if (Amat) PetscValidHeaderSpecific(Amat, MAT_CLASSID, 2); in KSPSetOperators()
545 if (Amat) PetscCheckSameComm(ksp, 1, Amat, 2); in KSPSetOperators()
548 PetscCall(PCSetOperators(ksp->pc, Amat, Pmat)); in KSPSetOperators()
577 PetscErrorCode KSPGetOperators(KSP ksp, Mat *Amat, Mat *Pmat) in KSPGetOperators() argument
582 PetscCall(PCGetOperators(ksp->pc, Amat, Pmat)); in KSPGetOperators()
/petsc/src/ksp/ksp/impls/cr/
H A Dcr.c17 Mat Amat, Pmat; in KSPSolve_CR() local
30 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_CR()
32 PetscCall(KSP_MatMult(ksp, Amat, X, R)); /* R <- A*X */ in KSPSolve_CR()
40 PetscCall(KSP_MatMult(ksp, Amat, P, AP)); /* AP <- A*P */ in KSPSolve_CR()
92 PetscCall(KSP_MatMult(ksp, Amat, RT, ART)); /* ART <- A*RT */ in KSPSolve_CR()
/petsc/src/ksp/ksp/impls/cg/
H A Dcg.c127 Mat Amat, Pmat; in KSPSolve_CG() local
150 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_CG()
154 PetscCall(KSP_MatMult(ksp, Amat, X, R)); /* r <- b - Ax */ in KSPSolve_CG()
257 PetscCall(KSP_MatMult(ksp, Amat, P, W)); /* w <- Ap */ in KSPSolve_CG()
371 Mat Amat, Pmat; in KSPSolve_CG_SingleReduction() local
394 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_CG_SingleReduction()
398 PetscCall(KSP_MatMult(ksp, Amat, X, R)); /* r <- b - Ax */ in KSPSolve_CG_SingleReduction()
416 PetscCall(KSP_MatMult(ksp, Amat, Z, S)); in KSPSolve_CG_SingleReduction()
437 PetscCall(KSP_MatMult(ksp, Amat, Z, S)); in KSPSolve_CG_SingleReduction()
471 PetscCall(KSP_MatMult(ksp, Amat, P, W)); /* w <- Ap */ in KSPSolve_CG_SingleReduction()
[all …]
/petsc/src/ksp/ksp/impls/lcd/
H A Dlcd.c35 Mat Amat, Pmat; in KSPSolve_LCD() local
50 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_LCD()
54 PetscCall(KSP_MatMult(ksp, Amat, X, Z)); /* z <- b - Ax */ in KSPSolve_LCD()
77 PetscCall(KSP_MatMult(ksp, Amat, lcd->P[it], Z)); in KSPSolve_LCD()
101 PetscCall(KSP_MatMult(ksp, Amat, lcd->P[it + 1], Z)); in KSPSolve_LCD()
/petsc/src/ksp/ksp/impls/cheby/
H A Dcheby.c395 Mat Amat, Pmat; in KSPSolve_Chebyshev_FirstKind() local
402 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_Chebyshev_FirstKind()
432 PetscCall(KSP_MatMult(ksp, Amat, sol_orig, r)); /* r = b - A*p[km1] */ in KSPSolve_Chebyshev_FirstKind()
475 PetscCall(KSP_MatMult(ksp, Amat, p[k], r)); /* r = b - Ap[k] */ in KSPSolve_Chebyshev_FirstKind()
520 PetscCall(KSP_MatMult(ksp, Amat, p[k], r)); /* r = b - Ap[k] */ in KSPSolve_Chebyshev_FirstKind()
564 Mat Amat, Pmat; in KSPSolve_Chebyshev_FourthKind() local
572 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_Chebyshev_FourthKind()
588 PetscCall(KSP_MatMult(ksp, Amat, x, r)); /* r = b - A*x */ in KSPSolve_Chebyshev_FourthKind()
633 PetscCall(KSP_MatMult(ksp, Amat, d, Br)); /* r = r - Ad */ in KSPSolve_Chebyshev_FourthKind()
677 PetscCall(KSP_MatMult(ksp, Amat, x, r)); /* r = b - Ax */ in KSPSolve_Chebyshev_FourthKind()
[all …]
/petsc/src/ksp/ksp/utils/lmvm/dense/
H A Dcd_utils.c143 PETSC_INTERN PetscErrorCode MatUpperTriangularSolveInPlace(Mat B, Mat Amat, Vec X, PetscBool hermit… in MatUpperTriangularSolveInPlace() argument
158 PetscCall(MatDenseGetArrayReadAndMemType(Amat, &A, &memtype_r)); in MatUpperTriangularSolveInPlace()
160 PetscCall(MatDenseRestoreArrayReadAndMemType(Amat, &A)); in MatUpperTriangularSolveInPlace()
165 PetscCall(MatDenseGetLDA(Amat, &lda)); in MatUpperTriangularSolveInPlace()
168 PetscCall(MatDenseRestoreArrayReadAndMemType(Amat, &A)); in MatUpperTriangularSolveInPlace()
/petsc/src/ksp/ksp/impls/lsqr/
H A Dlsqr.c74 Mat Amat, Pmat; in KSPSolve_LSQR() local
82 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_LSQR()
103 PetscCall(KSP_MatMult(ksp, Amat, X, U)); /* u <- b - Ax */ in KSPSolve_LSQR()
123 PetscCall(KSP_MatMultHermitianTranspose(ksp, Amat, U, V)); in KSPSolve_LSQR()
148 if (lsqr->exact_norm) PetscCall(MatNorm(Amat, NORM_FROBENIUS, &lsqr->anorm)); in KSPSolve_LSQR()
157 PetscCall(KSP_MatMult(ksp, Amat, V, U1)); in KSPSolve_LSQR()
159 PetscCall(KSP_MatMult(ksp, Amat, Z, U1)); in KSPSolve_LSQR()
169 PetscCall(KSP_MatMultHermitianTranspose(ksp, Amat, U1, V1)); in KSPSolve_LSQR()
235 PetscCall(MatGetSize(Amat, &size1, &size2)); in KSPSolve_LSQR()

1234