Lines Matching refs:tp
222 PetscScalar *x, x0, x1, x2, x3, x4, x5, x6, *t, *tp; in MatSolve_SeqSBAIJ_7_inplace() local
232 tp = t; in MatSolve_SeqSBAIJ_7_inplace()
235 tp[0] = b[idx]; in MatSolve_SeqSBAIJ_7_inplace()
236 tp[1] = b[idx + 1]; in MatSolve_SeqSBAIJ_7_inplace()
237 tp[2] = b[idx + 2]; in MatSolve_SeqSBAIJ_7_inplace()
238 tp[3] = b[idx + 3]; in MatSolve_SeqSBAIJ_7_inplace()
239 tp[4] = b[idx + 4]; in MatSolve_SeqSBAIJ_7_inplace()
240 tp[5] = b[idx + 5]; in MatSolve_SeqSBAIJ_7_inplace()
241 tp[6] = b[idx + 6]; in MatSolve_SeqSBAIJ_7_inplace()
242 tp += 7; in MatSolve_SeqSBAIJ_7_inplace()
248 tp = t + k * 7; in MatSolve_SeqSBAIJ_7_inplace()
249 x0 = tp[0]; in MatSolve_SeqSBAIJ_7_inplace()
250 x1 = tp[1]; in MatSolve_SeqSBAIJ_7_inplace()
251 x2 = tp[2]; in MatSolve_SeqSBAIJ_7_inplace()
252 x3 = tp[3]; in MatSolve_SeqSBAIJ_7_inplace()
253 x4 = tp[4]; in MatSolve_SeqSBAIJ_7_inplace()
254 x5 = tp[5]; in MatSolve_SeqSBAIJ_7_inplace()
255 x6 = tp[6]; in MatSolve_SeqSBAIJ_7_inplace()
257 tp = t + (*vj) * 7; in MatSolve_SeqSBAIJ_7_inplace()
259 tp[0] += v[0] * x0 + v[1] * x1 + v[2] * x2 + v[3] * x3 + v[4] * x4 + v[5] * x5 + v[6] * x6; in MatSolve_SeqSBAIJ_7_inplace()
260 … tp[1] += v[7] * x0 + v[8] * x1 + v[9] * x2 + v[10] * x3 + v[11] * x4 + v[12] * x5 + v[13] * x6; in MatSolve_SeqSBAIJ_7_inplace()
261 … tp[2] += v[14] * x0 + v[15] * x1 + v[16] * x2 + v[17] * x3 + v[18] * x4 + v[19] * x5 + v[20] * x6; in MatSolve_SeqSBAIJ_7_inplace()
262 … tp[3] += v[21] * x0 + v[22] * x1 + v[23] * x2 + v[24] * x3 + v[25] * x4 + v[26] * x5 + v[27] * x6; in MatSolve_SeqSBAIJ_7_inplace()
263 … tp[4] += v[28] * x0 + v[29] * x1 + v[30] * x2 + v[31] * x3 + v[32] * x4 + v[33] * x5 + v[34] * x6; in MatSolve_SeqSBAIJ_7_inplace()
264 … tp[5] += v[35] * x0 + v[36] * x1 + v[37] * x2 + v[38] * x3 + v[39] * x4 + v[40] * x5 + v[41] * x6; in MatSolve_SeqSBAIJ_7_inplace()
265 … tp[6] += v[42] * x0 + v[43] * x1 + v[44] * x2 + v[45] * x3 + v[46] * x4 + v[47] * x5 + v[48] * x6; in MatSolve_SeqSBAIJ_7_inplace()
267 tp = t + (*vj) * 7; in MatSolve_SeqSBAIJ_7_inplace()
273 tp = t + k * 7; in MatSolve_SeqSBAIJ_7_inplace()
274 tp[0] = d[0] * x0 + d[7] * x1 + d[14] * x2 + d[21] * x3 + d[28] * x4 + d[35] * x5 + d[42] * x6; in MatSolve_SeqSBAIJ_7_inplace()
275 tp[1] = d[1] * x0 + d[8] * x1 + d[15] * x2 + d[22] * x3 + d[29] * x4 + d[36] * x5 + d[43] * x6; in MatSolve_SeqSBAIJ_7_inplace()
276 tp[2] = d[2] * x0 + d[9] * x1 + d[16] * x2 + d[23] * x3 + d[30] * x4 + d[37] * x5 + d[44] * x6; in MatSolve_SeqSBAIJ_7_inplace()
277 tp[3] = d[3] * x0 + d[10] * x1 + d[17] * x2 + d[24] * x3 + d[31] * x4 + d[38] * x5 + d[45] * x6; in MatSolve_SeqSBAIJ_7_inplace()
278 tp[4] = d[4] * x0 + d[11] * x1 + d[18] * x2 + d[25] * x3 + d[32] * x4 + d[39] * x5 + d[46] * x6; in MatSolve_SeqSBAIJ_7_inplace()
279 tp[5] = d[5] * x0 + d[12] * x1 + d[19] * x2 + d[26] * x3 + d[33] * x4 + d[40] * x5 + d[47] * x6; in MatSolve_SeqSBAIJ_7_inplace()
280 tp[6] = d[6] * x0 + d[13] * x1 + d[20] * x2 + d[27] * x3 + d[34] * x4 + d[41] * x5 + d[48] * x6; in MatSolve_SeqSBAIJ_7_inplace()
287 tp = t + k * 7; in MatSolve_SeqSBAIJ_7_inplace()
288 x0 = tp[0]; in MatSolve_SeqSBAIJ_7_inplace()
289 x1 = tp[1]; in MatSolve_SeqSBAIJ_7_inplace()
290 x2 = tp[2]; in MatSolve_SeqSBAIJ_7_inplace()
291 x3 = tp[3]; in MatSolve_SeqSBAIJ_7_inplace()
292 x4 = tp[4]; in MatSolve_SeqSBAIJ_7_inplace()
293 x5 = tp[5]; in MatSolve_SeqSBAIJ_7_inplace()
294 x6 = tp[6]; /* xk */ in MatSolve_SeqSBAIJ_7_inplace()
297 tp = t + (*vj) * 7; in MatSolve_SeqSBAIJ_7_inplace()
300 …x0 += v[0] * tp[0] + v[7] * tp[1] + v[14] * tp[2] + v[21] * tp[3] + v[28] * tp[4] + v[35] * tp[5] … in MatSolve_SeqSBAIJ_7_inplace()
301 …x1 += v[1] * tp[0] + v[8] * tp[1] + v[15] * tp[2] + v[22] * tp[3] + v[29] * tp[4] + v[36] * tp[5] … in MatSolve_SeqSBAIJ_7_inplace()
302 …x2 += v[2] * tp[0] + v[9] * tp[1] + v[16] * tp[2] + v[23] * tp[3] + v[30] * tp[4] + v[37] * tp[5] … in MatSolve_SeqSBAIJ_7_inplace()
303 …x3 += v[3] * tp[0] + v[10] * tp[1] + v[17] * tp[2] + v[24] * tp[3] + v[31] * tp[4] + v[38] * tp[5]… in MatSolve_SeqSBAIJ_7_inplace()
304 …x4 += v[4] * tp[0] + v[11] * tp[1] + v[18] * tp[2] + v[25] * tp[3] + v[32] * tp[4] + v[39] * tp[5]… in MatSolve_SeqSBAIJ_7_inplace()
305 …x5 += v[5] * tp[0] + v[12] * tp[1] + v[19] * tp[2] + v[26] * tp[3] + v[33] * tp[4] + v[40] * tp[5]… in MatSolve_SeqSBAIJ_7_inplace()
306 …x6 += v[6] * tp[0] + v[13] * tp[1] + v[20] * tp[2] + v[27] * tp[3] + v[34] * tp[4] + v[41] * tp[5]… in MatSolve_SeqSBAIJ_7_inplace()
308 tp = t + (*vj) * 7; in MatSolve_SeqSBAIJ_7_inplace()
311 tp = t + k * 7; in MatSolve_SeqSBAIJ_7_inplace()
312 tp[0] = x0; in MatSolve_SeqSBAIJ_7_inplace()
313 tp[1] = x1; in MatSolve_SeqSBAIJ_7_inplace()
314 tp[2] = x2; in MatSolve_SeqSBAIJ_7_inplace()
315 tp[3] = x3; in MatSolve_SeqSBAIJ_7_inplace()
316 tp[4] = x4; in MatSolve_SeqSBAIJ_7_inplace()
317 tp[5] = x5; in MatSolve_SeqSBAIJ_7_inplace()
318 tp[6] = x6; in MatSolve_SeqSBAIJ_7_inplace()
502 PetscScalar *x, x0, x1, x2, x3, x4, x5, *t, *tp; in MatSolve_SeqSBAIJ_6_inplace() local
512 tp = t; in MatSolve_SeqSBAIJ_6_inplace()
515 tp[0] = b[idx]; in MatSolve_SeqSBAIJ_6_inplace()
516 tp[1] = b[idx + 1]; in MatSolve_SeqSBAIJ_6_inplace()
517 tp[2] = b[idx + 2]; in MatSolve_SeqSBAIJ_6_inplace()
518 tp[3] = b[idx + 3]; in MatSolve_SeqSBAIJ_6_inplace()
519 tp[4] = b[idx + 4]; in MatSolve_SeqSBAIJ_6_inplace()
520 tp[5] = b[idx + 5]; in MatSolve_SeqSBAIJ_6_inplace()
521 tp += 6; in MatSolve_SeqSBAIJ_6_inplace()
527 tp = t + k * 6; in MatSolve_SeqSBAIJ_6_inplace()
528 x0 = tp[0]; in MatSolve_SeqSBAIJ_6_inplace()
529 x1 = tp[1]; in MatSolve_SeqSBAIJ_6_inplace()
530 x2 = tp[2]; in MatSolve_SeqSBAIJ_6_inplace()
531 x3 = tp[3]; in MatSolve_SeqSBAIJ_6_inplace()
532 x4 = tp[4]; in MatSolve_SeqSBAIJ_6_inplace()
533 x5 = tp[5]; in MatSolve_SeqSBAIJ_6_inplace()
535 tp = t + (*vj) * 6; in MatSolve_SeqSBAIJ_6_inplace()
537 tp[0] += v[0] * x0 + v[1] * x1 + v[2] * x2 + v[3] * x3 + v[4] * x4 + v[5] * x5; in MatSolve_SeqSBAIJ_6_inplace()
538 tp[1] += v[6] * x0 + v[7] * x1 + v[8] * x2 + v[9] * x3 + v[10] * x4 + v[11] * x5; in MatSolve_SeqSBAIJ_6_inplace()
539 tp[2] += v[12] * x0 + v[13] * x1 + v[14] * x2 + v[15] * x3 + v[16] * x4 + v[17] * x5; in MatSolve_SeqSBAIJ_6_inplace()
540 tp[3] += v[18] * x0 + v[19] * x1 + v[20] * x2 + v[21] * x3 + v[22] * x4 + v[23] * x5; in MatSolve_SeqSBAIJ_6_inplace()
541 tp[4] += v[24] * x0 + v[25] * x1 + v[26] * x2 + v[27] * x3 + v[28] * x4 + v[29] * x5; in MatSolve_SeqSBAIJ_6_inplace()
542 tp[5] += v[30] * x0 + v[31] * x1 + v[32] * x2 + v[33] * x3 + v[34] * x4 + v[35] * x5; in MatSolve_SeqSBAIJ_6_inplace()
544 tp = t + (*vj) * 6; in MatSolve_SeqSBAIJ_6_inplace()
550 tp = t + k * 6; in MatSolve_SeqSBAIJ_6_inplace()
551 tp[0] = d[0] * x0 + d[6] * x1 + d[12] * x2 + d[18] * x3 + d[24] * x4 + d[30] * x5; in MatSolve_SeqSBAIJ_6_inplace()
552 tp[1] = d[1] * x0 + d[7] * x1 + d[13] * x2 + d[19] * x3 + d[25] * x4 + d[31] * x5; in MatSolve_SeqSBAIJ_6_inplace()
553 tp[2] = d[2] * x0 + d[8] * x1 + d[14] * x2 + d[20] * x3 + d[26] * x4 + d[32] * x5; in MatSolve_SeqSBAIJ_6_inplace()
554 tp[3] = d[3] * x0 + d[9] * x1 + d[15] * x2 + d[21] * x3 + d[27] * x4 + d[33] * x5; in MatSolve_SeqSBAIJ_6_inplace()
555 tp[4] = d[4] * x0 + d[10] * x1 + d[16] * x2 + d[22] * x3 + d[28] * x4 + d[34] * x5; in MatSolve_SeqSBAIJ_6_inplace()
556 tp[5] = d[5] * x0 + d[11] * x1 + d[17] * x2 + d[23] * x3 + d[29] * x4 + d[35] * x5; in MatSolve_SeqSBAIJ_6_inplace()
563 tp = t + k * 6; in MatSolve_SeqSBAIJ_6_inplace()
564 x0 = tp[0]; in MatSolve_SeqSBAIJ_6_inplace()
565 x1 = tp[1]; in MatSolve_SeqSBAIJ_6_inplace()
566 x2 = tp[2]; in MatSolve_SeqSBAIJ_6_inplace()
567 x3 = tp[3]; in MatSolve_SeqSBAIJ_6_inplace()
568 x4 = tp[4]; in MatSolve_SeqSBAIJ_6_inplace()
569 x5 = tp[5]; /* xk */ in MatSolve_SeqSBAIJ_6_inplace()
572 tp = t + (*vj) * 6; in MatSolve_SeqSBAIJ_6_inplace()
575 … x0 += v[0] * tp[0] + v[6] * tp[1] + v[12] * tp[2] + v[18] * tp[3] + v[24] * tp[4] + v[30] * tp[5]; in MatSolve_SeqSBAIJ_6_inplace()
576 … x1 += v[1] * tp[0] + v[7] * tp[1] + v[13] * tp[2] + v[19] * tp[3] + v[25] * tp[4] + v[31] * tp[5]; in MatSolve_SeqSBAIJ_6_inplace()
577 … x2 += v[2] * tp[0] + v[8] * tp[1] + v[14] * tp[2] + v[20] * tp[3] + v[26] * tp[4] + v[32] * tp[5]; in MatSolve_SeqSBAIJ_6_inplace()
578 … x3 += v[3] * tp[0] + v[9] * tp[1] + v[15] * tp[2] + v[21] * tp[3] + v[27] * tp[4] + v[33] * tp[5]; in MatSolve_SeqSBAIJ_6_inplace()
579 …x4 += v[4] * tp[0] + v[10] * tp[1] + v[16] * tp[2] + v[22] * tp[3] + v[28] * tp[4] + v[34] * tp[5]; in MatSolve_SeqSBAIJ_6_inplace()
580 …x5 += v[5] * tp[0] + v[11] * tp[1] + v[17] * tp[2] + v[23] * tp[3] + v[29] * tp[4] + v[35] * tp[5]; in MatSolve_SeqSBAIJ_6_inplace()
582 tp = t + (*vj) * 6; in MatSolve_SeqSBAIJ_6_inplace()
585 tp = t + k * 6; in MatSolve_SeqSBAIJ_6_inplace()
586 tp[0] = x0; in MatSolve_SeqSBAIJ_6_inplace()
587 tp[1] = x1; in MatSolve_SeqSBAIJ_6_inplace()
588 tp[2] = x2; in MatSolve_SeqSBAIJ_6_inplace()
589 tp[3] = x3; in MatSolve_SeqSBAIJ_6_inplace()
590 tp[4] = x4; in MatSolve_SeqSBAIJ_6_inplace()
591 tp[5] = x5; in MatSolve_SeqSBAIJ_6_inplace()
768 PetscScalar *x, x0, x1, x2, x3, x4, *t, *tp; in MatSolve_SeqSBAIJ_5_inplace() local
778 tp = t; in MatSolve_SeqSBAIJ_5_inplace()
781 tp[0] = b[idx]; in MatSolve_SeqSBAIJ_5_inplace()
782 tp[1] = b[idx + 1]; in MatSolve_SeqSBAIJ_5_inplace()
783 tp[2] = b[idx + 2]; in MatSolve_SeqSBAIJ_5_inplace()
784 tp[3] = b[idx + 3]; in MatSolve_SeqSBAIJ_5_inplace()
785 tp[4] = b[idx + 4]; in MatSolve_SeqSBAIJ_5_inplace()
786 tp += 5; in MatSolve_SeqSBAIJ_5_inplace()
792 tp = t + k * 5; in MatSolve_SeqSBAIJ_5_inplace()
793 x0 = tp[0]; in MatSolve_SeqSBAIJ_5_inplace()
794 x1 = tp[1]; in MatSolve_SeqSBAIJ_5_inplace()
795 x2 = tp[2]; in MatSolve_SeqSBAIJ_5_inplace()
796 x3 = tp[3]; in MatSolve_SeqSBAIJ_5_inplace()
797 x4 = tp[4]; in MatSolve_SeqSBAIJ_5_inplace()
800 tp = t + (*vj) * 5; in MatSolve_SeqSBAIJ_5_inplace()
802 tp[0] += v[0] * x0 + v[1] * x1 + v[2] * x2 + v[3] * x3 + v[4] * x4; in MatSolve_SeqSBAIJ_5_inplace()
803 tp[1] += v[5] * x0 + v[6] * x1 + v[7] * x2 + v[8] * x3 + v[9] * x4; in MatSolve_SeqSBAIJ_5_inplace()
804 tp[2] += v[10] * x0 + v[11] * x1 + v[12] * x2 + v[13] * x3 + v[14] * x4; in MatSolve_SeqSBAIJ_5_inplace()
805 tp[3] += v[15] * x0 + v[16] * x1 + v[17] * x2 + v[18] * x3 + v[19] * x4; in MatSolve_SeqSBAIJ_5_inplace()
806 tp[4] += v[20] * x0 + v[21] * x1 + v[22] * x2 + v[23] * x3 + v[24] * x4; in MatSolve_SeqSBAIJ_5_inplace()
808 tp = t + (*vj) * 5; in MatSolve_SeqSBAIJ_5_inplace()
814 tp = t + k * 5; in MatSolve_SeqSBAIJ_5_inplace()
815 tp[0] = diag[0] * x0 + diag[5] * x1 + diag[10] * x2 + diag[15] * x3 + diag[20] * x4; in MatSolve_SeqSBAIJ_5_inplace()
816 tp[1] = diag[1] * x0 + diag[6] * x1 + diag[11] * x2 + diag[16] * x3 + diag[21] * x4; in MatSolve_SeqSBAIJ_5_inplace()
817 tp[2] = diag[2] * x0 + diag[7] * x1 + diag[12] * x2 + diag[17] * x3 + diag[22] * x4; in MatSolve_SeqSBAIJ_5_inplace()
818 tp[3] = diag[3] * x0 + diag[8] * x1 + diag[13] * x2 + diag[18] * x3 + diag[23] * x4; in MatSolve_SeqSBAIJ_5_inplace()
819 tp[4] = diag[4] * x0 + diag[9] * x1 + diag[14] * x2 + diag[19] * x3 + diag[24] * x4; in MatSolve_SeqSBAIJ_5_inplace()
826 tp = t + k * 5; in MatSolve_SeqSBAIJ_5_inplace()
827 x0 = tp[0]; in MatSolve_SeqSBAIJ_5_inplace()
828 x1 = tp[1]; in MatSolve_SeqSBAIJ_5_inplace()
829 x2 = tp[2]; in MatSolve_SeqSBAIJ_5_inplace()
830 x3 = tp[3]; in MatSolve_SeqSBAIJ_5_inplace()
831 x4 = tp[4]; /* xk */ in MatSolve_SeqSBAIJ_5_inplace()
834 tp = t + (*vj) * 5; in MatSolve_SeqSBAIJ_5_inplace()
837 x0 += v[0] * tp[0] + v[5] * tp[1] + v[10] * tp[2] + v[15] * tp[3] + v[20] * tp[4]; in MatSolve_SeqSBAIJ_5_inplace()
838 x1 += v[1] * tp[0] + v[6] * tp[1] + v[11] * tp[2] + v[16] * tp[3] + v[21] * tp[4]; in MatSolve_SeqSBAIJ_5_inplace()
839 x2 += v[2] * tp[0] + v[7] * tp[1] + v[12] * tp[2] + v[17] * tp[3] + v[22] * tp[4]; in MatSolve_SeqSBAIJ_5_inplace()
840 x3 += v[3] * tp[0] + v[8] * tp[1] + v[13] * tp[2] + v[18] * tp[3] + v[23] * tp[4]; in MatSolve_SeqSBAIJ_5_inplace()
841 x4 += v[4] * tp[0] + v[9] * tp[1] + v[14] * tp[2] + v[19] * tp[3] + v[24] * tp[4]; in MatSolve_SeqSBAIJ_5_inplace()
843 tp = t + (*vj) * 5; in MatSolve_SeqSBAIJ_5_inplace()
846 tp = t + k * 5; in MatSolve_SeqSBAIJ_5_inplace()
847 tp[0] = x0; in MatSolve_SeqSBAIJ_5_inplace()
848 tp[1] = x1; in MatSolve_SeqSBAIJ_5_inplace()
849 tp[2] = x2; in MatSolve_SeqSBAIJ_5_inplace()
850 tp[3] = x3; in MatSolve_SeqSBAIJ_5_inplace()
851 tp[4] = x4; in MatSolve_SeqSBAIJ_5_inplace()
1022 PetscScalar *x, x0, x1, x2, x3, *t, *tp; in MatSolve_SeqSBAIJ_4_inplace() local
1032 tp = t; in MatSolve_SeqSBAIJ_4_inplace()
1035 tp[0] = b[idx]; in MatSolve_SeqSBAIJ_4_inplace()
1036 tp[1] = b[idx + 1]; in MatSolve_SeqSBAIJ_4_inplace()
1037 tp[2] = b[idx + 2]; in MatSolve_SeqSBAIJ_4_inplace()
1038 tp[3] = b[idx + 3]; in MatSolve_SeqSBAIJ_4_inplace()
1039 tp += 4; in MatSolve_SeqSBAIJ_4_inplace()
1045 tp = t + k * 4; in MatSolve_SeqSBAIJ_4_inplace()
1046 x0 = tp[0]; in MatSolve_SeqSBAIJ_4_inplace()
1047 x1 = tp[1]; in MatSolve_SeqSBAIJ_4_inplace()
1048 x2 = tp[2]; in MatSolve_SeqSBAIJ_4_inplace()
1049 x3 = tp[3]; in MatSolve_SeqSBAIJ_4_inplace()
1052 tp = t + (*vj) * 4; in MatSolve_SeqSBAIJ_4_inplace()
1054 tp[0] += v[0] * x0 + v[1] * x1 + v[2] * x2 + v[3] * x3; in MatSolve_SeqSBAIJ_4_inplace()
1055 tp[1] += v[4] * x0 + v[5] * x1 + v[6] * x2 + v[7] * x3; in MatSolve_SeqSBAIJ_4_inplace()
1056 tp[2] += v[8] * x0 + v[9] * x1 + v[10] * x2 + v[11] * x3; in MatSolve_SeqSBAIJ_4_inplace()
1057 tp[3] += v[12] * x0 + v[13] * x1 + v[14] * x2 + v[15] * x3; in MatSolve_SeqSBAIJ_4_inplace()
1059 tp = t + (*vj) * 4; in MatSolve_SeqSBAIJ_4_inplace()
1065 tp = t + k * 4; in MatSolve_SeqSBAIJ_4_inplace()
1066 tp[0] = diag[0] * x0 + diag[4] * x1 + diag[8] * x2 + diag[12] * x3; in MatSolve_SeqSBAIJ_4_inplace()
1067 tp[1] = diag[1] * x0 + diag[5] * x1 + diag[9] * x2 + diag[13] * x3; in MatSolve_SeqSBAIJ_4_inplace()
1068 tp[2] = diag[2] * x0 + diag[6] * x1 + diag[10] * x2 + diag[14] * x3; in MatSolve_SeqSBAIJ_4_inplace()
1069 tp[3] = diag[3] * x0 + diag[7] * x1 + diag[11] * x2 + diag[15] * x3; in MatSolve_SeqSBAIJ_4_inplace()
1076 tp = t + k * 4; in MatSolve_SeqSBAIJ_4_inplace()
1077 x0 = tp[0]; in MatSolve_SeqSBAIJ_4_inplace()
1078 x1 = tp[1]; in MatSolve_SeqSBAIJ_4_inplace()
1079 x2 = tp[2]; in MatSolve_SeqSBAIJ_4_inplace()
1080 x3 = tp[3]; /* xk */ in MatSolve_SeqSBAIJ_4_inplace()
1083 tp = t + (*vj) * 4; in MatSolve_SeqSBAIJ_4_inplace()
1086 x0 += v[0] * tp[0] + v[4] * tp[1] + v[8] * tp[2] + v[12] * tp[3]; in MatSolve_SeqSBAIJ_4_inplace()
1087 x1 += v[1] * tp[0] + v[5] * tp[1] + v[9] * tp[2] + v[13] * tp[3]; in MatSolve_SeqSBAIJ_4_inplace()
1088 x2 += v[2] * tp[0] + v[6] * tp[1] + v[10] * tp[2] + v[14] * tp[3]; in MatSolve_SeqSBAIJ_4_inplace()
1089 x3 += v[3] * tp[0] + v[7] * tp[1] + v[11] * tp[2] + v[15] * tp[3]; in MatSolve_SeqSBAIJ_4_inplace()
1091 tp = t + (*vj) * 4; in MatSolve_SeqSBAIJ_4_inplace()
1094 tp = t + k * 4; in MatSolve_SeqSBAIJ_4_inplace()
1095 tp[0] = x0; in MatSolve_SeqSBAIJ_4_inplace()
1096 tp[1] = x1; in MatSolve_SeqSBAIJ_4_inplace()
1097 tp[2] = x2; in MatSolve_SeqSBAIJ_4_inplace()
1098 tp[3] = x3; in MatSolve_SeqSBAIJ_4_inplace()
1262 PetscScalar *x, x0, x1, x2, *t, *tp; in MatSolve_SeqSBAIJ_3_inplace() local
1272 tp = t; in MatSolve_SeqSBAIJ_3_inplace()
1275 tp[0] = b[idx]; in MatSolve_SeqSBAIJ_3_inplace()
1276 tp[1] = b[idx + 1]; in MatSolve_SeqSBAIJ_3_inplace()
1277 tp[2] = b[idx + 2]; in MatSolve_SeqSBAIJ_3_inplace()
1278 tp += 3; in MatSolve_SeqSBAIJ_3_inplace()
1284 tp = t + k * 3; in MatSolve_SeqSBAIJ_3_inplace()
1285 x0 = tp[0]; in MatSolve_SeqSBAIJ_3_inplace()
1286 x1 = tp[1]; in MatSolve_SeqSBAIJ_3_inplace()
1287 x2 = tp[2]; in MatSolve_SeqSBAIJ_3_inplace()
1290 tp = t + (*vj) * 3; in MatSolve_SeqSBAIJ_3_inplace()
1292 tp[0] += v[0] * x0 + v[1] * x1 + v[2] * x2; in MatSolve_SeqSBAIJ_3_inplace()
1293 tp[1] += v[3] * x0 + v[4] * x1 + v[5] * x2; in MatSolve_SeqSBAIJ_3_inplace()
1294 tp[2] += v[6] * x0 + v[7] * x1 + v[8] * x2; in MatSolve_SeqSBAIJ_3_inplace()
1296 tp = t + (*vj) * 3; in MatSolve_SeqSBAIJ_3_inplace()
1302 tp = t + k * 3; in MatSolve_SeqSBAIJ_3_inplace()
1303 tp[0] = diag[0] * x0 + diag[3] * x1 + diag[6] * x2; in MatSolve_SeqSBAIJ_3_inplace()
1304 tp[1] = diag[1] * x0 + diag[4] * x1 + diag[7] * x2; in MatSolve_SeqSBAIJ_3_inplace()
1305 tp[2] = diag[2] * x0 + diag[5] * x1 + diag[8] * x2; in MatSolve_SeqSBAIJ_3_inplace()
1312 tp = t + k * 3; in MatSolve_SeqSBAIJ_3_inplace()
1313 x0 = tp[0]; in MatSolve_SeqSBAIJ_3_inplace()
1314 x1 = tp[1]; in MatSolve_SeqSBAIJ_3_inplace()
1315 x2 = tp[2]; /* xk */ in MatSolve_SeqSBAIJ_3_inplace()
1318 tp = t + (*vj) * 3; in MatSolve_SeqSBAIJ_3_inplace()
1321 x0 += v[0] * tp[0] + v[3] * tp[1] + v[6] * tp[2]; in MatSolve_SeqSBAIJ_3_inplace()
1322 x1 += v[1] * tp[0] + v[4] * tp[1] + v[7] * tp[2]; in MatSolve_SeqSBAIJ_3_inplace()
1323 x2 += v[2] * tp[0] + v[5] * tp[1] + v[8] * tp[2]; in MatSolve_SeqSBAIJ_3_inplace()
1325 tp = t + (*vj) * 3; in MatSolve_SeqSBAIJ_3_inplace()
1328 tp = t + k * 3; in MatSolve_SeqSBAIJ_3_inplace()
1329 tp[0] = x0; in MatSolve_SeqSBAIJ_3_inplace()
1330 tp[1] = x1; in MatSolve_SeqSBAIJ_3_inplace()
1331 tp[2] = x2; in MatSolve_SeqSBAIJ_3_inplace()