Lines Matching refs:k
16 PetscInt i__2, i__3, kp1, j, k, l, ll, i, ipvt[2], k3; in PetscKernel_A_gets_inverse_A_2() local
28 k = 1; in PetscKernel_A_gets_inverse_A_2()
29 kp1 = k + 1; in PetscKernel_A_gets_inverse_A_2()
30 k3 = 2 * k; in PetscKernel_A_gets_inverse_A_2()
31 k4 = k3 + k; in PetscKernel_A_gets_inverse_A_2()
34 i__2 = 3 - k; in PetscKernel_A_gets_inverse_A_2()
45 l += k - 1; in PetscKernel_A_gets_inverse_A_2()
46 ipvt[k - 1] = l; in PetscKernel_A_gets_inverse_A_2()
50 …(allowzeropivot, PETSC_COMM_SELF, PETSC_ERR_MAT_LU_ZRPVT, "Zero pivot, row %" PetscInt_FMT, k - 1); in PetscKernel_A_gets_inverse_A_2()
51 PetscCall(PetscInfo(NULL, "Zero pivot, row %" PetscInt_FMT "\n", k - 1)); in PetscKernel_A_gets_inverse_A_2()
59 if (l != k) { in PetscKernel_A_gets_inverse_A_2()
67 i__2 = 2 - k; in PetscKernel_A_gets_inverse_A_2()
76 if (l != k) { in PetscKernel_A_gets_inverse_A_2()
77 a[l + j3] = a[k + j3]; in PetscKernel_A_gets_inverse_A_2()
78 a[k + j3] = stmp; in PetscKernel_A_gets_inverse_A_2()
81 i__3 = 2 - k; in PetscKernel_A_gets_inverse_A_2()
82 ay = &a[1 + k + j3]; in PetscKernel_A_gets_inverse_A_2()
95 for (k = 1; k <= 2; ++k) { in PetscKernel_A_gets_inverse_A_2()
96 k3 = 2 * k; in PetscKernel_A_gets_inverse_A_2()
97 k4 = k3 + k; in PetscKernel_A_gets_inverse_A_2()
100 i__2 = k - 1; in PetscKernel_A_gets_inverse_A_2()
103 kp1 = k + 1; in PetscKernel_A_gets_inverse_A_2()
108 stmp = a[k + j3]; in PetscKernel_A_gets_inverse_A_2()
109 a[k + j3] = 0.0; in PetscKernel_A_gets_inverse_A_2()
111 for (ll = 0; ll < k; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_2()
116 k = 1; in PetscKernel_A_gets_inverse_A_2()
117 k3 = 2 * k; in PetscKernel_A_gets_inverse_A_2()
118 kp1 = k + 1; in PetscKernel_A_gets_inverse_A_2()
131 l = ipvt[k - 1]; in PetscKernel_A_gets_inverse_A_2()
132 if (l != k) { in PetscKernel_A_gets_inverse_A_2()
148 PetscInt i__2, i__3, kp1, j, k, l, ll, i, ipvt[9], kb, k3; in PetscKernel_A_gets_inverse_A_9() local
159 for (k = 1; k <= 8; ++k) { in PetscKernel_A_gets_inverse_A_9()
160 kp1 = k + 1; in PetscKernel_A_gets_inverse_A_9()
161 k3 = 9 * k; in PetscKernel_A_gets_inverse_A_9()
162 k4 = k3 + k; in PetscKernel_A_gets_inverse_A_9()
165 i__2 = 10 - k; in PetscKernel_A_gets_inverse_A_9()
176 l += k - 1; in PetscKernel_A_gets_inverse_A_9()
177 ipvt[k - 1] = l; in PetscKernel_A_gets_inverse_A_9()
181 …(allowzeropivot, PETSC_COMM_SELF, PETSC_ERR_MAT_LU_ZRPVT, "Zero pivot, row %" PetscInt_FMT, k - 1); in PetscKernel_A_gets_inverse_A_9()
182 PetscCall(PetscInfo(NULL, "Zero pivot, row %" PetscInt_FMT "\n", k - 1)); in PetscKernel_A_gets_inverse_A_9()
190 if (l != k) { in PetscKernel_A_gets_inverse_A_9()
198 i__2 = 9 - k; in PetscKernel_A_gets_inverse_A_9()
207 if (l != k) { in PetscKernel_A_gets_inverse_A_9()
208 a[l + j3] = a[k + j3]; in PetscKernel_A_gets_inverse_A_9()
209 a[k + j3] = stmp; in PetscKernel_A_gets_inverse_A_9()
212 i__3 = 9 - k; in PetscKernel_A_gets_inverse_A_9()
213 ay = &a[1 + k + j3]; in PetscKernel_A_gets_inverse_A_9()
226 for (k = 1; k <= 9; ++k) { in PetscKernel_A_gets_inverse_A_9()
227 k3 = 9 * k; in PetscKernel_A_gets_inverse_A_9()
228 k4 = k3 + k; in PetscKernel_A_gets_inverse_A_9()
231 i__2 = k - 1; in PetscKernel_A_gets_inverse_A_9()
234 kp1 = k + 1; in PetscKernel_A_gets_inverse_A_9()
239 stmp = a[k + j3]; in PetscKernel_A_gets_inverse_A_9()
240 a[k + j3] = 0.0; in PetscKernel_A_gets_inverse_A_9()
242 for (ll = 0; ll < k; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_9()
248 k = 9 - kb; in PetscKernel_A_gets_inverse_A_9()
249 k3 = 9 * k; in PetscKernel_A_gets_inverse_A_9()
250 kp1 = k + 1; in PetscKernel_A_gets_inverse_A_9()
270 l = ipvt[k - 1]; in PetscKernel_A_gets_inverse_A_9()
271 if (l != k) { in PetscKernel_A_gets_inverse_A_9()
319 PetscInt i__2, i__3, kp1, j, k, l, ll, i, kb, k3; in PetscKernel_A_gets_inverse_A_15() local
330 for (k = 1; k <= 14; ++k) { in PetscKernel_A_gets_inverse_A_15()
331 kp1 = k + 1; in PetscKernel_A_gets_inverse_A_15()
332 k3 = 15 * k; in PetscKernel_A_gets_inverse_A_15()
333 k4 = k3 + k; in PetscKernel_A_gets_inverse_A_15()
336 i__2 = 16 - k; in PetscKernel_A_gets_inverse_A_15()
347 l += k - 1; in PetscKernel_A_gets_inverse_A_15()
348 ipvt[k - 1] = l; in PetscKernel_A_gets_inverse_A_15()
352 …(allowzeropivot, PETSC_COMM_SELF, PETSC_ERR_MAT_LU_ZRPVT, "Zero pivot, row %" PetscInt_FMT, k - 1); in PetscKernel_A_gets_inverse_A_15()
353 PetscCall(PetscInfo(NULL, "Zero pivot, row %" PetscInt_FMT "\n", k - 1)); in PetscKernel_A_gets_inverse_A_15()
361 if (l != k) { in PetscKernel_A_gets_inverse_A_15()
369 i__2 = 15 - k; in PetscKernel_A_gets_inverse_A_15()
378 if (l != k) { in PetscKernel_A_gets_inverse_A_15()
379 a[l + j3] = a[k + j3]; in PetscKernel_A_gets_inverse_A_15()
380 a[k + j3] = stmp; in PetscKernel_A_gets_inverse_A_15()
383 i__3 = 15 - k; in PetscKernel_A_gets_inverse_A_15()
384 ay = &a[1 + k + j3]; in PetscKernel_A_gets_inverse_A_15()
397 for (k = 1; k <= 15; ++k) { in PetscKernel_A_gets_inverse_A_15()
398 k3 = 15 * k; in PetscKernel_A_gets_inverse_A_15()
399 k4 = k3 + k; in PetscKernel_A_gets_inverse_A_15()
402 i__2 = k - 1; in PetscKernel_A_gets_inverse_A_15()
405 kp1 = k + 1; in PetscKernel_A_gets_inverse_A_15()
410 stmp = a[k + j3]; in PetscKernel_A_gets_inverse_A_15()
411 a[k + j3] = 0.0; in PetscKernel_A_gets_inverse_A_15()
413 for (ll = 0; ll < k; ll++) ay[ll] += stmp * ax[ll]; in PetscKernel_A_gets_inverse_A_15()
419 k = 15 - kb; in PetscKernel_A_gets_inverse_A_15()
420 k3 = 15 * k; in PetscKernel_A_gets_inverse_A_15()
421 kp1 = k + 1; in PetscKernel_A_gets_inverse_A_15()
447 l = ipvt[k - 1]; in PetscKernel_A_gets_inverse_A_15()
448 if (l != k) { in PetscKernel_A_gets_inverse_A_15()