Home
last modified time | relevance | path

Searched refs:Pmat (Results 1 – 25 of 57) sorted by relevance

123

/petsc/src/ksp/ksp/tutorials/
H A Dex19.c16 Mat A, Pmat, Aseq, AA; /* linear system matrix */ in main() local
105 PetscCall(MatCreate(PETSC_COMM_WORLD, &Pmat)); in main()
106 PetscCall(MatSetSizes(Pmat, Iend - Istart, Iend - Istart, n * m, n * m)); in main()
107 PetscCall(MatSetFromOptions(Pmat)); in main()
108 PetscCall(MatSeqAIJSetPreallocation(Pmat, 5, NULL)); in main()
109 PetscCall(MatMPIAIJSetPreallocation(Pmat, 5, NULL, 3, NULL)); in main()
110 PetscCall(MatAssemblyBegin(Pmat, MAT_FINAL_ASSEMBLY)); in main()
111 PetscCall(MatAssemblyEnd(Pmat, MAT_FINAL_ASSEMBLY)); in main()
112 PetscCall(MatMPIAIJGetSeqAIJ(Pmat, &AAseq, NULL, NULL)); in main()
117 Pmat = AA; in main()
[all …]
H A Dex54.c9 Mat Amat, Pmat; in main() local
48 PetscCall(MatCreate(comm, &Pmat)); in main()
49 PetscCall(MatSetSizes(Pmat, PETSC_DECIDE, PETSC_DECIDE, M, M)); in main()
50 PetscCall(MatSetType(Pmat, MATMPIAIJ)); in main()
51 PetscCall(MatSetFromOptions(Pmat)); in main()
52 PetscCall(MatSeqAIJSetPreallocation(Pmat, 81, NULL)); in main()
53 PetscCall(MatMPIAIJSetPreallocation(Pmat, 81, NULL, 57, NULL)); in main()
115 PetscCall(MatSetValues(Pmat, 4, idx, 4, idx, (const PetscScalar *)DD, ADD_VALUES)); in main()
134 PetscCall(MatAssemblyBegin(Pmat, MAT_FINAL_ASSEMBLY)); in main()
135 PetscCall(MatAssemblyEnd(Pmat, MAT_FINAL_ASSEMBLY)); in main()
[all …]
/petsc/src/ts/tutorials/
H A Dextchemfield.c275 static PetscErrorCode FormDiffusionJacobian(TS ts, PetscReal t, Vec X, Mat Amat, Mat Pmat, void *pt… in FormDiffusionJacobian() argument
302 PetscCall(MatSetValuesStencil(Pmat, 1, &row, 3, col, values, ADD_VALUES)); in FormDiffusionJacobian()
305 PetscCall(MatAssemblyBegin(Pmat, MAT_FINAL_ASSEMBLY)); in FormDiffusionJacobian()
306 PetscCall(MatAssemblyEnd(Pmat, MAT_FINAL_ASSEMBLY)); in FormDiffusionJacobian()
341 static PetscErrorCode FormRHSJacobian(TS ts, PetscReal t, Vec X, Mat Amat, Mat Pmat, void *ptr) in FormRHSJacobian() argument
351 PetscCall(MatZeroEntries(Pmat)); in FormRHSJacobian()
352 PetscCall(MatSetOption(Pmat, MAT_ROW_ORIENTED, PETSC_FALSE)); in FormRHSJacobian()
353 PetscCall(MatSetOption(Pmat, MAT_IGNORE_ZERO_ENTRIES, PETSC_TRUE)); in FormRHSJacobian()
365 PetscCall(MatSetValues(Pmat, M, user->rows, M, user->rows, user->Jdense, INSERT_VALUES)); in FormRHSJacobian()
368 PetscCall(MatAssemblyBegin(Pmat, MAT_FINAL_ASSEMBLY)); in FormRHSJacobian()
[all …]
H A Dextchem.c247 static PetscErrorCode FormRHSJacobian(TS ts, PetscReal t, Vec X, Mat Amat, Mat Pmat, void *ptr) in FormRHSJacobian() argument
263 PetscCall(MatSetOption(Pmat, MAT_ROW_ORIENTED, PETSC_FALSE)); in FormRHSJacobian()
264 PetscCall(MatSetOption(Pmat, MAT_IGNORE_ZERO_ENTRIES, PETSC_TRUE)); in FormRHSJacobian()
265 PetscCall(MatZeroEntries(Pmat)); in FormRHSJacobian()
266 PetscCall(MatSetValues(Pmat, M, user->rows, M, user->rows, user->Jdense, INSERT_VALUES)); in FormRHSJacobian()
268 PetscCall(MatAssemblyBegin(Pmat, MAT_FINAL_ASSEMBLY)); in FormRHSJacobian()
269 PetscCall(MatAssemblyEnd(Pmat, MAT_FINAL_ASSEMBLY)); in FormRHSJacobian()
270 if (Amat != Pmat) { in FormRHSJacobian()
/petsc/src/ts/interface/
H A Dtsrhssplit.c216 PetscErrorCode TSRHSSplitSetIJacobian(TS ts, const char splitname[], Mat Amat, Mat Pmat, TSIJacobia… in TSRHSSplitSetIJacobian() argument
224 if (Pmat) PetscValidHeaderSpecific(Pmat, MAT_CLASSID, 4); in TSRHSSplitSetIJacobian()
226 if (Pmat) PetscCheckSameComm(ts, 1, Pmat, 4); in TSRHSSplitSetIJacobian()
243 PetscCall(TSSetIJacobian(isplit->ts, Amat, Pmat, ijac, ctx)); in TSRHSSplitSetIJacobian()
H A Dts.c1143 PetscErrorCode TSSetRHSJacobian(TS ts, Mat Amat, Mat Pmat, TSRHSJacobianFn *f, PetscCtx ctx) in TSSetRHSJacobian() argument
1152 if (Pmat) PetscValidHeaderSpecific(Pmat, MAT_CLASSID, 3); in TSSetRHSJacobian()
1154 if (Pmat) PetscCheckSameComm(ts, 1, Pmat, 3); in TSSetRHSJacobian()
1160 if (!ijacobian) PetscCall(SNESSetJacobian(snes, Amat, Pmat, SNESTSFormJacobian, ts)); in TSSetRHSJacobian()
1166 if (Pmat) { in TSSetRHSJacobian()
1167 PetscCall(PetscObjectReference((PetscObject)Pmat)); in TSSetRHSJacobian()
1169 ts->Brhs = Pmat; in TSSetRHSJacobian()
1317 PetscErrorCode TSSetIJacobian(TS ts, Mat Amat, Mat Pmat, TSIJacobianFn *f, PetscCtx ctx) in TSSetIJacobian() argument
1325 if (Pmat) PetscValidHeaderSpecific(Pmat, MAT_CLASSID, 3); in TSSetIJacobian()
1327 if (Pmat) PetscCheckSameComm(ts, 1, Pmat, 3); in TSSetIJacobian()
[all …]
/petsc/src/ksp/ksp/interface/
H A Ditcreate.c539 PetscErrorCode KSPSetOperators(KSP ksp, Mat Amat, Mat Pmat) in KSPSetOperators() argument
544 if (Pmat) PetscValidHeaderSpecific(Pmat, MAT_CLASSID, 3); in KSPSetOperators()
546 if (Pmat) PetscCheckSameComm(ksp, 1, Pmat, 3); 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()
H A Ditres.c37 Mat Amat, Pmat; in KSPInitialResidual() local
45 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPInitialResidual()
/petsc/src/ksp/pc/interface/
H A Dprecon.c1269 PetscErrorCode PCSetOperators(PC pc, Mat Amat, Mat Pmat) in PCSetOperators() argument
1276 if (Pmat) PetscValidHeaderSpecific(Pmat, MAT_CLASSID, 3); in PCSetOperators()
1278 if (Pmat) PetscCheckSameComm(pc, 1, Pmat, 3); in PCSetOperators()
1279 if (pc->setupcalled && pc->mat && pc->pmat && Amat && Pmat) { in PCSetOperators()
1283 PetscCall(MatGetLocalSize(Pmat, &m1, &n1)); in PCSetOperators()
1288 if (Pmat != pc->pmat) { in PCSetOperators()
1297 if (Pmat) PetscCall(PetscObjectReference((PetscObject)Pmat)); in PCSetOperators()
1300 pc->pmat = Pmat; in PCSetOperators()
1416 PetscErrorCode PCGetOperators(PC pc, Mat *Amat, Mat *Pmat) in PCGetOperators() argument
1422 if (pc->pmat && !Pmat) { /* Pmat has been set, but user did not request it, so use for Amat */ in PCGetOperators()
[all …]
/petsc/src/ksp/ksp/tutorials/output/
H A Dex43_2_mumps.out18 using Amat (not Pmat) as operator for blocks
19 using Amat (not Pmat) as operator for diagonal blocks
20 using Amat (not Pmat) as operator for off-diagonal blocks
/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()
564 Mat Amat, Pmat; in KSPSolve_Chebyshev_FourthKind() local
572 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_Chebyshev_FourthKind()
752 Mat Pmat, Amat; in KSPSetUp_Chebyshev() local
793 PetscCall(KSPGetOperators(ksp, &Amat, &Pmat)); in KSPSetUp_Chebyshev()
794 PetscCall(MatIsSPDKnown(Pmat, &isset, &flg)); in KSPSetUp_Chebyshev()
803 PetscCall(PetscObjectGetId((PetscObject)Pmat, &pmatid)); in KSPSetUp_Chebyshev()
805 PetscCall(PetscObjectStateGet((PetscObject)Pmat, &pmatstate)); in KSPSetUp_Chebyshev()
/petsc/src/ksp/pc/impls/amgx/
H A Damgx.cxx218 Mat Pmat = pc->pmat; in PCSetUp_AMGX() local
222 …PetscCall(PetscObjectTypeCompareAny((PetscObject)Pmat, &is_dev_ptrs, MATAIJCUSPARSE, MATSEQAIJCUSP… in PCSetUp_AMGX()
229 PetscCall(MatMPIAIJGetLocalMat(Pmat, MAT_REUSE_MATRIX, &amgx->localA)); in PCSetUp_AMGX()
231 PetscCall(MatMPIAIJGetLocalMat(Pmat, MAT_INITIAL_MATRIX, &amgx->localA)); in PCSetUp_AMGX()
236 amgx->localA = Pmat; in PCSetUp_AMGX()
288 PetscCall(MatGetBlockSize(Pmat, &amgx->bSize)); in PCSetUp_AMGX()
/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()
/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()
/petsc/src/ksp/ksp/impls/rich/
H A Drich.c26 Mat Amat, Pmat; in KSPSolve_Richardson() local
37 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_Richardson()
59 PetscCall(MatGetNullSpace(Pmat, &nullsp)); in KSPSolve_Richardson()
/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()
/petsc/src/ksp/ksp/impls/gmres/fgmres/
H A Dfgmres.c44 Mat Amat, Pmat; in KSPFGMRESResidual() local
47 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPFGMRESResidual()
64 Mat Amat, Pmat; in KSPFGMRESCycle() local
125 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPFGMRESCycle()
/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()
/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()
/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()
/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()
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A Dpetscts.pxi395 cdef Mat Pmat = ref_Mat(P)
400 jacobian(Ts, toReal(t), Xvec, Jmat, Pmat, *args, **kargs)
455 cdef Mat Pmat = ref_Mat(P)
460 jacobian(Ts, toReal(t), Xvec, XDvec, toReal(a), Jmat, Pmat, *args, **kargs)
521 cdef Mat Pmat = ref_Mat(P)
526 jacobian(Ts, toReal(t), Xvec, XDvec, XDDvec, toReal(v), toReal(a), Jmat, Pmat, *args, **kargs)
/petsc/src/ksp/ksp/impls/symmlq/
H A Dsymmlq.c22 Mat Amat, Pmat; in KSPSolve_SYMMLQ() local
41 PetscCall(PCGetOperators(ksp->pc, &Amat, &Pmat)); in KSPSolve_SYMMLQ()
/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()
/petsc/src/ksp/pc/impls/gamg/
H A Dgamg.c542 Mat Pmat = pc->pmat; in PCSetUp_GAMG() local
625 PetscCall(MatGetBlockSize(Pmat, &bs)); in PCSetUp_GAMG()
626 PetscCall(MatGetLocalSize(Pmat, &qq, NULL)); in PCSetUp_GAMG()
636 PetscCall(pc_gamg->ops->createdefaultdata(pc, Pmat)); in PCSetUp_GAMG()
649 PetscCall(MatGetBlockSize(Pmat, &bs)); in PCSetUp_GAMG()
650 PetscCall(MatGetSize(Pmat, &M, NULL)); in PCSetUp_GAMG()
653 PetscCall(MatGetInfo(Pmat, MAT_GLOBAL_SUM, &info)); /* global reduction */ in PCSetUp_GAMG()
661 …for (level = 0, Aarr[0] = Pmat, nactivepe = size; level < (pc_gamg->Nlevels - 1) && (level == 0 ||… in PCSetUp_GAMG()
678 PetscCall(MatGetOwnershipRange(Pmat, &Istart, &Iend)); in PCSetUp_GAMG()
679 PetscCall(MatGetLocalSize(Pmat, &prol_m, NULL)); // rows m x n in PCSetUp_GAMG()
[all …]

123