11   PetscInt    ncols;  in globalKMat_3d()  local
17         ncols = 0;  in globalKMat_3d()
25         vals[ncols++] = -6.; //ncols=1  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()
30         vals[ncols++] = 1.; //ncols=2  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()
35         vals[ncols++] = 1.; //ncols=3  in globalKMat_3d()
37         col[ncols].i  = i;  in globalKMat_3d()
38         col[ncols].j  = j - 1;  in globalKMat_3d()
39         col[ncols].k  = k;  in globalKMat_3d()
40         vals[ncols++] = 1.; //ncols=4  in globalKMat_3d()
42         col[ncols].i  = i;  in globalKMat_3d()
43         col[ncols].j  = j + 1;  in globalKMat_3d()
44         col[ncols].k  = k;  in globalKMat_3d()
45         vals[ncols++] = 1.; //ncols=5  in globalKMat_3d()
47         col[ncols].i  = i;  in globalKMat_3d()
48         col[ncols].j  = j;  in globalKMat_3d()
49         col[ncols].k  = k + 1;  in globalKMat_3d()
50         vals[ncols++] = 1.; //ncols=6  in globalKMat_3d()
52         col[ncols].i  = i;  in globalKMat_3d()
53         col[ncols].j  = j;  in globalKMat_3d()
54         col[ncols].k  = k - 1;  in globalKMat_3d()
55         vals[ncols++] = 1.; //ncols=7  in globalKMat_3d()
57         PetscCall(MatSetValuesStencil(K, 1, &row, ncols, col, vals, ADD_VALUES));  in globalKMat_3d()
70   PetscInt    ncols;  in globalKMat_2d()  local
75       ncols = 0;  in globalKMat_2d()
81       vals[ncols++] = -4.; //ncols=1  in globalKMat_2d()
83       col[ncols].i  = i - 1;  in globalKMat_2d()
84       col[ncols].j  = j;  in globalKMat_2d()
85       vals[ncols++] = 1.; //ncols=2  in globalKMat_2d()
87       col[ncols].i  = i;  in globalKMat_2d()
88       col[ncols].j  = j - 1;  in globalKMat_2d()
89       vals[ncols++] = 1.; //ncols=3  in globalKMat_2d()
91       col[ncols].i  = i + 1;  in globalKMat_2d()
92       col[ncols].j  = j;  in globalKMat_2d()
93       vals[ncols++] = 1.; //ncols=4  in globalKMat_2d()
95       col[ncols].i  = i;  in globalKMat_2d()
96       col[ncols].j  = j + 1;  in globalKMat_2d()
97       vals[ncols++] = 1.; //ncols=5  in globalKMat_2d()
99       PetscCall(MatSetValuesStencil(K, 1, &row, ncols, col, vals, ADD_VALUES));  in globalKMat_2d()