| /petsc/src/mat/impls/submat/ |
| H A D | submat.c | 13 Mat_SubVirtual *Na = (Mat_SubVirtual *)N->data; in MatScale_SubMatrix() local
16 PetscCall(MatScale(Na->A, a)); in MatScale_SubMatrix()
22 Mat_SubVirtual *Na = (Mat_SubVirtual *)N->data; in MatShift_SubMatrix() local
25 PetscCall(MatShift(Na->A, a)); in MatShift_SubMatrix()
31 Mat_SubVirtual *Na = (Mat_SubVirtual *)N->data; in MatDiagonalScale_SubMatrix() local
35 PetscCall(VecZeroEntries(Na->rwork)); in MatDiagonalScale_SubMatrix()
36 PetscCall(VecScatterBegin(Na->rprolong, right, Na->rwork, INSERT_VALUES, SCATTER_FORWARD)); in MatDiagonalScale_SubMatrix()
37 PetscCall(VecScatterEnd(Na->rprolong, right, Na->rwork, INSERT_VALUES, SCATTER_FORWARD)); in MatDiagonalScale_SubMatrix()
40 PetscCall(VecZeroEntries(Na->lwork)); in MatDiagonalScale_SubMatrix()
41 PetscCall(VecScatterBegin(Na->lrestrict, left, Na->lwork, INSERT_VALUES, SCATTER_REVERSE)); in MatDiagonalScale_SubMatrix()
[all …]
|
| /petsc/src/mat/impls/lrc/ |
| H A D | lrc.c | 15 Mat_LRC *Na = (Mat_LRC *)N->data; in MatMult_LRC_kernel() local
20 U = transpose ? Na->V : Na->U; in MatMult_LRC_kernel()
21 V = transpose ? Na->U : Na->V; in MatMult_LRC_kernel()
24 PetscCall(MatMultHermitianTranspose(V, x, Na->work1)); in MatMult_LRC_kernel()
25 if (Na->c) PetscCall(VecPointwiseMult(Na->work1, Na->c, Na->work1)); in MatMult_LRC_kernel()
26 if (Na->A) { in MatMult_LRC_kernel()
28 PetscCall(MatMultTranspose(Na->A, x, y)); in MatMult_LRC_kernel()
30 PetscCall(MatMult(Na->A, x, y)); in MatMult_LRC_kernel()
32 PetscCall(MatMultAdd(U, Na->work1, y, y)); in MatMult_LRC_kernel()
34 PetscCall(MatMult(U, Na->work1, y)); in MatMult_LRC_kernel()
[all …]
|
| /petsc/src/ksp/ksp/utils/schurm/ |
| H A D | schurm.c | 7 Mat_SchurComplement *Na = (Mat_SchurComplement *)N->data; in MatCreateVecs_SchurComplement() local
10 if (Na->D) { in MatCreateVecs_SchurComplement()
11 PetscCall(MatCreateVecs(Na->D, right, left)); in MatCreateVecs_SchurComplement()
14 if (right) PetscCall(MatCreateVecs(Na->B, right, NULL)); in MatCreateVecs_SchurComplement()
15 if (left) PetscCall(MatCreateVecs(Na->C, NULL, left)); in MatCreateVecs_SchurComplement()
21 Mat_SchurComplement *Na = (Mat_SchurComplement *)N->data; in MatView_SchurComplement() local
25 if (Na->D) { in MatView_SchurComplement()
28 PetscCall(MatView(Na->D, viewer)); in MatView_SchurComplement()
35 PetscCall(MatView(Na->C, viewer)); in MatView_SchurComplement()
37 …e with the additional option -%sksp_view\n", ((PetscObject)Na->ksp)->prefix ? ((PetscObject)Na->ks… in MatView_SchurComplement()
[all …]
|
| /petsc/src/mat/impls/normal/ |
| H A D | normmh.c | 90 Mat_NormalHermitian *Na; in MatMult_NormalHermitian() local
93 PetscCall(MatShellGetContext(N, &Na)); in MatMult_NormalHermitian()
94 PetscCall(MatMult(Na->A, x, Na->w)); in MatMult_NormalHermitian()
95 PetscCall(MatMultHermitianTranspose(Na->A, Na->w, y)); in MatMult_NormalHermitian()
101 Mat_NormalHermitian *Na; in MatDestroy_NormalHermitian() local
104 PetscCall(MatShellGetContext(N, &Na)); in MatDestroy_NormalHermitian()
105 PetscCall(MatDestroy(&Na->A)); in MatDestroy_NormalHermitian()
106 PetscCall(MatDestroy(&Na->D)); in MatDestroy_NormalHermitian()
107 PetscCall(VecDestroy(&Na->w)); in MatDestroy_NormalHermitian()
108 PetscCall(PetscFree(Na)); in MatDestroy_NormalHermitian()
[all …]
|
| H A D | normm.c | 107 Mat_Normal *Na; in MatMult_Normal() local
110 PetscCall(MatShellGetContext(N, &Na)); in MatMult_Normal()
111 PetscCall(MatMult(Na->A, x, Na->w)); in MatMult_Normal()
112 PetscCall(MatMultTranspose(Na->A, Na->w, y)); in MatMult_Normal()
118 Mat_Normal *Na; in MatDestroy_Normal() local
121 PetscCall(MatShellGetContext(N, &Na)); in MatDestroy_Normal()
122 PetscCall(MatDestroy(&Na->A)); in MatDestroy_Normal()
123 PetscCall(MatDestroy(&Na->D)); in MatDestroy_Normal()
124 PetscCall(VecDestroy(&Na->w)); in MatDestroy_Normal()
125 PetscCall(PetscFree(Na)); in MatDestroy_Normal()
[all …]
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex82.c | 12 PetscInt i, N = 10, Na = N; in main() local
42 …PetscCall(KSPComputeRitz(ksp, harmonic ? PETSC_FALSE : PETSC_TRUE, PETSC_TRUE, &Na, S, tetar, teta… in main()
44 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "%% Number of Ritz pairs %" PetscInt_FMT "\n", Na)); in main()
45 for (i = 0; i < Na; i++) { in main()
|
| /petsc/src/sys/objects/ |
| H A D | options.c | 87 int Na; /* number or aliases */ member
1152 for (i = 0; i < options->Na; i++) { in PetscOptionsClear()
1156 options->Na = 0; in PetscOptionsClear()
1215 if (options->Na == options->Naalloc) { in PetscOptionsSetAlias()
1221 for (int i = 0; i < options->Na; ++i) { in PetscOptionsSetAlias()
1233 options->aliases1[options->Na] = (char *)malloc((len + 1) * sizeof(char)); in PetscOptionsSetAlias()
1234 PetscCall(PetscStrncpy(options->aliases1[options->Na], newname, len + 1)); in PetscOptionsSetAlias()
1236 options->aliases2[options->Na] = (char *)malloc((len + 1) * sizeof(char)); in PetscOptionsSetAlias()
1237 PetscCall(PetscStrncpy(options->aliases2[options->Na], oldname, len + 1)); in PetscOptionsSetAlias()
1238 ++options->Na; in PetscOptionsSetAlias()
[all …]
|
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 3379 PetscInt dim, qdim, d, Na, o, Nq, q, qp; in PetscQuadratureComputePermutations() local
3385 Na = DMPolytopeTypeGetNumArrangements(ct); in PetscQuadratureComputePermutations()
3386 PetscCall(PetscMalloc1(Na, perm)); in PetscQuadratureComputePermutations()
3387 if (Np) *Np = Na; in PetscQuadratureComputePermutations()
3388 Na /= 2; in PetscQuadratureComputePermutations()
3389 for (o = -Na; o < Na; ++o) { in PetscQuadratureComputePermutations()
3411 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, Nq, idx, PETSC_OWN_POINTER, &(*perm)[o + Na])); in PetscQuadratureComputePermutations()
3412 PetscCall(ISGetInfo((*perm)[o + Na], IS_PERMUTATION, IS_LOCAL, PETSC_TRUE, &flg)); in PetscQuadratureComputePermutations()
3414 PetscCall(ISSetPermutation((*perm)[o + Na])); in PetscQuadratureComputePermutations()
3416 if (!Na) (*perm)[0] = NULL; in PetscQuadratureComputePermutations()
|
| H A D | dtds.c | 655 const PetscInt Na = DMPolytopeTypeGetNumArrangements((DMPolytopeType)c); in PetscDSDestroy() local
657 for (PetscInt o = 0; o < Na; ++o) PetscCall(ISDestroy(&(*ds)->quadPerm[c][o])); in PetscDSDestroy()
3989 PetscInt Na, Nq; in PetscDSPermuteQuadPoint() local
3996 Na = DMPolytopeTypeGetNumArrangements(ct) / 2; in PetscDSPermuteQuadPoint()
3997 …Na && ornt < Na, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Orientation %" PetscInt_FMT " of %s i… in PetscDSPermuteQuadPoint()
3999 permIS = ds->quadPerm[(PetscInt)ct][ornt + Na]; in PetscDSPermuteQuadPoint()
|
| /petsc/src/dm/impls/plex/tests/ |
| H A D | ex5.c | 1248 PetscInt Na, a; in TestAssembly() local
1254 PetscCall(PetscQuadratureComputePermutations(quad, &Na, &perm)); in TestAssembly()
1255 for (a = 0; a < Na; ++a) PetscCall(ISDestroy(&perm[a])); in TestAssembly()
|