| /petsc/src/ts/tutorials/hybrid/ |
| H A D | ex1adj.c | 197 Mat Ap; /* dfdp */ in main() local
238 PetscCall(MatCreate(PETSC_COMM_WORLD, &Ap)); in main()
239 PetscCall(MatSetSizes(Ap, n, 1, PETSC_DETERMINE, PETSC_DETERMINE)); in main()
240 PetscCall(MatSetType(Ap, MATDENSE)); in main()
241 PetscCall(MatSetFromOptions(Ap)); in main()
242 PetscCall(MatSetUp(Ap)); in main()
243 PetscCall(MatZeroEntries(Ap)); /* initialize to zeros */ in main()
257 PetscCall(TSSetRHSJacobianP(ts, Ap, RHSJacobianP, &app)); in main()
302 PetscCall(MatCreateVecs(Ap, &mu[0], NULL)); in main()
303 PetscCall(MatCreateVecs(Ap, &mu[1], NULL)); in main()
[all …]
|
| H A D | ex1fwd.c | 209 static PetscErrorCode RHSJacobianP(TS ts, PetscReal t, Vec X, Mat Ap, PetscCtx ctx) in RHSJacobianP() argument
220 Mat Ap; /* Jacobian dfdp */ in main() local
260 PetscCall(MatCreate(PETSC_COMM_WORLD, &Ap)); in main()
261 PetscCall(MatSetSizes(Ap, n, 3, PETSC_DETERMINE, PETSC_DETERMINE)); in main()
262 PetscCall(MatSetType(Ap, MATDENSE)); in main()
263 PetscCall(MatSetFromOptions(Ap)); in main()
264 PetscCall(MatSetUp(Ap)); in main()
265 PetscCall(MatZeroEntries(Ap)); in main()
292 PetscCall(TSSetRHSJacobianP(ts, Ap, RHSJacobianP, &app)); in main()
319 PetscCall(MatDestroy(&Ap)); in main()
|
| H A D | ex1fd.c | 190 Mat Ap; /* dfdp */ in main() local
230 PetscCall(MatCreate(PETSC_COMM_WORLD, &Ap)); in main()
231 PetscCall(MatSetSizes(Ap, n, 1, PETSC_DETERMINE, PETSC_DETERMINE)); in main()
232 PetscCall(MatSetType(Ap, MATDENSE)); in main()
233 PetscCall(MatSetFromOptions(Ap)); in main()
234 PetscCall(MatSetUp(Ap)); in main()
235 PetscCall(MatZeroEntries(Ap)); /* initialize to zeros */ in main()
317 PetscCall(MatDestroy(&Ap)); in main()
|
| /petsc/src/ksp/pc/impls/galerkin/ |
| H A D | galerkin.c | 43 Mat Ap; in PCSetUp_Galerkin() local
45 PetscCall((*jac->computeasub)(pc, pc->pmat, NULL, &Ap, jac->computeasub_ctx)); in PCSetUp_Galerkin()
46 PetscCall(KSPSetOperators(jac->ksp, Ap, Ap)); in PCSetUp_Galerkin()
47 PetscCall(MatDestroy(&Ap)); in PCSetUp_Galerkin()
49 PetscCall(KSPGetOperators(jac->ksp, NULL, &Ap)); in PCSetUp_Galerkin()
50 PetscCall((*jac->computeasub)(pc, pc->pmat, Ap, NULL, jac->computeasub_ctx)); in PCSetUp_Galerkin()
231 …SetComputeSubmatrix(PC pc, PetscErrorCode (*computeAsub)(PC pc, Mat A, Mat Ap, Mat *cAp, PetscCtx … in PCGalerkinSetComputeSubmatrix() argument
|
| /petsc/src/ksp/pc/impls/pbjacobi/kokkos/ |
| H A D | pbjacobi_kok.kokkos.cxx | 49 const PetscScalar *Ap, *xp; in PCApplyOrTranspose_PBJacobi_Kokkos()
55 Ap = &Av(bs2 * k + i * (transpose ? bs : 1)); /* Ap points to the first entry of i-th row */ in PCApplyOrTranspose_PBJacobi_Kokkos()
61 yp[i] += Ap[0] * xp[j]; in PCApplyOrTranspose_PBJacobi_Kokkos()
62 Ap += (transpose ? 1 : bs); /* block is in column major order */ in PCApplyOrTranspose_PBJacobi_Kokkos()
|
| /petsc/src/ksp/pc/impls/pbjacobi/cuda/ |
| H A D | pbjacobi_cuda.cu | 14 const PetscScalar *Ap, *xp; in MatMultBatched() local
20 Ap = &A[bs2 * k + i * (transpose ? bs : 1)]; /* Ap points to the first entry of i-th row */ in MatMultBatched()
26 yp[i] += Ap[0] * xp[j]; in MatMultBatched()
27 Ap += (transpose ? 1 : bs); /* block is in column major order */ in MatMultBatched()
|
| /petsc/src/ksp/ksp/utils/schurm/ |
| H A D | schurm.c | 124 PetscCall(MatDestroy(&Na->Ap)); in MatDestroy_SchurComplement()
251 Na->Ap = Ap00; in MatSchurComplementSetSubMatrices()
325 PetscCall(KSPSetOperators(Na->ksp, Na->A, Na->Ap)); in MatSchurComplementSetKSP()
394 PetscCall(MatDestroy(&Na->Ap)); in MatSchurComplementUpdateSubMatrices()
400 Na->Ap = Ap00; in MatSchurComplementUpdateSubMatrices()
443 if (Ap00) *Ap00 = Na->Ap; in MatSchurComplementGetSubMatrices()
603 Mat A = NULL, Ap = NULL, B = NULL, C = NULL, D = NULL; in MatGetSchurComplement_Basic() local
624 PetscCall(MatSchurComplementGetSubMatrices(*S, &A, &Ap, &B, &C, &D)); in MatGetSchurComplement_Basic()
625 …PetscCheck(A && Ap && B && C, PetscObjectComm((PetscObject)mat), PETSC_ERR_ARG_WRONGSTATE, "Attemp… in MatGetSchurComplement_Basic()
626 …PetscCheck(A == Ap, PetscObjectComm((PetscObject)mat), PETSC_ERR_ARG_WRONGSTATE, "Matrix for const… in MatGetSchurComplement_Basic()
[all …]
|
| H A D | schurm.h | 8 Mat A, Ap, B, C, D; member
|
| /petsc/src/ksp/ksp/impls/fetidp/ |
| H A D | fetidp.c | 511 Mat A, Ap; in KSPFETIDPSetUpOperators() local
536 PetscCall(KSPGetOperators(ksp, &A, &Ap)); in KSPFETIDPSetUpOperators()
930 PetscCall(MatGetNearNullSpace(Ap, &nnsp)); in KSPFETIDPSetUpOperators()
931 if (!nnsp) PetscCall(MatGetNullSpace(Ap, &nnsp)); in KSPFETIDPSetUpOperators()
1004 …if (A != Ap) { /* user has provided a different Pmat, this always supersedes the setter (TODO: is … in KSPFETIDPSetUpOperators()
1006 PetscCall(KSPFETIDPSetPressureOperator(ksp, Ap)); in KSPFETIDPSetUpOperators()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex59.c | 790 Mat A, Ap; in ComputeKSPFETIDP() local
793 PetscCall(KSPGetOperators(ksp_bddc, &A, &Ap)); in ComputeKSPFETIDP()
794 PetscCall(KSPSetOperators(temp_ksp, A, Ap)); in ComputeKSPFETIDP()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexcreate.c | 4011 …6, R = 17, S = 18, T = 19, U = 20, V = 21, W = 22, X = 23, Y = 24, Z = 25, Ap = 26, Bp = 27, Cp = … in DMPlexCreateTPSMesh_Internal() local
4035 {R, O + 56 * 1, U + 56 * 1, Ap }, in DMPlexCreateTPSMesh_Internal()
4036 {P, R, Ap, W }, in DMPlexCreateTPSMesh_Internal()
4044 {Ap, U + 56 * 1, Dp + 56 * 1, Gp }, in DMPlexCreateTPSMesh_Internal()
4045 {W, Ap, Gp, Cp }, in DMPlexCreateTPSMesh_Internal()
|
| /petsc/doc/manual/ |
| H A D | ksp.md | 2103 For notational simplicity define $\text{ksp}(A,Ap)$ to mean
2105 $A$ and preconditioner built from matrix $Ap$.
|