| /petsc/src/dm/impls/stag/tests/ |
| H A D | ex30.c | 147 DMStagStencil row, col[9]; in CreateMat() local
159 col[0].i = ex; in CreateMat()
160 col[0].j = ey; in CreateMat()
161 col[0].k = ez; in CreateMat()
162 col[0].loc = LEFT; in CreateMat()
163 col[0].c = 0; in CreateMat()
166 col[1].i = ex; in CreateMat()
167 col[1].j = ey + 1; in CreateMat()
168 col[1].k = ez; in CreateMat()
169 col[1].loc = LEFT; in CreateMat()
[all …]
|
| H A D | ex15.c | 150 DMStagStencil row, col[5]; in CreateSystem1d() local
157 col[0].i = e; in CreateSystem1d()
158 col[0].loc = DMSTAG_ELEMENT; in CreateSystem1d()
159 col[0].c = 0; in CreateSystem1d()
161 col[1].i = e - 1; in CreateSystem1d()
162 col[1].loc = DMSTAG_ELEMENT; in CreateSystem1d()
163 col[1].c = 0; in CreateSystem1d()
168 col[2].i = e; in CreateSystem1d()
169 col[2].loc = DMSTAG_LEFT; in CreateSystem1d()
170 col[2].c = 0; in CreateSystem1d()
[all …]
|
| H A D | ex18.c | 341 DMStagStencil row, col[7]; in CreateSystem() local
351 col[0].i = ex; in CreateSystem()
352 col[0].j = ey; in CreateSystem()
353 col[0].loc = DOWN; in CreateSystem()
354 col[0].c = 0; in CreateSystem()
356 col[1].i = ex; in CreateSystem()
357 col[1].j = ey - 1; in CreateSystem()
358 col[1].loc = DOWN; in CreateSystem()
359 col[1].c = 0; in CreateSystem()
361 col[2].i = ex; in CreateSystem()
[all …]
|
| /petsc/src/dm/impls/stag/tutorials/ |
| H A D | ex3.c | 265 DMStagStencil row, col[9]; in CreateSystem() local
277 col[0].i = ex; in CreateSystem()
278 col[0].j = ey; in CreateSystem()
279 col[0].k = ez; in CreateSystem()
280 col[0].loc = LEFT; in CreateSystem()
281 col[0].c = 0; in CreateSystem()
284 col[1].i = ex; in CreateSystem()
285 col[1].j = ey + 1; in CreateSystem()
286 col[1].k = ez; in CreateSystem()
287 col[1].loc = LEFT; in CreateSystem()
[all …]
|
| H A D | ex2.c | 223 DMStagStencil row, col[7]; in CreateSystem() local
233 col[0].i = ex; in CreateSystem()
234 col[0].j = ey; in CreateSystem()
235 col[0].loc = DOWN; in CreateSystem()
236 col[0].c = 0; in CreateSystem()
238 col[1].i = ex; in CreateSystem()
239 col[1].j = ey - 1; in CreateSystem()
240 col[1].loc = DOWN; in CreateSystem()
241 col[1].c = 0; in CreateSystem()
243 col[2].i = ex; in CreateSystem()
[all …]
|
| H A D | ex4.c | 652 DMStagStencil row, col[11]; in CreateSystem2d() local
701 col[count] = row; in CreateSystem2d()
707 col[count].i = ex; in CreateSystem2d()
708 col[count].j = ey - 1; in CreateSystem2d()
709 col[count].loc = DMSTAG_DOWN; in CreateSystem2d()
710 col[count].c = 0; in CreateSystem2d()
713 col[count].i = ex; in CreateSystem2d()
714 col[count].j = ey + 1; in CreateSystem2d()
715 col[count].loc = DMSTAG_DOWN; in CreateSystem2d()
716 col[count].c = 0; in CreateSystem2d()
[all …]
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex17.c | 122 PetscInt i, j, Ii, J, col[5], Istart, Iend; in FormTestMatrix() local
131 col[0] = i - 1; in FormTestMatrix()
132 col[1] = i; in FormTestMatrix()
133 col[2] = i + 1; in FormTestMatrix()
134 PetscCall(MatSetValues(A, 1, &i, 3, col, val, INSERT_VALUES)); in FormTestMatrix()
137 col[0] = n - 2; in FormTestMatrix()
138 col[1] = n - 1; in FormTestMatrix()
139 PetscCall(MatSetValues(A, 1, &i, 2, col, val, INSERT_VALUES)); in FormTestMatrix()
141 col[0] = 0; in FormTestMatrix()
142 col[1] = 1; in FormTestMatrix()
[all …]
|
| H A D | ex32.c | 129 MatStencil row, col; in ComputeMatrix() local
171 col.i = i; in ComputeMatrix()
172 col.j = j; in ComputeMatrix()
173 col.k = k; in ComputeMatrix()
174 PetscCall(MatSetValuesBlockedStencil(B, 1, &row, 1, &col, v, INSERT_VALUES)); in ComputeMatrix()
177 col.i = i - 1; in ComputeMatrix()
178 col.j = j; in ComputeMatrix()
179 col.k = k; in ComputeMatrix()
180 PetscCall(MatSetValuesBlockedStencil(B, 1, &row, 1, &col, v_neighbor, INSERT_VALUES)); in ComputeMatrix()
181 col.i = i + 1; in ComputeMatrix()
[all …]
|
| H A D | ex55.c | 8 PetscInt n = 10, i, col[3]; in main() local
30 col[0] = i - 1; in main()
31 col[1] = i; in main()
32 col[2] = i + 1; in main()
33 PetscCall(MatSetValues(A, 1, &i, 3, col, value, INSERT_VALUES)); in main()
36 col[0] = n - 2; in main()
37 col[1] = n - 1; in main()
38 PetscCall(MatSetValues(A, 1, &i, 2, col, value, INSERT_VALUES)); in main()
40 col[0] = 0; in main()
41 col[1] = 1; in main()
[all …]
|
| /petsc/src/dm/tests/ |
| H A D | noflux_check.c | 9 MatStencil row, col[7]; in globalKMat_3d() local
22 col[0].i = i; in globalKMat_3d()
23 col[0].j = j; in globalKMat_3d()
24 col[0].k = k; in globalKMat_3d()
27 col[ncols].i = i - 1; in globalKMat_3d()
28 col[ncols].j = j; in globalKMat_3d()
29 col[ncols].k = k; in globalKMat_3d()
32 col[ncols].i = i + 1; in globalKMat_3d()
33 col[ncols].j = j; in globalKMat_3d()
34 col[ncols].k = k; in globalKMat_3d()
[all …]
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex34.c | 137 MatStencil row, col[7]; in ComputeMatrix() local
167 col[num].i = i; in ComputeMatrix()
168 col[num].j = j; in ComputeMatrix()
169 col[num].k = k - 1; in ComputeMatrix()
170 col[num].c = d; in ComputeMatrix()
176 col[num].i = i; in ComputeMatrix()
177 col[num].j = j - 1; in ComputeMatrix()
178 col[num].k = k; in ComputeMatrix()
179 col[num].c = d; in ComputeMatrix()
185 col[num].i = i - 1; in ComputeMatrix()
[all …]
|
| H A D | ex22f.F90 | 47 MatStencil row(1), col(7)
73 col(1)%i = i
74 col(1)%j = j
75 col(1)%k = k - 1
77 col(2)%i = i
78 col(2)%j = j - 1
79 col(2)%k = k
81 col(3)%i = i - 1
82 col(3)%j = j
83 col(3)%k = k
[all …]
|
| H A D | ex28.c | 47 PetscInt mx, col[2], xs, xm, i; in ComputeInitialSolution() local
56 col[0] = 2 * i; in ComputeInitialSolution()
57 col[1] = 2 * i + 1; in ComputeInitialSolution()
59 PetscCall(VecSetValues(x, 2, col, val, INSERT_VALUES)); in ComputeInitialSolution()
87 MatStencil row, col[7]; in ComputeMatrix() local
93 PetscCall(PetscArrayzero(col, 7)); in ComputeMatrix()
104 col[0].i = i; in ComputeMatrix()
105 col[0].c = 0; in ComputeMatrix()
107 col[1].i = i - 1; in ComputeMatrix()
108 col[1].c = 1; in ComputeMatrix()
[all …]
|
| H A D | ex50.c | 106 MatStencil row, col[5]; in ComputeJacobian() local
129 col[num].i = i; in ComputeJacobian()
130 col[num].j = j - 1; in ComputeJacobian()
136 col[num].i = i - 1; in ComputeJacobian()
137 col[num].j = j; in ComputeJacobian()
143 col[num].i = i + 1; in ComputeJacobian()
144 col[num].j = j; in ComputeJacobian()
150 col[num].i = i; in ComputeJacobian()
151 col[num].j = j + 1; in ComputeJacobian()
156 col[num].i = i; in ComputeJacobian()
[all …]
|
| H A D | ex32.c | 122 MatStencil row, col[5]; in ComputeMatrix() local
148 col[num].i = i; in ComputeMatrix()
149 col[num].j = j - 1; in ComputeMatrix()
155 col[num].i = i - 1; in ComputeMatrix()
156 col[num].j = j; in ComputeMatrix()
162 col[num].i = i + 1; in ComputeMatrix()
163 col[num].j = j; in ComputeMatrix()
169 col[num].i = i; in ComputeMatrix()
170 col[num].j = j + 1; in ComputeMatrix()
175 col[num].i = i; in ComputeMatrix()
[all …]
|
| H A D | ex45.c | 122 MatStencil row, col[7]; in ComputeMatrix() local
147 col[0].i = i; in ComputeMatrix()
148 col[0].j = j; in ComputeMatrix()
149 col[0].k = k - 1; in ComputeMatrix()
151 col[1].i = i; in ComputeMatrix()
152 col[1].j = j - 1; in ComputeMatrix()
153 col[1].k = k; in ComputeMatrix()
155 col[2].i = i - 1; in ComputeMatrix()
156 col[2].j = j; in ComputeMatrix()
157 col[2].k = k; in ComputeMatrix()
[all …]
|
| H A D | ex53.c | 15 PetscInt i, n = 10, col[3], its; in main() local
46 col[0] = i - 1; in main()
47 col[1] = i; in main()
48 col[2] = i + 1; in main()
49 PetscCall(MatSetValues(A, 1, &i, 3, col, value, INSERT_VALUES)); in main()
52 col[0] = n - 2; in main()
53 col[1] = n - 1; in main()
54 PetscCall(MatSetValues(A, 1, &i, 2, col, value, INSERT_VALUES)); in main()
56 col[0] = 0; in main()
57 col[1] = 1; in main()
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | ij.c | 33 PetscInt *work, *ia, *ja, *j, i, nz, row, col; in MatToSymmetricIJ_SeqAIJ() local
47 col = *j++ + shiftin; in MatToSymmetricIJ_SeqAIJ()
49 if (col > row) break; in MatToSymmetricIJ_SeqAIJ()
51 if (col < row) break; in MatToSymmetricIJ_SeqAIJ()
53 if (col != row) ia[row + 1]++; in MatToSymmetricIJ_SeqAIJ()
54 ia[col + 1]++; in MatToSymmetricIJ_SeqAIJ()
75 col = *j++ + shiftin; in MatToSymmetricIJ_SeqAIJ()
77 if (col > row) break; in MatToSymmetricIJ_SeqAIJ()
79 if (col < row) break; in MatToSymmetricIJ_SeqAIJ()
81 if (col != row) ja[work[col]++] = row + shiftout; in MatToSymmetricIJ_SeqAIJ()
[all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex18.c | 301 MatStencil col[5], row; in FormJacobian() local
361 col[0].j = j - 1; in FormJacobian()
362 col[0].i = i; in FormJacobian()
364 col[1].j = j; in FormJacobian()
365 col[1].i = i - 1; in FormJacobian()
367 col[2].j = row.j = j; in FormJacobian()
368 col[2].i = row.i = i; in FormJacobian()
370 col[3].j = j; in FormJacobian()
371 col[3].i = i + 1; in FormJacobian()
373 col[4].j = j + 1; in FormJacobian()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex77.c | 10 PetscInt n, mbs = 16, bs = 1, nz = 3, prob = 2, i, j, col[3], row, Ii, J, n1; in main() local
46 col[0] = i - 1; in main()
47 col[1] = i; in main()
48 col[2] = i + 1; in main()
49 PetscCall(MatSetValues(A, 1, &i, 3, col, value, INSERT_VALUES)); in main()
50 PetscCall(MatSetValues(sA, 1, &i, 3, col, value, INSERT_VALUES)); in main()
53 col[0] = 0; in main()
54 col[1] = n - 2; in main()
55 col[2] = n - 1; in main()
60 PetscCall(MatSetValues(A, 1, &i, 3, col, value, INSERT_VALUES)); in main()
[all …]
|
| H A D | ex161.c | 9 PetscInt row, col, m, n; in main() local
29 col = 0; in main()
31 PetscCall(MatSetValues(A, 1, &row, 1, &col, &val, ADD_VALUES)); in main()
33 col = 3; in main()
35 PetscCall(MatSetValues(A, 1, &row, 1, &col, &val, ADD_VALUES)); in main()
37 col = 2; in main()
39 PetscCall(MatSetValues(A, 1, &row, 1, &col, &val, ADD_VALUES)); in main()
41 col = 0; in main()
43 PetscCall(MatSetValues(A, 1, &row, 1, &col, &val, ADD_VALUES)); in main()
57 col = 0; in main()
[all …]
|
| H A D | ex76.c | 10 PetscInt n, mbs = 16, bs = 1, nz = 3, prob = 1, i, j, col[3], block, row, Ii, J, n1, lvl; in main() local
50 col[0] = i - 1; in main()
51 col[1] = i; in main()
52 col[2] = i + 1; in main()
53 PetscCall(MatSetValues(A, 1, &i, 3, col, value, INSERT_VALUES)); in main()
56 col[0] = 0; in main()
57 col[1] = n - 2; in main()
58 col[2] = n - 1; in main()
63 PetscCall(MatSetValues(A, 1, &i, 3, col, value, INSERT_VALUES)); in main()
66 col[0] = 0; in main()
[all …]
|
| H A D | ex52.c | 8 PetscInt m = 2, bs = 1, M, row, col, start, end, i, j, k; in main() local
31 …for (col = start; col < end; col++, data += 1) PetscCall(MatSetValues(A, 1, &row, 1, &col, &data, … in main()
39 for (col = 0; col < M; col++) PetscCall(MatSetValues(A, 1, &row, 1, &col, &data, ADD_VALUES)); in main()
48 col = start / bs; in main()
58 PetscCall(MatSetValuesBlocked(A, 1, &row, 1, &col, bval, INSERT_VALUES)); in main()
|
| /petsc/src/ksp/pc/tests/ |
| H A D | ex4.c | 11 PetscInt n = 5, i, col[3]; in main() local
26 col[0] = i - 1; in main()
27 col[1] = i; in main()
28 col[2] = i + 1; in main()
29 PetscCall(MatSetValues(mat, 1, &i, 3, col, value, INSERT_VALUES)); in main()
32 col[0] = n - 2; in main()
33 col[1] = n - 1; in main()
34 PetscCall(MatSetValues(mat, 1, &i, 2, col, value, INSERT_VALUES)); in main()
36 col[0] = 0; in main()
37 col[1] = 1; in main()
[all …]
|
| /petsc/src/ts/tutorials/power_grid/ |
| H A D | ex6.c | 313 MatStencil row, col[5]; in IJacobian() local
332 col[nc].i = i; in IJacobian()
333 col[nc].j = j; in IJacobian()
340 col[nc].i = i + 1; in IJacobian()
341 col[nc].j = j; in IJacobian()
343 col[nc].i = i; in IJacobian()
344 col[nc].j = j + 1; in IJacobian()
346 col[nc].i = i; in IJacobian()
347 col[nc].j = j; in IJacobian()
350 col[nc].i = i + 1; in IJacobian()
[all …]
|