| /petsc/src/mat/tests/ |
| H A D | ex230.c | 42 PetscInt ii, jj; in ex1_nonsquare_bs1() local
45 ii = 3; in ex1_nonsquare_bs1()
47 PetscCall(MatSetValues(preallocator, 1, &ii, 1, &jj, &vv, INSERT_VALUES)); in ex1_nonsquare_bs1()
49 ii = 7; in ex1_nonsquare_bs1()
51 PetscCall(MatSetValues(preallocator, 1, &ii, 1, &jj, &vv, INSERT_VALUES)); in ex1_nonsquare_bs1()
53 ii = 9; in ex1_nonsquare_bs1()
55 PetscCall(MatSetValue(preallocator, ii, jj, vv, INSERT_VALUES)); in ex1_nonsquare_bs1()
78 PetscInt ii, jj; in ex1_nonsquare_bs1() local
81 ii = 3; in ex1_nonsquare_bs1()
84 PetscCall(MatSetValue(A, ii, jj, vv, INSERT_VALUES)); in ex1_nonsquare_bs1()
[all …]
|
| H A D | ex259.c | 37 PetscInt ii, jj; in ex1_nonsquare_bs1() local
40 ii = 3; in ex1_nonsquare_bs1()
42 PetscCall(MatSetValue(A, ii, jj, vv, INSERT_VALUES)); in ex1_nonsquare_bs1()
44 ii = 7; in ex1_nonsquare_bs1()
46 PetscCall(MatSetValue(A, ii, jj, vv, INSERT_VALUES)); in ex1_nonsquare_bs1()
47 PetscCall(MatSetValue(A, jj, ii, vv, INSERT_VALUES)); in ex1_nonsquare_bs1()
49 ii = 9; in ex1_nonsquare_bs1()
51 PetscCall(MatSetValue(A, ii, jj, vv, INSERT_VALUES)); in ex1_nonsquare_bs1()
52 PetscCall(MatSetValue(A, jj, ii, vv, INSERT_VALUES)); in ex1_nonsquare_bs1()
64 PetscInt ii, jj; in ex1_nonsquare_bs1() local
[all …]
|
| H A D | ex137.c | 8 PetscInt m = 4, bs = 1, ii[5], jj[7]; in main() local
19 ii[0] = 0; in main()
20 ii[1] = 2; in main()
21 ii[2] = 5; in main()
22 ii[3] = 7; in main()
23 ii[4] = 7; in main()
42 ii[0] = 0; in main()
43 ii[1] = 2; in main()
44 ii[2] = 4; in main()
45 ii[3] = 6; in main()
[all …]
|
| H A D | ex260.c | 10 PetscInt *ii, *jd, *jo; in main() local
20 PetscCall(PetscMalloc4(n + 1, &ii, n * nzr, &jd, n * nzr, &jo, n * nzr, &a)); in main()
21 ii[0] = 0; in main()
23 ii[i + 1] = ii[i] + nzr; in main()
24 for (PetscInt j = 0; j < nzr; j++) jd[ii[i] + j] = j; in main()
26 for (PetscInt j = 0; j < nzr; j++) jo[ii[i] + j] = n + j - 1; in main()
28 for (PetscInt j = 0; j < nzr; j++) jo[ii[i] + j] = j; in main()
31 …thSplitArrays(PETSC_COMM_WORLD, n, n, PETSC_DETERMINE, PETSC_DETERMINE, ii, jd, a, ii, jo, a, &A)); in main()
34 PetscCall(PetscFree4(ii, jd, jo, a)); in main()
|
| H A D | ex195.c | 17 PetscInt m, M, n, istart, iend, ii, i, J, j, K = 10; in main() local
36 for (ii = istart; ii < iend; ii++) { in main()
38 i = ii / n; in main()
39 j = ii - i * n; in main()
41 J = ii - n; in main()
42 PetscCall(MatSetValues(A1, 1, &ii, 1, &J, &v, INSERT_VALUES)); in main()
45 J = ii + n; in main()
46 PetscCall(MatSetValues(A1, 1, &ii, 1, &J, &v, INSERT_VALUES)); in main()
49 J = ii - 1; in main()
50 PetscCall(MatSetValues(A1, 1, &ii, 1, &J, &v, INSERT_VALUES)); in main()
[all …]
|
| /petsc/src/sys/tests/ |
| H A D | ex63.c | 18 PetscInt64 ii, jj, count = 0; in Prime() local
25 for (ii = 2; ii < n + 1; ii++) { in Prime()
27 for (jj = 2; jj <= PetscFloorReal(PetscSqrtReal(ii)); jj++) { in Prime()
28 if (ii % jj == 0) { in Prime()
34 prime[count] = ii; in Prime()
41 for (ii = 1; ii < count + 1; ii++) (*set)[ii] = prime[ii - 1]; in Prime()
50 PetscInt ii; in PrintSet() local
54 for (ii = 1; ii <= (PetscInt)set[0]; ii++) { in PrintSet()
56 PetscCall(PetscSynchronizedPrintf(comm, text, set[ii])); in PrintSet()
66 PetscInt ii; in AssertSetsEqual() local
[all …]
|
| H A D | ex62.c | 25 PetscInt ii; in Fibonnaci() local
34 for (ii = 0; ii < n; ii++) (*set)[ii + 1] = fib[ii]; in Fibonnaci()
41 PetscInt64 ii; in Square() local
46 for (ii = 1; ii < n + 1; ii++) (*set)[ii] = ii * ii; in Square()
53 PetscInt64 ii; in Cube() local
58 for (ii = 1; ii < n + 1; ii++) (*set)[ii] = ii * ii * ii; in Cube()
65 PetscInt64 ii; in Sixth() local
70 for (ii = 1; ii < n + 1; ii++) (*set)[ii] = ii * ii * ii * ii * ii * ii; in Sixth()
81 for (PetscInt64 ii = 1; ii <= set[0]; ii++) { in PrintSet() local
83 PetscCall(PetscPrintf(PETSC_COMM_WORLD, text, set[ii])); in PrintSet()
[all …]
|
| /petsc/src/sys/objects/ |
| H A D | garbage.c | 95 PetscInt ii, jj = 0, counter = 0; in GarbageKeySortedIntersect_Private() local
106 for (ii = 0; ii < *lena; ii++) { in GarbageKeySortedIntersect_Private()
107 while (jj < lenb && seta[ii] > setb[jj]) jj++; in GarbageKeySortedIntersect_Private()
109 if (seta[ii] == setb[jj]) { in GarbageKeySortedIntersect_Private()
110 seta[counter] = seta[ii]; in GarbageKeySortedIntersect_Private()
137 PetscInt ii, max_entries; in GarbageKeyAllReduceIntersect_Private() local
151 for (ii = 1; ii < *entries + 1; ii++) sendset[ii] = set[ii - 1]; in GarbageKeyAllReduceIntersect_Private()
165 for (ii = 0; ii < *entries; ii++) set[ii] = recvset[ii + 1]; in GarbageKeyAllReduceIntersect_Private()
208 PetscInt ii, entries, offset; in PetscGarbageCleanup() local
233 for (ii = 0; ii < entries; ii++) { in PetscGarbageCleanup()
[all …]
|
| /petsc/src/dm/impls/plex/tutorials/ |
| H A D | ex18.c | 28 for (PetscInt ii = cStart; ii < cEnd; ++ii) { in surfArea() local
29 PetscCall(DMLabelGetValue(faceLabel, ii, &faceID)); in surfArea()
31 PetscCall(DMPlexComputeCellGeometryFVM(dm, ii, &vol, centroid, normal)); in surfArea()
39 for (PetscInt ii = fStart; ii < fEnd; ++ii) { in surfArea() local
40 PetscCall(DMLabelGetValue(faceLabel, ii, &faceID)); in surfArea()
42 PetscCall(DMPlexComputeCellGeometryFVM(dm, ii, &vol, centroid, normal)); in surfArea()
48 for (PetscInt ii = cStart; ii < cEnd; ++ii) { in surfArea() local
49 PetscCall(DMLabelGetValue(bodyLabel, ii, &bodyID)); in surfArea()
51 PetscCall(DMPlexComputeCellGeometryFVM(dm, ii, &vol, centroid, normal)); in surfArea()
151 for (PetscInt ii = 0; ii < numCntrlPnts; ++ii) { in main() local
[all …]
|
| H A D | ex19.c | 31 for (PetscInt ii = cStart; ii < cEnd; ++ii) { in surfArea() local
32 PetscCall(DMLabelGetValue(faceLabel, ii, &faceID)); in surfArea()
34 PetscCall(DMPlexComputeCellGeometryFVM(dm, ii, &vol, centroid, normal)); in surfArea()
42 for (PetscInt ii = fStart; ii < fEnd; ++ii) { in surfArea() local
43 PetscCall(DMLabelGetValue(faceLabel, ii, &faceID)); in surfArea()
45 PetscCall(DMPlexComputeCellGeometryFVM(dm, ii, &vol, centroid, normal)); in surfArea()
51 for (PetscInt ii = cStart; ii < cEnd; ++ii) { in surfArea() local
52 PetscCall(DMLabelGetValue(bodyLabel, ii, &bodyID)); in surfArea()
54 PetscCall(DMPlexComputeCellGeometryFVM(dm, ii, &vol, centroid, normal)); in surfArea()
215 for (PetscInt ii = 0; ii < numRows; ++ii) { in main() local
[all …]
|
| /petsc/src/dm/tests/ |
| H A D | ex28.c | 34 PetscInt func, ii; in main() local
97 for (ii = 0; ii < 3; ++ii) { in main()
98 PetscCall(DMGetGlobalVector(da1, &xxsplit[ii])); in main()
99 PetscCall(PetscObjectSetName((PetscObject)xxsplit[ii], "Real space split vector")); in main()
100 PetscCall(DMGetGlobalVector(da1, &yysplit[ii])); in main()
101 PetscCall(PetscObjectSetName((PetscObject)yysplit[ii], "FFTW frequency space split vector")); in main()
102 PetscCall(DMGetGlobalVector(da1, &zzsplit[ii])); in main()
103 PetscCall(PetscObjectSetName((PetscObject)zzsplit[ii], "FFTW reconstructed split vector")); in main()
147 for (ii = 0; ii < 3; ++ii) PetscCall(MatMult(AA, xxsplit[ii], zzsplit[ii])); in main()
151 for (ii = 0; ii < 3; ++ii) PetscCall(MatMult(AA, xxsplit[ii], yysplit[ii])); in main()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijsolvtran5.c | 10 PetscInt i, nz, idx, idt, ii, ic, ir, oidx; in MatSolveTranspose_SeqBAIJ_5_inplace() local
26 ii = 0; in MatSolveTranspose_SeqBAIJ_5_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_5_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_5_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_5_inplace()
32 t[ii + 3] = b[ic + 3]; in MatSolveTranspose_SeqBAIJ_5_inplace()
33 t[ii + 4] = b[ic + 4]; in MatSolveTranspose_SeqBAIJ_5_inplace()
34 ii += 5; in MatSolveTranspose_SeqBAIJ_5_inplace()
95 ii = 0; in MatSolveTranspose_SeqBAIJ_5_inplace()
98 x[ir] = t[ii]; in MatSolveTranspose_SeqBAIJ_5_inplace()
[all …]
|
| H A D | baijsolvtran7.c | 10 PetscInt i, nz, idx, idt, ii, ic, ir, oidx; in MatSolveTranspose_SeqBAIJ_7_inplace() local
26 ii = 0; in MatSolveTranspose_SeqBAIJ_7_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_7_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_7_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_7_inplace()
32 t[ii + 3] = b[ic + 3]; in MatSolveTranspose_SeqBAIJ_7_inplace()
33 t[ii + 4] = b[ic + 4]; in MatSolveTranspose_SeqBAIJ_7_inplace()
34 t[ii + 5] = b[ic + 5]; in MatSolveTranspose_SeqBAIJ_7_inplace()
35 t[ii + 6] = b[ic + 6]; in MatSolveTranspose_SeqBAIJ_7_inplace()
36 ii += 7; in MatSolveTranspose_SeqBAIJ_7_inplace()
[all …]
|
| H A D | baijsolvtran4.c | 10 PetscInt i, nz, idx, idt, ii, ic, ir, oidx; in MatSolveTranspose_SeqBAIJ_4_inplace() local
26 ii = 0; in MatSolveTranspose_SeqBAIJ_4_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_4_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_4_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_4_inplace()
32 t[ii + 3] = b[ic + 3]; in MatSolveTranspose_SeqBAIJ_4_inplace()
33 ii += 4; in MatSolveTranspose_SeqBAIJ_4_inplace()
88 ii = 0; in MatSolveTranspose_SeqBAIJ_4_inplace()
91 x[ir] = t[ii]; in MatSolveTranspose_SeqBAIJ_4_inplace()
92 x[ir + 1] = t[ii + 1]; in MatSolveTranspose_SeqBAIJ_4_inplace()
[all …]
|
| H A D | baijsolvtran6.c | 10 PetscInt i, nz, idx, idt, ii, ic, ir, oidx; in MatSolveTranspose_SeqBAIJ_6_inplace() local
26 ii = 0; in MatSolveTranspose_SeqBAIJ_6_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_6_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_6_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_6_inplace()
32 t[ii + 3] = b[ic + 3]; in MatSolveTranspose_SeqBAIJ_6_inplace()
33 t[ii + 4] = b[ic + 4]; in MatSolveTranspose_SeqBAIJ_6_inplace()
34 t[ii + 5] = b[ic + 5]; in MatSolveTranspose_SeqBAIJ_6_inplace()
35 ii += 6; in MatSolveTranspose_SeqBAIJ_6_inplace()
102 ii = 0; in MatSolveTranspose_SeqBAIJ_6_inplace()
[all …]
|
| H A D | baijsolvtran3.c | 10 PetscInt i, nz, idx, idt, ii, ic, ir, oidx; in MatSolveTranspose_SeqBAIJ_3_inplace() local
26 ii = 0; in MatSolveTranspose_SeqBAIJ_3_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_3_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_3_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_3_inplace()
32 ii += 3; in MatSolveTranspose_SeqBAIJ_3_inplace()
81 ii = 0; in MatSolveTranspose_SeqBAIJ_3_inplace()
84 x[ir] = t[ii]; in MatSolveTranspose_SeqBAIJ_3_inplace()
85 x[ir + 1] = t[ii + 1]; in MatSolveTranspose_SeqBAIJ_3_inplace()
86 x[ir + 2] = t[ii + 2]; in MatSolveTranspose_SeqBAIJ_3_inplace()
[all …]
|
| H A D | baijsolvtran2.c | 10 PetscInt i, nz, idx, idt, ii, ic, ir, oidx; in MatSolveTranspose_SeqBAIJ_2_inplace() local
26 ii = 0; in MatSolveTranspose_SeqBAIJ_2_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_2_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_2_inplace()
31 ii += 2; in MatSolveTranspose_SeqBAIJ_2_inplace()
74 ii = 0; in MatSolveTranspose_SeqBAIJ_2_inplace()
77 x[ir] = t[ii]; in MatSolveTranspose_SeqBAIJ_2_inplace()
78 x[ir + 1] = t[ii + 1]; in MatSolveTranspose_SeqBAIJ_2_inplace()
79 ii += 2; in MatSolveTranspose_SeqBAIJ_2_inplace()
96 PetscInt nz, idx, idt, j, i, oidx, ii, ic, ir; in MatSolveTranspose_SeqBAIJ_2() local
[all …]
|
| H A D | baij2.c | 243 const PetscInt *idx, *ii, *ridx = NULL; in MatMult_SeqBAIJ_1() local
252 ii = a->compressedrow.i; in MatMult_SeqBAIJ_1()
257 ii = a->i; in MatMult_SeqBAIJ_1()
261 n = ii[1] - ii[0]; in MatMult_SeqBAIJ_1()
262 v = a->a + ii[0]; in MatMult_SeqBAIJ_1()
263 idx = a->j + ii[0]; in MatMult_SeqBAIJ_1()
264 ii++; in MatMult_SeqBAIJ_1()
288 PetscInt mbs, i, *idx, *ii, j, n, *ridx = NULL; in MatMult_SeqBAIJ_2() local
299 ii = a->compressedrow.i; in MatMult_SeqBAIJ_2()
304 ii = a->i; in MatMult_SeqBAIJ_2()
[all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex16.c | 284 PetscInt ii, jj, kk, l; in DeformationGradient() local
292 for (ii = 0; ii < NB; ii++) { in DeformationGradient()
293 PetscInt idx = ii + jj * NB + kk * NB * NB; in DeformationGradient()
311 void DeformationGradientJacobian(PetscInt qi, PetscInt qj, PetscInt qk, PetscInt ii, PetscInt jj, P… in DeformationGradientJacobian() argument
315 PetscInt idx = ii + jj * NB + kk * NB * NB; in DeformationGradientJacobian()
372 PetscInt i, j, k, ii, jj, kk; in FormElements() local
383 for (ii = 0; ii < NQ; ii++) { in FormElements()
384 … PetscInt idx = ii + NQ * jj + NQ * NQ * kk + NEQ * i + NEQ * NB * j + NEQ * NB * NB * k; in FormElements()
385 bx = pts[ii]; in FormElements()
419 PetscInt ii, jj, kk; in GatherElementData() local
[all …]
|
| /petsc/src/ksp/ksp/impls/gmres/lgmres/ |
| H A D | lgmres.c | 105 PetscInt ii, jj; in KSPLGMRESCycle() local
173 for (ii = 0; ii < aug_dim; ii++) { in KSPLGMRESCycle()
174 if (lgmres->aug_order[ii] == order) { in KSPLGMRESCycle()
175 spot = ii; in KSPLGMRESCycle()
261 for (ii = 0; ii < aug_dim; ii++) { in KSPLGMRESCycle()
262 if (lgmres->aug_order[ii] == aug_dim) spot = ii; in KSPLGMRESCycle()
275 for (ii = 0; ii < aug_dim; ii++) AUG_ORDER(ii)++; in KSPLGMRESCycle()
284 for (ii = 0; ii < it_total + 1; ii++) { in KSPLGMRESCycle()
285 for (jj = 0; jj <= ii + 1 && jj < it_total + 1; jj++) avec[jj] += *HES(jj, ii) * *GRS(ii); in KSPLGMRESCycle()
303 PetscInt ii; /* LGMRES_MOD variable */ in KSPSolve_LGMRES() local
[all …]
|
| /petsc/src/mat/impls/aij/seq/ftn-kernels/ |
| H A D | fmult.F90 | 12 PetscInt, intent(in) :: n, ii(*), jj(0:*)
16 jend = ii(1)
19 jend = ii(i + 1)
29 PetscInt, intent(in) :: n, ii(*), jj(0:*)
37 jstart = ii(i)
38 jend = ii(i + 1)
|
| /petsc/src/vec/is/is/tests/ |
| H A D | ex2.c | 13 const PetscInt *ii; in main() local
31 PetscCall(ISGetIndices(is, &ii)); in main()
32 PetscCall(ISRestoreIndices(is, &ii)); in main()
40 PetscCall(ISGetIndices(is, &ii)); in main()
41 …for (i = 0; i < 10000; i++) PetscCheck(ii[i] == -8 + 3 * i, PETSC_COMM_SELF, PETSC_ERR_PLIB, "ISGe… in main()
42 PetscCall(ISRestoreIndices(is, &ii)); in main()
|
| H A D | ex2f.F90 | 14 PetscInt, pointer :: ii(:)
32 PetscCallA(ISGetIndices(is, ii, ierr))
33 PetscCallA(ISRestoreIndices(is, ii, ierr))
43 PetscCallA(ISGetIndices(is, ii, ierr))
45 … PetscCheckA(ii(i) == -11 + 3*i, PETSC_COMM_SELF, PETSC_ERR_PLIB, 'Wrong result from ISGetIndices')
47 PetscCallA(ISRestoreIndices(is, ii, ierr))
|
| /petsc/src/mat/graphops/coarsen/impls/hem/ |
| H A D | hem.c | 204 PetscInt ii; in PetscCDPrint() local
207 for (ii = 0; ii < ail->size; ii++) { in PetscCDPrint()
208 n2 = n = ail->array[ii]; in PetscCDPrint()
209 if (n) PetscCall(PetscSynchronizedPrintf(comm, "list %" PetscInt_FMT ":", ii + Istart)); in PetscCDPrint()
288 for (PetscInt ii = 0; ii < ail->size; ii++) { in PetscCDCount() local
289 PetscCDIntNd *n1 = ail->array[ii]; in PetscCDCount()
315 PetscInt ii, kk; in PetscCDGetNonemptyIS() local
319 for (ii = kk = 0; ii < ail->size; ii++) { in PetscCDGetNonemptyIS()
320 n = ail->array[ii]; in PetscCDGetNonemptyIS()
324 for (ii = kk = 0; ii < ail->size; ii++) { in PetscCDGetNonemptyIS()
[all …]
|
| /petsc/src/mat/tutorials/ |
| H A D | ex11.c | 17 PetscInt *ii, *jj; in main() local
28 PetscCall(PetscMalloc1(3, &ii)); in main()
30 ii[0] = 0; in main()
31 ii[1] = 3; in main()
32 ii[2] = 6; in main()
48 PetscCall(MatCreateMPIAdj(MPI_COMM_WORLD, ncells, Nvertices, ii, jj, NULL, &mesh)); in main()
|