Lines Matching refs:v
12 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_N_inplace() local
30 v = aa + bs2 * ai[i]; in MatSolve_SeqBAIJ_N_inplace()
36 PetscKernel_v_gets_v_minus_A_times_w(bs, s, v, t + bs * (*vi++)); in MatSolve_SeqBAIJ_N_inplace()
37 v += bs2; in MatSolve_SeqBAIJ_N_inplace()
43 v = aa + bs2 * (a->diag[i] + 1); in MatSolve_SeqBAIJ_N_inplace()
48 PetscKernel_v_gets_v_minus_A_times_w(bs, ls, v, t + bs * (*vi++)); in MatSolve_SeqBAIJ_N_inplace()
49 v += bs2; in MatSolve_SeqBAIJ_N_inplace()
70 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_7_inplace() local
95 v = aa + 49 * ai[i]; in MatSolve_SeqBAIJ_7_inplace()
115 s1 -= v[0] * x1 + v[7] * x2 + v[14] * x3 + v[21] * x4 + v[28] * x5 + v[35] * x6 + v[42] * x7; in MatSolve_SeqBAIJ_7_inplace()
116 s2 -= v[1] * x1 + v[8] * x2 + v[15] * x3 + v[22] * x4 + v[29] * x5 + v[36] * x6 + v[43] * x7; in MatSolve_SeqBAIJ_7_inplace()
117 s3 -= v[2] * x1 + v[9] * x2 + v[16] * x3 + v[23] * x4 + v[30] * x5 + v[37] * x6 + v[44] * x7; in MatSolve_SeqBAIJ_7_inplace()
118 s4 -= v[3] * x1 + v[10] * x2 + v[17] * x3 + v[24] * x4 + v[31] * x5 + v[38] * x6 + v[45] * x7; in MatSolve_SeqBAIJ_7_inplace()
119 s5 -= v[4] * x1 + v[11] * x2 + v[18] * x3 + v[25] * x4 + v[32] * x5 + v[39] * x6 + v[46] * x7; in MatSolve_SeqBAIJ_7_inplace()
120 s6 -= v[5] * x1 + v[12] * x2 + v[19] * x3 + v[26] * x4 + v[33] * x5 + v[40] * x6 + v[47] * x7; in MatSolve_SeqBAIJ_7_inplace()
121 s7 -= v[6] * x1 + v[13] * x2 + v[20] * x3 + v[27] * x4 + v[34] * x5 + v[41] * x6 + v[48] * x7; in MatSolve_SeqBAIJ_7_inplace()
122 v += 49; in MatSolve_SeqBAIJ_7_inplace()
135 v = aa + 49 * diag[i] + 49; in MatSolve_SeqBAIJ_7_inplace()
155 s1 -= v[0] * x1 + v[7] * x2 + v[14] * x3 + v[21] * x4 + v[28] * x5 + v[35] * x6 + v[42] * x7; in MatSolve_SeqBAIJ_7_inplace()
156 s2 -= v[1] * x1 + v[8] * x2 + v[15] * x3 + v[22] * x4 + v[29] * x5 + v[36] * x6 + v[43] * x7; in MatSolve_SeqBAIJ_7_inplace()
157 s3 -= v[2] * x1 + v[9] * x2 + v[16] * x3 + v[23] * x4 + v[30] * x5 + v[37] * x6 + v[44] * x7; in MatSolve_SeqBAIJ_7_inplace()
158 s4 -= v[3] * x1 + v[10] * x2 + v[17] * x3 + v[24] * x4 + v[31] * x5 + v[38] * x6 + v[45] * x7; in MatSolve_SeqBAIJ_7_inplace()
159 s5 -= v[4] * x1 + v[11] * x2 + v[18] * x3 + v[25] * x4 + v[32] * x5 + v[39] * x6 + v[46] * x7; in MatSolve_SeqBAIJ_7_inplace()
160 s6 -= v[5] * x1 + v[12] * x2 + v[19] * x3 + v[26] * x4 + v[33] * x5 + v[40] * x6 + v[47] * x7; in MatSolve_SeqBAIJ_7_inplace()
161 s7 -= v[6] * x1 + v[13] * x2 + v[20] * x3 + v[27] * x4 + v[34] * x5 + v[41] * x6 + v[48] * x7; in MatSolve_SeqBAIJ_7_inplace()
162 v += 49; in MatSolve_SeqBAIJ_7_inplace()
165 v = aa + 49 * diag[i]; in MatSolve_SeqBAIJ_7_inplace()
166 …x[idc] = t[idt] = v[0] * s1 + v[7] * s2 + v[14] * s3 + v[21] * s4 + v[28] * s5 + v[35] * s6 + v[42… in MatSolve_SeqBAIJ_7_inplace()
167 …x[1 + idc] = t[1 + idt] = v[1] * s1 + v[8] * s2 + v[15] * s3 + v[22] * s4 + v[29] * s5 + v[36] * s… in MatSolve_SeqBAIJ_7_inplace()
168 …x[2 + idc] = t[2 + idt] = v[2] * s1 + v[9] * s2 + v[16] * s3 + v[23] * s4 + v[30] * s5 + v[37] * s… in MatSolve_SeqBAIJ_7_inplace()
169 …x[3 + idc] = t[3 + idt] = v[3] * s1 + v[10] * s2 + v[17] * s3 + v[24] * s4 + v[31] * s5 + v[38] * … in MatSolve_SeqBAIJ_7_inplace()
170 …x[4 + idc] = t[4 + idt] = v[4] * s1 + v[11] * s2 + v[18] * s3 + v[25] * s4 + v[32] * s5 + v[39] * … in MatSolve_SeqBAIJ_7_inplace()
171 …x[5 + idc] = t[5 + idt] = v[5] * s1 + v[12] * s2 + v[19] * s3 + v[26] * s4 + v[33] * s5 + v[40] * … in MatSolve_SeqBAIJ_7_inplace()
172 …x[6 + idc] = t[6 + idt] = v[6] * s1 + v[13] * s2 + v[20] * s3 + v[27] * s4 + v[34] * s5 + v[41] * … in MatSolve_SeqBAIJ_7_inplace()
190 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_7() local
215 v = aa + 49 * ai[i]; in MatSolve_SeqBAIJ_7()
235 s1 -= v[0] * x1 + v[7] * x2 + v[14] * x3 + v[21] * x4 + v[28] * x5 + v[35] * x6 + v[42] * x7; in MatSolve_SeqBAIJ_7()
236 s2 -= v[1] * x1 + v[8] * x2 + v[15] * x3 + v[22] * x4 + v[29] * x5 + v[36] * x6 + v[43] * x7; in MatSolve_SeqBAIJ_7()
237 s3 -= v[2] * x1 + v[9] * x2 + v[16] * x3 + v[23] * x4 + v[30] * x5 + v[37] * x6 + v[44] * x7; in MatSolve_SeqBAIJ_7()
238 s4 -= v[3] * x1 + v[10] * x2 + v[17] * x3 + v[24] * x4 + v[31] * x5 + v[38] * x6 + v[45] * x7; in MatSolve_SeqBAIJ_7()
239 s5 -= v[4] * x1 + v[11] * x2 + v[18] * x3 + v[25] * x4 + v[32] * x5 + v[39] * x6 + v[46] * x7; in MatSolve_SeqBAIJ_7()
240 s6 -= v[5] * x1 + v[12] * x2 + v[19] * x3 + v[26] * x4 + v[33] * x5 + v[40] * x6 + v[47] * x7; in MatSolve_SeqBAIJ_7()
241 s7 -= v[6] * x1 + v[13] * x2 + v[20] * x3 + v[27] * x4 + v[34] * x5 + v[41] * x6 + v[48] * x7; in MatSolve_SeqBAIJ_7()
242 v += 49; in MatSolve_SeqBAIJ_7()
255 v = aa + 49 * (adiag[i + 1] + 1); in MatSolve_SeqBAIJ_7()
275 s1 -= v[0] * x1 + v[7] * x2 + v[14] * x3 + v[21] * x4 + v[28] * x5 + v[35] * x6 + v[42] * x7; in MatSolve_SeqBAIJ_7()
276 s2 -= v[1] * x1 + v[8] * x2 + v[15] * x3 + v[22] * x4 + v[29] * x5 + v[36] * x6 + v[43] * x7; in MatSolve_SeqBAIJ_7()
277 s3 -= v[2] * x1 + v[9] * x2 + v[16] * x3 + v[23] * x4 + v[30] * x5 + v[37] * x6 + v[44] * x7; in MatSolve_SeqBAIJ_7()
278 s4 -= v[3] * x1 + v[10] * x2 + v[17] * x3 + v[24] * x4 + v[31] * x5 + v[38] * x6 + v[45] * x7; in MatSolve_SeqBAIJ_7()
279 s5 -= v[4] * x1 + v[11] * x2 + v[18] * x3 + v[25] * x4 + v[32] * x5 + v[39] * x6 + v[46] * x7; in MatSolve_SeqBAIJ_7()
280 s6 -= v[5] * x1 + v[12] * x2 + v[19] * x3 + v[26] * x4 + v[33] * x5 + v[40] * x6 + v[47] * x7; in MatSolve_SeqBAIJ_7()
281 s7 -= v[6] * x1 + v[13] * x2 + v[20] * x3 + v[27] * x4 + v[34] * x5 + v[41] * x6 + v[48] * x7; in MatSolve_SeqBAIJ_7()
282 v += 49; in MatSolve_SeqBAIJ_7()
285 …x[idc] = t[idt] = v[0] * s1 + v[7] * s2 + v[14] * s3 + v[21] * s4 + v[28] * s5 + v[35] * s6 + v[42… in MatSolve_SeqBAIJ_7()
286 …x[1 + idc] = t[1 + idt] = v[1] * s1 + v[8] * s2 + v[15] * s3 + v[22] * s4 + v[29] * s5 + v[36] * s… in MatSolve_SeqBAIJ_7()
287 …x[2 + idc] = t[2 + idt] = v[2] * s1 + v[9] * s2 + v[16] * s3 + v[23] * s4 + v[30] * s5 + v[37] * s… in MatSolve_SeqBAIJ_7()
288 …x[3 + idc] = t[3 + idt] = v[3] * s1 + v[10] * s2 + v[17] * s3 + v[24] * s4 + v[31] * s5 + v[38] * … in MatSolve_SeqBAIJ_7()
289 …x[4 + idc] = t[4 + idt] = v[4] * s1 + v[11] * s2 + v[18] * s3 + v[25] * s4 + v[32] * s5 + v[39] * … in MatSolve_SeqBAIJ_7()
290 …x[5 + idc] = t[5 + idt] = v[5] * s1 + v[12] * s2 + v[19] * s3 + v[26] * s4 + v[33] * s5 + v[40] * … in MatSolve_SeqBAIJ_7()
291 …x[6 + idc] = t[6 + idt] = v[6] * s1 + v[13] * s2 + v[20] * s3 + v[27] * s4 + v[34] * s5 + v[41] * … in MatSolve_SeqBAIJ_7()
309 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_6_inplace() local
332 v = aa + 36 * ai[i]; in MatSolve_SeqBAIJ_6_inplace()
350 s1 -= v[0] * x1 + v[6] * x2 + v[12] * x3 + v[18] * x4 + v[24] * x5 + v[30] * x6; in MatSolve_SeqBAIJ_6_inplace()
351 s2 -= v[1] * x1 + v[7] * x2 + v[13] * x3 + v[19] * x4 + v[25] * x5 + v[31] * x6; in MatSolve_SeqBAIJ_6_inplace()
352 s3 -= v[2] * x1 + v[8] * x2 + v[14] * x3 + v[20] * x4 + v[26] * x5 + v[32] * x6; in MatSolve_SeqBAIJ_6_inplace()
353 s4 -= v[3] * x1 + v[9] * x2 + v[15] * x3 + v[21] * x4 + v[27] * x5 + v[33] * x6; in MatSolve_SeqBAIJ_6_inplace()
354 s5 -= v[4] * x1 + v[10] * x2 + v[16] * x3 + v[22] * x4 + v[28] * x5 + v[34] * x6; in MatSolve_SeqBAIJ_6_inplace()
355 s6 -= v[5] * x1 + v[11] * x2 + v[17] * x3 + v[23] * x4 + v[29] * x5 + v[35] * x6; in MatSolve_SeqBAIJ_6_inplace()
356 v += 36; in MatSolve_SeqBAIJ_6_inplace()
368 v = aa + 36 * diag[i] + 36; in MatSolve_SeqBAIJ_6_inplace()
386 s1 -= v[0] * x1 + v[6] * x2 + v[12] * x3 + v[18] * x4 + v[24] * x5 + v[30] * x6; in MatSolve_SeqBAIJ_6_inplace()
387 s2 -= v[1] * x1 + v[7] * x2 + v[13] * x3 + v[19] * x4 + v[25] * x5 + v[31] * x6; in MatSolve_SeqBAIJ_6_inplace()
388 s3 -= v[2] * x1 + v[8] * x2 + v[14] * x3 + v[20] * x4 + v[26] * x5 + v[32] * x6; in MatSolve_SeqBAIJ_6_inplace()
389 s4 -= v[3] * x1 + v[9] * x2 + v[15] * x3 + v[21] * x4 + v[27] * x5 + v[33] * x6; in MatSolve_SeqBAIJ_6_inplace()
390 s5 -= v[4] * x1 + v[10] * x2 + v[16] * x3 + v[22] * x4 + v[28] * x5 + v[34] * x6; in MatSolve_SeqBAIJ_6_inplace()
391 s6 -= v[5] * x1 + v[11] * x2 + v[17] * x3 + v[23] * x4 + v[29] * x5 + v[35] * x6; in MatSolve_SeqBAIJ_6_inplace()
392 v += 36; in MatSolve_SeqBAIJ_6_inplace()
395 v = aa + 36 * diag[i]; in MatSolve_SeqBAIJ_6_inplace()
396 x[idc] = t[idt] = v[0] * s1 + v[6] * s2 + v[12] * s3 + v[18] * s4 + v[24] * s5 + v[30] * s6; in MatSolve_SeqBAIJ_6_inplace()
397 …x[1 + idc] = t[1 + idt] = v[1] * s1 + v[7] * s2 + v[13] * s3 + v[19] * s4 + v[25] * s5 + v[31] * s… in MatSolve_SeqBAIJ_6_inplace()
398 …x[2 + idc] = t[2 + idt] = v[2] * s1 + v[8] * s2 + v[14] * s3 + v[20] * s4 + v[26] * s5 + v[32] * s… in MatSolve_SeqBAIJ_6_inplace()
399 …x[3 + idc] = t[3 + idt] = v[3] * s1 + v[9] * s2 + v[15] * s3 + v[21] * s4 + v[27] * s5 + v[33] * s… in MatSolve_SeqBAIJ_6_inplace()
400 …x[4 + idc] = t[4 + idt] = v[4] * s1 + v[10] * s2 + v[16] * s3 + v[22] * s4 + v[28] * s5 + v[34] * … in MatSolve_SeqBAIJ_6_inplace()
401 …x[5 + idc] = t[5 + idt] = v[5] * s1 + v[11] * s2 + v[17] * s3 + v[23] * s4 + v[29] * s5 + v[35] * … in MatSolve_SeqBAIJ_6_inplace()
419 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_6() local
442 v = aa + 36 * ai[i]; in MatSolve_SeqBAIJ_6()
460 s1 -= v[0] * x1 + v[6] * x2 + v[12] * x3 + v[18] * x4 + v[24] * x5 + v[30] * x6; in MatSolve_SeqBAIJ_6()
461 s2 -= v[1] * x1 + v[7] * x2 + v[13] * x3 + v[19] * x4 + v[25] * x5 + v[31] * x6; in MatSolve_SeqBAIJ_6()
462 s3 -= v[2] * x1 + v[8] * x2 + v[14] * x3 + v[20] * x4 + v[26] * x5 + v[32] * x6; in MatSolve_SeqBAIJ_6()
463 s4 -= v[3] * x1 + v[9] * x2 + v[15] * x3 + v[21] * x4 + v[27] * x5 + v[33] * x6; in MatSolve_SeqBAIJ_6()
464 s5 -= v[4] * x1 + v[10] * x2 + v[16] * x3 + v[22] * x4 + v[28] * x5 + v[34] * x6; in MatSolve_SeqBAIJ_6()
465 s6 -= v[5] * x1 + v[11] * x2 + v[17] * x3 + v[23] * x4 + v[29] * x5 + v[35] * x6; in MatSolve_SeqBAIJ_6()
466 v += 36; in MatSolve_SeqBAIJ_6()
478 v = aa + 36 * (adiag[i + 1] + 1); in MatSolve_SeqBAIJ_6()
496 s1 -= v[0] * x1 + v[6] * x2 + v[12] * x3 + v[18] * x4 + v[24] * x5 + v[30] * x6; in MatSolve_SeqBAIJ_6()
497 s2 -= v[1] * x1 + v[7] * x2 + v[13] * x3 + v[19] * x4 + v[25] * x5 + v[31] * x6; in MatSolve_SeqBAIJ_6()
498 s3 -= v[2] * x1 + v[8] * x2 + v[14] * x3 + v[20] * x4 + v[26] * x5 + v[32] * x6; in MatSolve_SeqBAIJ_6()
499 s4 -= v[3] * x1 + v[9] * x2 + v[15] * x3 + v[21] * x4 + v[27] * x5 + v[33] * x6; in MatSolve_SeqBAIJ_6()
500 s5 -= v[4] * x1 + v[10] * x2 + v[16] * x3 + v[22] * x4 + v[28] * x5 + v[34] * x6; in MatSolve_SeqBAIJ_6()
501 s6 -= v[5] * x1 + v[11] * x2 + v[17] * x3 + v[23] * x4 + v[29] * x5 + v[35] * x6; in MatSolve_SeqBAIJ_6()
502 v += 36; in MatSolve_SeqBAIJ_6()
505 x[idc] = t[idt] = v[0] * s1 + v[6] * s2 + v[12] * s3 + v[18] * s4 + v[24] * s5 + v[30] * s6; in MatSolve_SeqBAIJ_6()
506 …x[1 + idc] = t[1 + idt] = v[1] * s1 + v[7] * s2 + v[13] * s3 + v[19] * s4 + v[25] * s5 + v[31] * s… in MatSolve_SeqBAIJ_6()
507 …x[2 + idc] = t[2 + idt] = v[2] * s1 + v[8] * s2 + v[14] * s3 + v[20] * s4 + v[26] * s5 + v[32] * s… in MatSolve_SeqBAIJ_6()
508 …x[3 + idc] = t[3 + idt] = v[3] * s1 + v[9] * s2 + v[15] * s3 + v[21] * s4 + v[27] * s5 + v[33] * s… in MatSolve_SeqBAIJ_6()
509 …x[4 + idc] = t[4 + idt] = v[4] * s1 + v[10] * s2 + v[16] * s3 + v[22] * s4 + v[28] * s5 + v[34] * … in MatSolve_SeqBAIJ_6()
510 …x[5 + idc] = t[5 + idt] = v[5] * s1 + v[11] * s2 + v[17] * s3 + v[23] * s4 + v[29] * s5 + v[35] * … in MatSolve_SeqBAIJ_6()
528 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_5_inplace() local
550 v = aa + 25 * ai[i]; in MatSolve_SeqBAIJ_5_inplace()
566 s1 -= v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatSolve_SeqBAIJ_5_inplace()
567 s2 -= v[1] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatSolve_SeqBAIJ_5_inplace()
568 s3 -= v[2] * x1 + v[7] * x2 + v[12] * x3 + v[17] * x4 + v[22] * x5; in MatSolve_SeqBAIJ_5_inplace()
569 s4 -= v[3] * x1 + v[8] * x2 + v[13] * x3 + v[18] * x4 + v[23] * x5; in MatSolve_SeqBAIJ_5_inplace()
570 s5 -= v[4] * x1 + v[9] * x2 + v[14] * x3 + v[19] * x4 + v[24] * x5; in MatSolve_SeqBAIJ_5_inplace()
571 v += 25; in MatSolve_SeqBAIJ_5_inplace()
582 v = aa + 25 * diag[i] + 25; in MatSolve_SeqBAIJ_5_inplace()
598 s1 -= v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatSolve_SeqBAIJ_5_inplace()
599 s2 -= v[1] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatSolve_SeqBAIJ_5_inplace()
600 s3 -= v[2] * x1 + v[7] * x2 + v[12] * x3 + v[17] * x4 + v[22] * x5; in MatSolve_SeqBAIJ_5_inplace()
601 s4 -= v[3] * x1 + v[8] * x2 + v[13] * x3 + v[18] * x4 + v[23] * x5; in MatSolve_SeqBAIJ_5_inplace()
602 s5 -= v[4] * x1 + v[9] * x2 + v[14] * x3 + v[19] * x4 + v[24] * x5; in MatSolve_SeqBAIJ_5_inplace()
603 v += 25; in MatSolve_SeqBAIJ_5_inplace()
606 v = aa + 25 * diag[i]; in MatSolve_SeqBAIJ_5_inplace()
607 x[idc] = t[idt] = v[0] * s1 + v[5] * s2 + v[10] * s3 + v[15] * s4 + v[20] * s5; in MatSolve_SeqBAIJ_5_inplace()
608 x[1 + idc] = t[1 + idt] = v[1] * s1 + v[6] * s2 + v[11] * s3 + v[16] * s4 + v[21] * s5; in MatSolve_SeqBAIJ_5_inplace()
609 x[2 + idc] = t[2 + idt] = v[2] * s1 + v[7] * s2 + v[12] * s3 + v[17] * s4 + v[22] * s5; in MatSolve_SeqBAIJ_5_inplace()
610 x[3 + idc] = t[3 + idt] = v[3] * s1 + v[8] * s2 + v[13] * s3 + v[18] * s4 + v[23] * s5; in MatSolve_SeqBAIJ_5_inplace()
611 x[4 + idc] = t[4 + idt] = v[4] * s1 + v[9] * s2 + v[14] * s3 + v[19] * s4 + v[24] * s5; in MatSolve_SeqBAIJ_5_inplace()
629 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_5() local
651 v = aa + 25 * ai[i]; in MatSolve_SeqBAIJ_5()
667 s1 -= v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatSolve_SeqBAIJ_5()
668 s2 -= v[1] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatSolve_SeqBAIJ_5()
669 s3 -= v[2] * x1 + v[7] * x2 + v[12] * x3 + v[17] * x4 + v[22] * x5; in MatSolve_SeqBAIJ_5()
670 s4 -= v[3] * x1 + v[8] * x2 + v[13] * x3 + v[18] * x4 + v[23] * x5; in MatSolve_SeqBAIJ_5()
671 s5 -= v[4] * x1 + v[9] * x2 + v[14] * x3 + v[19] * x4 + v[24] * x5; in MatSolve_SeqBAIJ_5()
672 v += 25; in MatSolve_SeqBAIJ_5()
683 v = aa + 25 * (adiag[i + 1] + 1); in MatSolve_SeqBAIJ_5()
699 s1 -= v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatSolve_SeqBAIJ_5()
700 s2 -= v[1] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatSolve_SeqBAIJ_5()
701 s3 -= v[2] * x1 + v[7] * x2 + v[12] * x3 + v[17] * x4 + v[22] * x5; in MatSolve_SeqBAIJ_5()
702 s4 -= v[3] * x1 + v[8] * x2 + v[13] * x3 + v[18] * x4 + v[23] * x5; in MatSolve_SeqBAIJ_5()
703 s5 -= v[4] * x1 + v[9] * x2 + v[14] * x3 + v[19] * x4 + v[24] * x5; in MatSolve_SeqBAIJ_5()
704 v += 25; in MatSolve_SeqBAIJ_5()
707 x[idc] = t[idt] = v[0] * s1 + v[5] * s2 + v[10] * s3 + v[15] * s4 + v[20] * s5; in MatSolve_SeqBAIJ_5()
708 x[1 + idc] = t[1 + idt] = v[1] * s1 + v[6] * s2 + v[11] * s3 + v[16] * s4 + v[21] * s5; in MatSolve_SeqBAIJ_5()
709 x[2 + idc] = t[2 + idt] = v[2] * s1 + v[7] * s2 + v[12] * s3 + v[17] * s4 + v[22] * s5; in MatSolve_SeqBAIJ_5()
710 x[3 + idc] = t[3 + idt] = v[3] * s1 + v[8] * s2 + v[13] * s3 + v[18] * s4 + v[23] * s5; in MatSolve_SeqBAIJ_5()
711 x[4 + idc] = t[4 + idt] = v[4] * s1 + v[9] * s2 + v[14] * s3 + v[19] * s4 + v[24] * s5; in MatSolve_SeqBAIJ_5()
729 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_4_inplace() local
750 v = aa + 16 * ai[i]; in MatSolve_SeqBAIJ_4_inplace()
764 s1 -= v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatSolve_SeqBAIJ_4_inplace()
765 s2 -= v[1] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatSolve_SeqBAIJ_4_inplace()
766 s3 -= v[2] * x1 + v[6] * x2 + v[10] * x3 + v[14] * x4; in MatSolve_SeqBAIJ_4_inplace()
767 s4 -= v[3] * x1 + v[7] * x2 + v[11] * x3 + v[15] * x4; in MatSolve_SeqBAIJ_4_inplace()
768 v += 16; in MatSolve_SeqBAIJ_4_inplace()
778 v = aa + 16 * diag[i] + 16; in MatSolve_SeqBAIJ_4_inplace()
792 s1 -= v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatSolve_SeqBAIJ_4_inplace()
793 s2 -= v[1] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatSolve_SeqBAIJ_4_inplace()
794 s3 -= v[2] * x1 + v[6] * x2 + v[10] * x3 + v[14] * x4; in MatSolve_SeqBAIJ_4_inplace()
795 s4 -= v[3] * x1 + v[7] * x2 + v[11] * x3 + v[15] * x4; in MatSolve_SeqBAIJ_4_inplace()
796 v += 16; in MatSolve_SeqBAIJ_4_inplace()
799 v = aa + 16 * diag[i]; in MatSolve_SeqBAIJ_4_inplace()
800 x[idc] = t[idt] = v[0] * s1 + v[4] * s2 + v[8] * s3 + v[12] * s4; in MatSolve_SeqBAIJ_4_inplace()
801 x[1 + idc] = t[1 + idt] = v[1] * s1 + v[5] * s2 + v[9] * s3 + v[13] * s4; in MatSolve_SeqBAIJ_4_inplace()
802 x[2 + idc] = t[2 + idt] = v[2] * s1 + v[6] * s2 + v[10] * s3 + v[14] * s4; in MatSolve_SeqBAIJ_4_inplace()
803 x[3 + idc] = t[3 + idt] = v[3] * s1 + v[7] * s2 + v[11] * s3 + v[15] * s4; in MatSolve_SeqBAIJ_4_inplace()
821 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_4() local
842 v = aa + 16 * ai[i]; in MatSolve_SeqBAIJ_4()
856 s1 -= v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatSolve_SeqBAIJ_4()
857 s2 -= v[1] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatSolve_SeqBAIJ_4()
858 s3 -= v[2] * x1 + v[6] * x2 + v[10] * x3 + v[14] * x4; in MatSolve_SeqBAIJ_4()
859 s4 -= v[3] * x1 + v[7] * x2 + v[11] * x3 + v[15] * x4; in MatSolve_SeqBAIJ_4()
860 v += 16; in MatSolve_SeqBAIJ_4()
870 v = aa + 16 * (adiag[i + 1] + 1); in MatSolve_SeqBAIJ_4()
884 s1 -= v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatSolve_SeqBAIJ_4()
885 s2 -= v[1] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatSolve_SeqBAIJ_4()
886 s3 -= v[2] * x1 + v[6] * x2 + v[10] * x3 + v[14] * x4; in MatSolve_SeqBAIJ_4()
887 s4 -= v[3] * x1 + v[7] * x2 + v[11] * x3 + v[15] * x4; in MatSolve_SeqBAIJ_4()
888 v += 16; in MatSolve_SeqBAIJ_4()
891 x[idc] = t[idt] = v[0] * s1 + v[4] * s2 + v[8] * s3 + v[12] * s4; in MatSolve_SeqBAIJ_4()
892 x[1 + idc] = t[1 + idt] = v[1] * s1 + v[5] * s2 + v[9] * s3 + v[13] * s4; in MatSolve_SeqBAIJ_4()
893 x[2 + idc] = t[2 + idt] = v[2] * s1 + v[6] * s2 + v[10] * s3 + v[14] * s4; in MatSolve_SeqBAIJ_4()
894 x[3 + idc] = t[3 + idt] = v[3] * s1 + v[7] * s2 + v[11] * s3 + v[15] * s4; in MatSolve_SeqBAIJ_4()
912 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_3_inplace() local
932 v = aa + 9 * ai[i]; in MatSolve_SeqBAIJ_3_inplace()
944 s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3; in MatSolve_SeqBAIJ_3_inplace()
945 s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3; in MatSolve_SeqBAIJ_3_inplace()
946 s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3; in MatSolve_SeqBAIJ_3_inplace()
947 v += 9; in MatSolve_SeqBAIJ_3_inplace()
956 v = aa + 9 * diag[i] + 9; in MatSolve_SeqBAIJ_3_inplace()
968 s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3; in MatSolve_SeqBAIJ_3_inplace()
969 s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3; in MatSolve_SeqBAIJ_3_inplace()
970 s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3; in MatSolve_SeqBAIJ_3_inplace()
971 v += 9; in MatSolve_SeqBAIJ_3_inplace()
974 v = aa + 9 * diag[i]; in MatSolve_SeqBAIJ_3_inplace()
975 x[idc] = t[idt] = v[0] * s1 + v[3] * s2 + v[6] * s3; in MatSolve_SeqBAIJ_3_inplace()
976 x[1 + idc] = t[1 + idt] = v[1] * s1 + v[4] * s2 + v[7] * s3; in MatSolve_SeqBAIJ_3_inplace()
977 x[2 + idc] = t[2 + idt] = v[2] * s1 + v[5] * s2 + v[8] * s3; in MatSolve_SeqBAIJ_3_inplace()
994 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_3() local
1014 v = aa + 9 * ai[i]; in MatSolve_SeqBAIJ_3()
1026 s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3; in MatSolve_SeqBAIJ_3()
1027 s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3; in MatSolve_SeqBAIJ_3()
1028 s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3; in MatSolve_SeqBAIJ_3()
1029 v += 9; in MatSolve_SeqBAIJ_3()
1038 v = aa + 9 * (adiag[i + 1] + 1); in MatSolve_SeqBAIJ_3()
1050 s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3; in MatSolve_SeqBAIJ_3()
1051 s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3; in MatSolve_SeqBAIJ_3()
1052 s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3; in MatSolve_SeqBAIJ_3()
1053 v += 9; in MatSolve_SeqBAIJ_3()
1056 x[idc] = t[idt] = v[0] * s1 + v[3] * s2 + v[6] * s3; in MatSolve_SeqBAIJ_3()
1057 x[1 + idc] = t[1 + idt] = v[1] * s1 + v[4] * s2 + v[7] * s3; in MatSolve_SeqBAIJ_3()
1058 x[2 + idc] = t[2 + idt] = v[2] * s1 + v[5] * s2 + v[8] * s3; in MatSolve_SeqBAIJ_3()
1075 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_2_inplace() local
1094 v = aa + 4 * ai[i]; in MatSolve_SeqBAIJ_2_inplace()
1104 s1 -= v[0] * x1 + v[2] * x2; in MatSolve_SeqBAIJ_2_inplace()
1105 s2 -= v[1] * x1 + v[3] * x2; in MatSolve_SeqBAIJ_2_inplace()
1106 v += 4; in MatSolve_SeqBAIJ_2_inplace()
1114 v = aa + 4 * diag[i] + 4; in MatSolve_SeqBAIJ_2_inplace()
1124 s1 -= v[0] * x1 + v[2] * x2; in MatSolve_SeqBAIJ_2_inplace()
1125 s2 -= v[1] * x1 + v[3] * x2; in MatSolve_SeqBAIJ_2_inplace()
1126 v += 4; in MatSolve_SeqBAIJ_2_inplace()
1129 v = aa + 4 * diag[i]; in MatSolve_SeqBAIJ_2_inplace()
1130 x[idc] = t[idt] = v[0] * s1 + v[2] * s2; in MatSolve_SeqBAIJ_2_inplace()
1131 x[1 + idc] = t[1 + idt] = v[1] * s1 + v[3] * s2; in MatSolve_SeqBAIJ_2_inplace()
1148 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_2() local
1167 v = aa + 4 * ai[i]; in MatSolve_SeqBAIJ_2()
1177 s1 -= v[0] * x1 + v[2] * x2; in MatSolve_SeqBAIJ_2()
1178 s2 -= v[1] * x1 + v[3] * x2; in MatSolve_SeqBAIJ_2()
1179 v += 4; in MatSolve_SeqBAIJ_2()
1187 v = aa + 4 * (adiag[i + 1] + 1); in MatSolve_SeqBAIJ_2()
1197 s1 -= v[0] * x1 + v[2] * x2; in MatSolve_SeqBAIJ_2()
1198 s2 -= v[1] * x1 + v[3] * x2; in MatSolve_SeqBAIJ_2()
1199 v += 4; in MatSolve_SeqBAIJ_2()
1202 x[idc] = t[idt] = v[0] * s1 + v[2] * s2; in MatSolve_SeqBAIJ_2()
1203 x[1 + idc] = t[1 + idt] = v[1] * s1 + v[3] * s2; in MatSolve_SeqBAIJ_2()
1220 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_1_inplace() local
1239 v = aa + ai[i]; in MatSolve_SeqBAIJ_1_inplace()
1243 while (nz--) s1 -= (*v++) * t[*vi++]; in MatSolve_SeqBAIJ_1_inplace()
1248 v = aa + diag[i] + 1; in MatSolve_SeqBAIJ_1_inplace()
1252 while (nz--) s1 -= (*v++) * t[*vi++]; in MatSolve_SeqBAIJ_1_inplace()
1272 const MatScalar *aa = a->a, *v; in MatSolve_SeqBAIJ_1() local
1288 v = aa; in MatSolve_SeqBAIJ_1()
1293 PetscSparseDenseMinusDot(sum, tmp, v, vi, nz); in MatSolve_SeqBAIJ_1()
1295 v += nz; in MatSolve_SeqBAIJ_1()
1301 v = aa + adiag[i + 1] + 1; in MatSolve_SeqBAIJ_1()
1305 PetscSparseDenseMinusDot(sum, tmp, v, vi, nz); in MatSolve_SeqBAIJ_1()
1306 x[c[i]] = tmp[i] = sum * v[nz]; /* v[nz] = aa[adiag[i]] */ in MatSolve_SeqBAIJ_1()