11   PetscScalar       *yy, x0, x1, x2, x3, x4, x5, x6;  in PCApply_PBJacobi()  local
16   PetscCall(VecGetArray(y, &yy));  in PCApply_PBJacobi()
19     for (i = 0; i < m; i++) yy[i] = diag[i] * xx[i];  in PCApply_PBJacobi()
25       yy[2 * i]     = diag[0] * x0 + diag[2] * x1;  in PCApply_PBJacobi()
26       yy[2 * i + 1] = diag[1] * x0 + diag[3] * x1;  in PCApply_PBJacobi()
36       yy[3 * i]     = diag[0] * x0 + diag[3] * x1 + diag[6] * x2;  in PCApply_PBJacobi()
37       yy[3 * i + 1] = diag[1] * x0 + diag[4] * x1 + diag[7] * x2;  in PCApply_PBJacobi()
38       yy[3 * i + 2] = diag[2] * x0 + diag[5] * x1 + diag[8] * x2;  in PCApply_PBJacobi()
49       yy[4 * i]     = diag[0] * x0 + diag[4] * x1 + diag[8] * x2 + diag[12] * x3;  in PCApply_PBJacobi()
50       yy[4 * i + 1] = diag[1] * x0 + diag[5] * x1 + diag[9] * x2 + diag[13] * x3;  in PCApply_PBJacobi()
51       yy[4 * i + 2] = diag[2] * x0 + diag[6] * x1 + diag[10] * x2 + diag[14] * x3;  in PCApply_PBJacobi()
52       yy[4 * i + 3] = diag[3] * x0 + diag[7] * x1 + diag[11] * x2 + diag[15] * x3;  in PCApply_PBJacobi()
64       yy[5 * i]     = diag[0] * x0 + diag[5] * x1 + diag[10] * x2 + diag[15] * x3 + diag[20] * x4;  in PCApply_PBJacobi()
65       yy[5 * i + 1] = diag[1] * x0 + diag[6] * x1 + diag[11] * x2 + diag[16] * x3 + diag[21] * x4;  in PCApply_PBJacobi()
66       yy[5 * i + 2] = diag[2] * x0 + diag[7] * x1 + diag[12] * x2 + diag[17] * x3 + diag[22] * x4;  in PCApply_PBJacobi()
67       yy[5 * i + 3] = diag[3] * x0 + diag[8] * x1 + diag[13] * x2 + diag[18] * x3 + diag[23] * x4;  in PCApply_PBJacobi()
68       yy[5 * i + 4] = diag[4] * x0 + diag[9] * x1 + diag[14] * x2 + diag[19] * x3 + diag[24] * x4;  in PCApply_PBJacobi()
81 …yy[6 * i]     = diag[0] * x0 + diag[6] * x1 + diag[12] * x2 + diag[18] * x3 + diag[24] * x4 + diag…  in PCApply_PBJacobi()
82 …yy[6 * i + 1] = diag[1] * x0 + diag[7] * x1 + diag[13] * x2 + diag[19] * x3 + diag[25] * x4 + diag…  in PCApply_PBJacobi()
83 …yy[6 * i + 2] = diag[2] * x0 + diag[8] * x1 + diag[14] * x2 + diag[20] * x3 + diag[26] * x4 + diag…  in PCApply_PBJacobi()
84 …yy[6 * i + 3] = diag[3] * x0 + diag[9] * x1 + diag[15] * x2 + diag[21] * x3 + diag[27] * x4 + diag…  in PCApply_PBJacobi()
85 …yy[6 * i + 4] = diag[4] * x0 + diag[10] * x1 + diag[16] * x2 + diag[22] * x3 + diag[28] * x4 + dia…  in PCApply_PBJacobi()
86 …yy[6 * i + 5] = diag[5] * x0 + diag[11] * x1 + diag[17] * x2 + diag[23] * x3 + diag[29] * x4 + dia…  in PCApply_PBJacobi()
100 …yy[7 * i]     = diag[0] * x0 + diag[7] * x1 + diag[14] * x2 + diag[21] * x3 + diag[28] * x4 + diag…  in PCApply_PBJacobi()
101 …yy[7 * i + 1] = diag[1] * x0 + diag[8] * x1 + diag[15] * x2 + diag[22] * x3 + diag[29] * x4 + diag…  in PCApply_PBJacobi()
102 …yy[7 * i + 2] = diag[2] * x0 + diag[9] * x1 + diag[16] * x2 + diag[23] * x3 + diag[30] * x4 + diag…  in PCApply_PBJacobi()
103 …yy[7 * i + 3] = diag[3] * x0 + diag[10] * x1 + diag[17] * x2 + diag[24] * x3 + diag[31] * x4 + dia…  in PCApply_PBJacobi()
104 …yy[7 * i + 4] = diag[4] * x0 + diag[11] * x1 + diag[18] * x2 + diag[25] * x3 + diag[32] * x4 + dia…  in PCApply_PBJacobi()
105 …yy[7 * i + 5] = diag[5] * x0 + diag[12] * x1 + diag[19] * x2 + diag[26] * x3 + diag[33] * x4 + dia…  in PCApply_PBJacobi()
106 …yy[7 * i + 6] = diag[6] * x0 + diag[13] * x1 + diag[20] * x2 + diag[27] * x3 + diag[34] * x4 + dia…  in PCApply_PBJacobi()
115         yy[bs * i + ib] = rowsum;  in PCApply_PBJacobi()
121   PetscCall(VecRestoreArray(y, &yy));  in PCApply_PBJacobi()
133   PetscScalar       *yy, x0, x1, x2, x3, x4, x5, x6;  in PCApplyTranspose_PBJacobi()  local
138   PetscCall(VecGetArray(y, &yy));  in PCApplyTranspose_PBJacobi()
141     for (i = 0; i < m; i++) yy[i] = diag[i] * xx[i];  in PCApplyTranspose_PBJacobi()
147       yy[2 * i]     = diag[0] * x0 + diag[1] * x1;  in PCApplyTranspose_PBJacobi()
148       yy[2 * i + 1] = diag[2] * x0 + diag[3] * x1;  in PCApplyTranspose_PBJacobi()
158       yy[3 * i]     = diag[0] * x0 + diag[1] * x1 + diag[2] * x2;  in PCApplyTranspose_PBJacobi()
159       yy[3 * i + 1] = diag[3] * x0 + diag[4] * x1 + diag[5] * x2;  in PCApplyTranspose_PBJacobi()
160       yy[3 * i + 2] = diag[6] * x0 + diag[7] * x1 + diag[8] * x2;  in PCApplyTranspose_PBJacobi()
171       yy[4 * i]     = diag[0] * x0 + diag[1] * x1 + diag[2] * x2 + diag[3] * x3;  in PCApplyTranspose_PBJacobi()
172       yy[4 * i + 1] = diag[4] * x0 + diag[5] * x1 + diag[6] * x2 + diag[7] * x3;  in PCApplyTranspose_PBJacobi()
173       yy[4 * i + 2] = diag[8] * x0 + diag[9] * x1 + diag[10] * x2 + diag[11] * x3;  in PCApplyTranspose_PBJacobi()
174       yy[4 * i + 3] = diag[12] * x0 + diag[13] * x1 + diag[14] * x2 + diag[15] * x3;  in PCApplyTranspose_PBJacobi()
186       yy[5 * i]     = diag[0] * x0 + diag[1] * x1 + diag[2] * x2 + diag[3] * x3 + diag[4] * x4;  in PCApplyTranspose_PBJacobi()
187       yy[5 * i + 1] = diag[5] * x0 + diag[6] * x1 + diag[7] * x2 + diag[8] * x3 + diag[9] * x4;  in PCApplyTranspose_PBJacobi()
188       yy[5 * i + 2] = diag[10] * x0 + diag[11] * x1 + diag[12] * x2 + diag[13] * x3 + diag[14] * x4;  in PCApplyTranspose_PBJacobi()
189       yy[5 * i + 3] = diag[15] * x0 + diag[16] * x1 + diag[17] * x2 + diag[18] * x3 + diag[19] * x4;  in PCApplyTranspose_PBJacobi()
190       yy[5 * i + 4] = diag[20] * x0 + diag[21] * x1 + diag[22] * x2 + diag[23] * x3 + diag[24] * x4;  in PCApplyTranspose_PBJacobi()
203 …yy[6 * i]     = diag[0] * x0 + diag[1] * x1 + diag[2] * x2 + diag[3] * x3 + diag[4] * x4 + diag[5]…  in PCApplyTranspose_PBJacobi()
204 …yy[6 * i + 1] = diag[6] * x0 + diag[7] * x1 + diag[8] * x2 + diag[9] * x3 + diag[10] * x4 + diag[1…  in PCApplyTranspose_PBJacobi()
205 …yy[6 * i + 2] = diag[12] * x0 + diag[13] * x1 + diag[14] * x2 + diag[15] * x3 + diag[16] * x4 + di…  in PCApplyTranspose_PBJacobi()
206 …yy[6 * i + 3] = diag[18] * x0 + diag[19] * x1 + diag[20] * x2 + diag[21] * x3 + diag[22] * x4 + di…  in PCApplyTranspose_PBJacobi()
207 …yy[6 * i + 4] = diag[24] * x0 + diag[25] * x1 + diag[26] * x2 + diag[27] * x3 + diag[28] * x4 + di…  in PCApplyTranspose_PBJacobi()
208 …yy[6 * i + 5] = diag[30] * x0 + diag[31] * x1 + diag[32] * x2 + diag[33] * x3 + diag[34] * x4 + di…  in PCApplyTranspose_PBJacobi()
222 …yy[7 * i]     = diag[0] * x0 + diag[1] * x1 + diag[2] * x2 + diag[3] * x3 + diag[4] * x4 + diag[5]…  in PCApplyTranspose_PBJacobi()
223 …yy[7 * i + 1] = diag[7] * x0 + diag[8] * x1 + diag[9] * x2 + diag[10] * x3 + diag[11] * x4 + diag[…  in PCApplyTranspose_PBJacobi()
224 …yy[7 * i + 2] = diag[14] * x0 + diag[15] * x1 + diag[16] * x2 + diag[17] * x3 + diag[18] * x4 + di…  in PCApplyTranspose_PBJacobi()
225 …yy[7 * i + 3] = diag[21] * x0 + diag[22] * x1 + diag[23] * x2 + diag[24] * x3 + diag[25] * x4 + di…  in PCApplyTranspose_PBJacobi()
226 …yy[7 * i + 4] = diag[28] * x0 + diag[29] * x1 + diag[30] * x2 + diag[31] * x3 + diag[32] * x4 + di…  in PCApplyTranspose_PBJacobi()
227 …yy[7 * i + 5] = diag[35] * x0 + diag[36] * x1 + diag[37] * x2 + diag[38] * x3 + diag[39] * x4 + di…  in PCApplyTranspose_PBJacobi()
228 …yy[7 * i + 6] = diag[42] * x0 + diag[43] * x1 + diag[44] * x2 + diag[45] * x3 + diag[46] * x4 + di…  in PCApplyTranspose_PBJacobi()
237         yy[bs * i + ib] = rowsum;  in PCApplyTranspose_PBJacobi()
243   PetscCall(VecRestoreArray(y, &yy));  in PCApplyTranspose_PBJacobi()