Lines Matching refs:v
292 const MatScalar *v; in MatMult_SeqSBAIJ_2() local
303 v = a->a; in MatMult_SeqSBAIJ_2()
314 z[2 * i] += v[0] * x1 + v[2] * x2; in MatMult_SeqSBAIJ_2()
315 z[2 * i + 1] += v[2] * x1 + v[3] * x2; in MatMult_SeqSBAIJ_2()
316 v += 4; in MatMult_SeqSBAIJ_2()
320 PetscPrefetchBlock(v + 4 * n, 4 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMult_SeqSBAIJ_2()
324 z[cval] += v[0] * x1 + v[1] * x2; in MatMult_SeqSBAIJ_2()
325 z[cval + 1] += v[2] * x1 + v[3] * x2; in MatMult_SeqSBAIJ_2()
327 z[2 * i] += v[0] * x[cval] + v[2] * x[cval + 1]; in MatMult_SeqSBAIJ_2()
328 z[2 * i + 1] += v[1] * x[cval] + v[3] * x[cval + 1]; in MatMult_SeqSBAIJ_2()
329 v += 4; in MatMult_SeqSBAIJ_2()
346 const MatScalar *v; in MatMult_SeqSBAIJ_3() local
357 v = a->a; in MatMult_SeqSBAIJ_3()
369 z[3 * i] += v[0] * x1 + v[3] * x2 + v[6] * x3; in MatMult_SeqSBAIJ_3()
370 z[3 * i + 1] += v[3] * x1 + v[4] * x2 + v[7] * x3; in MatMult_SeqSBAIJ_3()
371 z[3 * i + 2] += v[6] * x1 + v[7] * x2 + v[8] * x3; in MatMult_SeqSBAIJ_3()
372 v += 9; in MatMult_SeqSBAIJ_3()
376 PetscPrefetchBlock(v + 9 * n, 9 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMult_SeqSBAIJ_3()
380 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3; in MatMult_SeqSBAIJ_3()
381 z[cval + 1] += v[3] * x1 + v[4] * x2 + v[5] * x3; in MatMult_SeqSBAIJ_3()
382 z[cval + 2] += v[6] * x1 + v[7] * x2 + v[8] * x3; in MatMult_SeqSBAIJ_3()
384 z[3 * i] += v[0] * x[cval] + v[3] * x[cval + 1] + v[6] * x[cval + 2]; in MatMult_SeqSBAIJ_3()
385 z[3 * i + 1] += v[1] * x[cval] + v[4] * x[cval + 1] + v[7] * x[cval + 2]; in MatMult_SeqSBAIJ_3()
386 z[3 * i + 2] += v[2] * x[cval] + v[5] * x[cval + 1] + v[8] * x[cval + 2]; in MatMult_SeqSBAIJ_3()
387 v += 9; in MatMult_SeqSBAIJ_3()
404 const MatScalar *v; in MatMult_SeqSBAIJ_4() local
415 v = a->a; in MatMult_SeqSBAIJ_4()
428 z[4 * i] += v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatMult_SeqSBAIJ_4()
429 z[4 * i + 1] += v[4] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatMult_SeqSBAIJ_4()
430 z[4 * i + 2] += v[8] * x1 + v[9] * x2 + v[10] * x3 + v[14] * x4; in MatMult_SeqSBAIJ_4()
431 z[4 * i + 3] += v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4; in MatMult_SeqSBAIJ_4()
432 v += 16; in MatMult_SeqSBAIJ_4()
436 … PetscPrefetchBlock(v + 16 * n, 16 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMult_SeqSBAIJ_4()
440 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4; in MatMult_SeqSBAIJ_4()
441 z[cval + 1] += v[4] * x1 + v[5] * x2 + v[6] * x3 + v[7] * x4; in MatMult_SeqSBAIJ_4()
442 z[cval + 2] += v[8] * x1 + v[9] * x2 + v[10] * x3 + v[11] * x4; in MatMult_SeqSBAIJ_4()
443 z[cval + 3] += v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4; in MatMult_SeqSBAIJ_4()
445 z[4 * i] += v[0] * x[cval] + v[4] * x[cval + 1] + v[8] * x[cval + 2] + v[12] * x[cval + 3]; in MatMult_SeqSBAIJ_4()
446 … z[4 * i + 1] += v[1] * x[cval] + v[5] * x[cval + 1] + v[9] * x[cval + 2] + v[13] * x[cval + 3]; in MatMult_SeqSBAIJ_4()
447 … z[4 * i + 2] += v[2] * x[cval] + v[6] * x[cval + 1] + v[10] * x[cval + 2] + v[14] * x[cval + 3]; in MatMult_SeqSBAIJ_4()
448 … z[4 * i + 3] += v[3] * x[cval] + v[7] * x[cval + 1] + v[11] * x[cval + 2] + v[15] * x[cval + 3]; in MatMult_SeqSBAIJ_4()
449 v += 16; in MatMult_SeqSBAIJ_4()
466 const MatScalar *v; in MatMult_SeqSBAIJ_5() local
477 v = a->a; in MatMult_SeqSBAIJ_5()
491 z[5 * i] += v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatMult_SeqSBAIJ_5()
492 z[5 * i + 1] += v[5] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatMult_SeqSBAIJ_5()
493 z[5 * i + 2] += v[10] * x1 + v[11] * x2 + v[12] * x3 + v[17] * x4 + v[22] * x5; in MatMult_SeqSBAIJ_5()
494 z[5 * i + 3] += v[15] * x1 + v[16] * x2 + v[17] * x3 + v[18] * x4 + v[23] * x5; in MatMult_SeqSBAIJ_5()
495 z[5 * i + 4] += v[20] * x1 + v[21] * x2 + v[22] * x3 + v[23] * x4 + v[24] * x5; in MatMult_SeqSBAIJ_5()
496 v += 25; in MatMult_SeqSBAIJ_5()
500 … PetscPrefetchBlock(v + 25 * n, 25 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMult_SeqSBAIJ_5()
504 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5; in MatMult_SeqSBAIJ_5()
505 z[cval + 1] += v[5] * x1 + v[6] * x2 + v[7] * x3 + v[8] * x4 + v[9] * x5; in MatMult_SeqSBAIJ_5()
506 z[cval + 2] += v[10] * x1 + v[11] * x2 + v[12] * x3 + v[13] * x4 + v[14] * x5; in MatMult_SeqSBAIJ_5()
507 z[cval + 3] += v[15] * x1 + v[16] * x2 + v[17] * x3 + v[18] * x4 + v[19] * x5; in MatMult_SeqSBAIJ_5()
508 z[cval + 4] += v[20] * x1 + v[21] * x2 + v[22] * x3 + v[23] * x4 + v[24] * x5; in MatMult_SeqSBAIJ_5()
510 …z[5 * i] += v[0] * x[cval] + v[5] * x[cval + 1] + v[10] * x[cval + 2] + v[15] * x[cval + 3] + v[20… in MatMult_SeqSBAIJ_5()
511 …z[5 * i + 1] += v[1] * x[cval] + v[6] * x[cval + 1] + v[11] * x[cval + 2] + v[16] * x[cval + 3] + … in MatMult_SeqSBAIJ_5()
512 …z[5 * i + 2] += v[2] * x[cval] + v[7] * x[cval + 1] + v[12] * x[cval + 2] + v[17] * x[cval + 3] + … in MatMult_SeqSBAIJ_5()
513 …z[5 * i + 3] += v[3] * x[cval] + v[8] * x[cval + 1] + v[13] * x[cval + 2] + v[18] * x[cval + 3] + … in MatMult_SeqSBAIJ_5()
514 …z[5 * i + 4] += v[4] * x[cval] + v[9] * x[cval + 1] + v[14] * x[cval + 2] + v[19] * x[cval + 3] + … in MatMult_SeqSBAIJ_5()
515 v += 25; in MatMult_SeqSBAIJ_5()
532 const MatScalar *v; in MatMult_SeqSBAIJ_6() local
543 v = a->a; in MatMult_SeqSBAIJ_6()
558 z[6 * i] += v[0] * x1 + v[6] * x2 + v[12] * x3 + v[18] * x4 + v[24] * x5 + v[30] * x6; in MatMult_SeqSBAIJ_6()
559 z[6 * i + 1] += v[6] * x1 + v[7] * x2 + v[13] * x3 + v[19] * x4 + v[25] * x5 + v[31] * x6; in MatMult_SeqSBAIJ_6()
560 z[6 * i + 2] += v[12] * x1 + v[13] * x2 + v[14] * x3 + v[20] * x4 + v[26] * x5 + v[32] * x6; in MatMult_SeqSBAIJ_6()
561 z[6 * i + 3] += v[18] * x1 + v[19] * x2 + v[20] * x3 + v[21] * x4 + v[27] * x5 + v[33] * x6; in MatMult_SeqSBAIJ_6()
562 z[6 * i + 4] += v[24] * x1 + v[25] * x2 + v[26] * x3 + v[27] * x4 + v[28] * x5 + v[34] * x6; in MatMult_SeqSBAIJ_6()
563 z[6 * i + 5] += v[30] * x1 + v[31] * x2 + v[32] * x3 + v[33] * x4 + v[34] * x5 + v[35] * x6; in MatMult_SeqSBAIJ_6()
564 v += 36; in MatMult_SeqSBAIJ_6()
568 … PetscPrefetchBlock(v + 36 * n, 36 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMult_SeqSBAIJ_6()
572 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5 + v[5] * x6; in MatMult_SeqSBAIJ_6()
573 z[cval + 1] += v[6] * x1 + v[7] * x2 + v[8] * x3 + v[9] * x4 + v[10] * x5 + v[11] * x6; in MatMult_SeqSBAIJ_6()
574 z[cval + 2] += v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4 + v[16] * x5 + v[17] * x6; in MatMult_SeqSBAIJ_6()
575 z[cval + 3] += v[18] * x1 + v[19] * x2 + v[20] * x3 + v[21] * x4 + v[22] * x5 + v[23] * x6; in MatMult_SeqSBAIJ_6()
576 z[cval + 4] += v[24] * x1 + v[25] * x2 + v[26] * x3 + v[27] * x4 + v[28] * x5 + v[29] * x6; in MatMult_SeqSBAIJ_6()
577 z[cval + 5] += v[30] * x1 + v[31] * x2 + v[32] * x3 + v[33] * x4 + v[34] * x5 + v[35] * x6; in MatMult_SeqSBAIJ_6()
579 …z[6 * i] += v[0] * x[cval] + v[6] * x[cval + 1] + v[12] * x[cval + 2] + v[18] * x[cval + 3] + v[24… in MatMult_SeqSBAIJ_6()
580 …z[6 * i + 1] += v[1] * x[cval] + v[7] * x[cval + 1] + v[13] * x[cval + 2] + v[19] * x[cval + 3] + … in MatMult_SeqSBAIJ_6()
581 …z[6 * i + 2] += v[2] * x[cval] + v[8] * x[cval + 1] + v[14] * x[cval + 2] + v[20] * x[cval + 3] + … in MatMult_SeqSBAIJ_6()
582 …z[6 * i + 3] += v[3] * x[cval] + v[9] * x[cval + 1] + v[15] * x[cval + 2] + v[21] * x[cval + 3] + … in MatMult_SeqSBAIJ_6()
583 …z[6 * i + 4] += v[4] * x[cval] + v[10] * x[cval + 1] + v[16] * x[cval + 2] + v[22] * x[cval + 3] +… in MatMult_SeqSBAIJ_6()
584 …z[6 * i + 5] += v[5] * x[cval] + v[11] * x[cval + 1] + v[17] * x[cval + 2] + v[23] * x[cval + 3] +… in MatMult_SeqSBAIJ_6()
585 v += 36; in MatMult_SeqSBAIJ_6()
602 const MatScalar *v; in MatMult_SeqSBAIJ_7() local
613 v = a->a; in MatMult_SeqSBAIJ_7()
629 …z[7 * i] += v[0] * x1 + v[7] * x2 + v[14] * x3 + v[21] * x4 + v[28] * x5 + v[35] * x6 + v[42] * x7; in MatMult_SeqSBAIJ_7()
630 …z[7 * i + 1] += v[7] * x1 + v[8] * x2 + v[15] * x3 + v[22] * x4 + v[29] * x5 + v[36] * x6 + v[43] … in MatMult_SeqSBAIJ_7()
631 …z[7 * i + 2] += v[14] * x1 + v[15] * x2 + v[16] * x3 + v[23] * x4 + v[30] * x5 + v[37] * x6 + v[44… in MatMult_SeqSBAIJ_7()
632 …z[7 * i + 3] += v[21] * x1 + v[22] * x2 + v[23] * x3 + v[24] * x4 + v[31] * x5 + v[38] * x6 + v[45… in MatMult_SeqSBAIJ_7()
633 …z[7 * i + 4] += v[28] * x1 + v[29] * x2 + v[30] * x3 + v[31] * x4 + v[32] * x5 + v[39] * x6 + v[46… in MatMult_SeqSBAIJ_7()
634 …z[7 * i + 5] += v[35] * x1 + v[36] * x2 + v[37] * x3 + v[38] * x4 + v[39] * x5 + v[40] * x6 + v[47… in MatMult_SeqSBAIJ_7()
635 …z[7 * i + 6] += v[42] * x1 + v[43] * x2 + v[44] * x3 + v[45] * x4 + v[46] * x5 + v[47] * x6 + v[48… in MatMult_SeqSBAIJ_7()
636 v += 49; in MatMult_SeqSBAIJ_7()
640 … PetscPrefetchBlock(v + 49 * n, 49 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMult_SeqSBAIJ_7()
644 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5 + v[5] * x6 + v[6] * x7; in MatMult_SeqSBAIJ_7()
645 …z[cval + 1] += v[7] * x1 + v[8] * x2 + v[9] * x3 + v[10] * x4 + v[11] * x5 + v[12] * x6 + v[13] * … in MatMult_SeqSBAIJ_7()
646 …z[cval + 2] += v[14] * x1 + v[15] * x2 + v[16] * x3 + v[17] * x4 + v[18] * x5 + v[19] * x6 + v[20]… in MatMult_SeqSBAIJ_7()
647 …z[cval + 3] += v[21] * x1 + v[22] * x2 + v[23] * x3 + v[24] * x4 + v[25] * x5 + v[26] * x6 + v[27]… in MatMult_SeqSBAIJ_7()
648 …z[cval + 4] += v[28] * x1 + v[29] * x2 + v[30] * x3 + v[31] * x4 + v[32] * x5 + v[33] * x6 + v[34]… in MatMult_SeqSBAIJ_7()
649 …z[cval + 5] += v[35] * x1 + v[36] * x2 + v[37] * x3 + v[38] * x4 + v[39] * x5 + v[40] * x6 + v[41]… in MatMult_SeqSBAIJ_7()
650 …z[cval + 6] += v[42] * x1 + v[43] * x2 + v[44] * x3 + v[45] * x4 + v[46] * x5 + v[47] * x6 + v[48]… in MatMult_SeqSBAIJ_7()
652 …z[7 * i] += v[0] * x[cval] + v[7] * x[cval + 1] + v[14] * x[cval + 2] + v[21] * x[cval + 3] + v[28… in MatMult_SeqSBAIJ_7()
653 … * i + 1] += v[1] * x[cval] + v[8] * x[cval + 1] + v[15] * x[cval + 2] + v[22] * x[cval + 3] + v[2… in MatMult_SeqSBAIJ_7()
654 … * i + 2] += v[2] * x[cval] + v[9] * x[cval + 1] + v[16] * x[cval + 2] + v[23] * x[cval + 3] + v[3… in MatMult_SeqSBAIJ_7()
655 …* i + 3] += v[3] * x[cval] + v[10] * x[cval + 1] + v[17] * x[cval + 2] + v[24] * x[cval + 3] + v[3… in MatMult_SeqSBAIJ_7()
656 …* i + 4] += v[4] * x[cval] + v[11] * x[cval + 1] + v[18] * x[cval + 2] + v[25] * x[cval + 3] + v[3… in MatMult_SeqSBAIJ_7()
657 …* i + 5] += v[5] * x[cval] + v[12] * x[cval + 1] + v[19] * x[cval + 2] + v[26] * x[cval + 3] + v[3… in MatMult_SeqSBAIJ_7()
658 …* i + 6] += v[6] * x[cval] + v[13] * x[cval + 1] + v[20] * x[cval + 2] + v[27] * x[cval + 3] + v[3… in MatMult_SeqSBAIJ_7()
659 v += 49; in MatMult_SeqSBAIJ_7()
678 const MatScalar *v; in MatMult_SeqSBAIJ_N() local
693 v = a->a; in MatMult_SeqSBAIJ_N()
713 PetscKernel_w_gets_w_plus_Ar_times_v(bs, ncols, work, v, z); in MatMult_SeqSBAIJ_N()
719 v += bs2; in MatMult_SeqSBAIJ_N()
727 PetscKernel_w_gets_w_plus_trans_Ar_times_v(bs, ncols, x, v, workt); in MatMult_SeqSBAIJ_N()
735 v += n * bs2; in MatMult_SeqSBAIJ_N()
751 const MatScalar *v; in MatMultAdd_SeqSBAIJ_1() local
762 v = a->a; in MatMultAdd_SeqSBAIJ_1()
772 z[i] += *v++ * x[*ib++]; in MatMultAdd_SeqSBAIJ_1()
778 z[cval] += PetscConj(*v) * x1; /* (strict lower triangular part of A)*x */ in MatMultAdd_SeqSBAIJ_1()
779 z[i] += *v++ * x[*ib++]; /* (strict upper triangular part of A)*x */ in MatMultAdd_SeqSBAIJ_1()
784 z[cval] += *v * x1; /* (strict lower triangular part of A)*x */ in MatMultAdd_SeqSBAIJ_1()
785 z[i] += *v++ * x[*ib++]; /* (strict upper triangular part of A)*x */ in MatMultAdd_SeqSBAIJ_1()
804 const MatScalar *v; in MatMultAdd_SeqSBAIJ_2() local
815 v = a->a; in MatMultAdd_SeqSBAIJ_2()
826 z[2 * i] += v[0] * x1 + v[2] * x2; in MatMultAdd_SeqSBAIJ_2()
827 z[2 * i + 1] += v[2] * x1 + v[3] * x2; in MatMultAdd_SeqSBAIJ_2()
828 v += 4; in MatMultAdd_SeqSBAIJ_2()
832 PetscPrefetchBlock(v + 4 * n, 4 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMultAdd_SeqSBAIJ_2()
836 z[cval] += v[0] * x1 + v[1] * x2; in MatMultAdd_SeqSBAIJ_2()
837 z[cval + 1] += v[2] * x1 + v[3] * x2; in MatMultAdd_SeqSBAIJ_2()
839 z[2 * i] += v[0] * x[cval] + v[2] * x[cval + 1]; in MatMultAdd_SeqSBAIJ_2()
840 z[2 * i + 1] += v[1] * x[cval] + v[3] * x[cval + 1]; in MatMultAdd_SeqSBAIJ_2()
841 v += 4; in MatMultAdd_SeqSBAIJ_2()
858 const MatScalar *v; in MatMultAdd_SeqSBAIJ_3() local
869 v = a->a; in MatMultAdd_SeqSBAIJ_3()
881 z[3 * i] += v[0] * x1 + v[3] * x2 + v[6] * x3; in MatMultAdd_SeqSBAIJ_3()
882 z[3 * i + 1] += v[3] * x1 + v[4] * x2 + v[7] * x3; in MatMultAdd_SeqSBAIJ_3()
883 z[3 * i + 2] += v[6] * x1 + v[7] * x2 + v[8] * x3; in MatMultAdd_SeqSBAIJ_3()
884 v += 9; in MatMultAdd_SeqSBAIJ_3()
888 PetscPrefetchBlock(v + 9 * n, 9 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMultAdd_SeqSBAIJ_3()
892 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3; in MatMultAdd_SeqSBAIJ_3()
893 z[cval + 1] += v[3] * x1 + v[4] * x2 + v[5] * x3; in MatMultAdd_SeqSBAIJ_3()
894 z[cval + 2] += v[6] * x1 + v[7] * x2 + v[8] * x3; in MatMultAdd_SeqSBAIJ_3()
896 z[3 * i] += v[0] * x[cval] + v[3] * x[cval + 1] + v[6] * x[cval + 2]; in MatMultAdd_SeqSBAIJ_3()
897 z[3 * i + 1] += v[1] * x[cval] + v[4] * x[cval + 1] + v[7] * x[cval + 2]; in MatMultAdd_SeqSBAIJ_3()
898 z[3 * i + 2] += v[2] * x[cval] + v[5] * x[cval + 1] + v[8] * x[cval + 2]; in MatMultAdd_SeqSBAIJ_3()
899 v += 9; in MatMultAdd_SeqSBAIJ_3()
917 const MatScalar *v; in MatMultAdd_SeqSBAIJ_4() local
928 v = a->a; in MatMultAdd_SeqSBAIJ_4()
941 z[4 * i] += v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatMultAdd_SeqSBAIJ_4()
942 z[4 * i + 1] += v[4] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatMultAdd_SeqSBAIJ_4()
943 z[4 * i + 2] += v[8] * x1 + v[9] * x2 + v[10] * x3 + v[14] * x4; in MatMultAdd_SeqSBAIJ_4()
944 z[4 * i + 3] += v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4; in MatMultAdd_SeqSBAIJ_4()
945 v += 16; in MatMultAdd_SeqSBAIJ_4()
949 … PetscPrefetchBlock(v + 16 * n, 16 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMultAdd_SeqSBAIJ_4()
953 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4; in MatMultAdd_SeqSBAIJ_4()
954 z[cval + 1] += v[4] * x1 + v[5] * x2 + v[6] * x3 + v[7] * x4; in MatMultAdd_SeqSBAIJ_4()
955 z[cval + 2] += v[8] * x1 + v[9] * x2 + v[10] * x3 + v[11] * x4; in MatMultAdd_SeqSBAIJ_4()
956 z[cval + 3] += v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4; in MatMultAdd_SeqSBAIJ_4()
958 z[4 * i] += v[0] * x[cval] + v[4] * x[cval + 1] + v[8] * x[cval + 2] + v[12] * x[cval + 3]; in MatMultAdd_SeqSBAIJ_4()
959 … z[4 * i + 1] += v[1] * x[cval] + v[5] * x[cval + 1] + v[9] * x[cval + 2] + v[13] * x[cval + 3]; in MatMultAdd_SeqSBAIJ_4()
960 … z[4 * i + 2] += v[2] * x[cval] + v[6] * x[cval + 1] + v[10] * x[cval + 2] + v[14] * x[cval + 3]; in MatMultAdd_SeqSBAIJ_4()
961 … z[4 * i + 3] += v[3] * x[cval] + v[7] * x[cval + 1] + v[11] * x[cval + 2] + v[15] * x[cval + 3]; in MatMultAdd_SeqSBAIJ_4()
962 v += 16; in MatMultAdd_SeqSBAIJ_4()
980 const MatScalar *v; in MatMultAdd_SeqSBAIJ_5() local
991 v = a->a; in MatMultAdd_SeqSBAIJ_5()
1005 z[5 * i] += v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatMultAdd_SeqSBAIJ_5()
1006 z[5 * i + 1] += v[5] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatMultAdd_SeqSBAIJ_5()
1007 z[5 * i + 2] += v[10] * x1 + v[11] * x2 + v[12] * x3 + v[17] * x4 + v[22] * x5; in MatMultAdd_SeqSBAIJ_5()
1008 z[5 * i + 3] += v[15] * x1 + v[16] * x2 + v[17] * x3 + v[18] * x4 + v[23] * x5; in MatMultAdd_SeqSBAIJ_5()
1009 z[5 * i + 4] += v[20] * x1 + v[21] * x2 + v[22] * x3 + v[23] * x4 + v[24] * x5; in MatMultAdd_SeqSBAIJ_5()
1010 v += 25; in MatMultAdd_SeqSBAIJ_5()
1014 … PetscPrefetchBlock(v + 25 * n, 25 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMultAdd_SeqSBAIJ_5()
1018 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5; in MatMultAdd_SeqSBAIJ_5()
1019 z[cval + 1] += v[5] * x1 + v[6] * x2 + v[7] * x3 + v[8] * x4 + v[9] * x5; in MatMultAdd_SeqSBAIJ_5()
1020 z[cval + 2] += v[10] * x1 + v[11] * x2 + v[12] * x3 + v[13] * x4 + v[14] * x5; in MatMultAdd_SeqSBAIJ_5()
1021 z[cval + 3] += v[15] * x1 + v[16] * x2 + v[17] * x3 + v[18] * x4 + v[19] * x5; in MatMultAdd_SeqSBAIJ_5()
1022 z[cval + 4] += v[20] * x1 + v[21] * x2 + v[22] * x3 + v[23] * x4 + v[24] * x5; in MatMultAdd_SeqSBAIJ_5()
1024 …z[5 * i] += v[0] * x[cval] + v[5] * x[cval + 1] + v[10] * x[cval + 2] + v[15] * x[cval + 3] + v[20… in MatMultAdd_SeqSBAIJ_5()
1025 …z[5 * i + 1] += v[1] * x[cval] + v[6] * x[cval + 1] + v[11] * x[cval + 2] + v[16] * x[cval + 3] + … in MatMultAdd_SeqSBAIJ_5()
1026 …z[5 * i + 2] += v[2] * x[cval] + v[7] * x[cval + 1] + v[12] * x[cval + 2] + v[17] * x[cval + 3] + … in MatMultAdd_SeqSBAIJ_5()
1027 …z[5 * i + 3] += v[3] * x[cval] + v[8] * x[cval + 1] + v[13] * x[cval + 2] + v[18] * x[cval + 3] + … in MatMultAdd_SeqSBAIJ_5()
1028 …z[5 * i + 4] += v[4] * x[cval] + v[9] * x[cval + 1] + v[14] * x[cval + 2] + v[19] * x[cval + 3] + … in MatMultAdd_SeqSBAIJ_5()
1029 v += 25; in MatMultAdd_SeqSBAIJ_5()
1047 const MatScalar *v; in MatMultAdd_SeqSBAIJ_6() local
1058 v = a->a; in MatMultAdd_SeqSBAIJ_6()
1073 z[6 * i] += v[0] * x1 + v[6] * x2 + v[12] * x3 + v[18] * x4 + v[24] * x5 + v[30] * x6; in MatMultAdd_SeqSBAIJ_6()
1074 z[6 * i + 1] += v[6] * x1 + v[7] * x2 + v[13] * x3 + v[19] * x4 + v[25] * x5 + v[31] * x6; in MatMultAdd_SeqSBAIJ_6()
1075 z[6 * i + 2] += v[12] * x1 + v[13] * x2 + v[14] * x3 + v[20] * x4 + v[26] * x5 + v[32] * x6; in MatMultAdd_SeqSBAIJ_6()
1076 z[6 * i + 3] += v[18] * x1 + v[19] * x2 + v[20] * x3 + v[21] * x4 + v[27] * x5 + v[33] * x6; in MatMultAdd_SeqSBAIJ_6()
1077 z[6 * i + 4] += v[24] * x1 + v[25] * x2 + v[26] * x3 + v[27] * x4 + v[28] * x5 + v[34] * x6; in MatMultAdd_SeqSBAIJ_6()
1078 z[6 * i + 5] += v[30] * x1 + v[31] * x2 + v[32] * x3 + v[33] * x4 + v[34] * x5 + v[35] * x6; in MatMultAdd_SeqSBAIJ_6()
1079 v += 36; in MatMultAdd_SeqSBAIJ_6()
1083 … PetscPrefetchBlock(v + 36 * n, 36 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMultAdd_SeqSBAIJ_6()
1087 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5 + v[5] * x6; in MatMultAdd_SeqSBAIJ_6()
1088 z[cval + 1] += v[6] * x1 + v[7] * x2 + v[8] * x3 + v[9] * x4 + v[10] * x5 + v[11] * x6; in MatMultAdd_SeqSBAIJ_6()
1089 z[cval + 2] += v[12] * x1 + v[13] * x2 + v[14] * x3 + v[15] * x4 + v[16] * x5 + v[17] * x6; in MatMultAdd_SeqSBAIJ_6()
1090 z[cval + 3] += v[18] * x1 + v[19] * x2 + v[20] * x3 + v[21] * x4 + v[22] * x5 + v[23] * x6; in MatMultAdd_SeqSBAIJ_6()
1091 z[cval + 4] += v[24] * x1 + v[25] * x2 + v[26] * x3 + v[27] * x4 + v[28] * x5 + v[29] * x6; in MatMultAdd_SeqSBAIJ_6()
1092 z[cval + 5] += v[30] * x1 + v[31] * x2 + v[32] * x3 + v[33] * x4 + v[34] * x5 + v[35] * x6; in MatMultAdd_SeqSBAIJ_6()
1094 …z[6 * i] += v[0] * x[cval] + v[6] * x[cval + 1] + v[12] * x[cval + 2] + v[18] * x[cval + 3] + v[24… in MatMultAdd_SeqSBAIJ_6()
1095 …z[6 * i + 1] += v[1] * x[cval] + v[7] * x[cval + 1] + v[13] * x[cval + 2] + v[19] * x[cval + 3] + … in MatMultAdd_SeqSBAIJ_6()
1096 …z[6 * i + 2] += v[2] * x[cval] + v[8] * x[cval + 1] + v[14] * x[cval + 2] + v[20] * x[cval + 3] + … in MatMultAdd_SeqSBAIJ_6()
1097 …z[6 * i + 3] += v[3] * x[cval] + v[9] * x[cval + 1] + v[15] * x[cval + 2] + v[21] * x[cval + 3] + … in MatMultAdd_SeqSBAIJ_6()
1098 …z[6 * i + 4] += v[4] * x[cval] + v[10] * x[cval + 1] + v[16] * x[cval + 2] + v[22] * x[cval + 3] +… in MatMultAdd_SeqSBAIJ_6()
1099 …z[6 * i + 5] += v[5] * x[cval] + v[11] * x[cval + 1] + v[17] * x[cval + 2] + v[23] * x[cval + 3] +… in MatMultAdd_SeqSBAIJ_6()
1100 v += 36; in MatMultAdd_SeqSBAIJ_6()
1118 const MatScalar *v; in MatMultAdd_SeqSBAIJ_7() local
1129 v = a->a; in MatMultAdd_SeqSBAIJ_7()
1145 …z[7 * i] += v[0] * x1 + v[7] * x2 + v[14] * x3 + v[21] * x4 + v[28] * x5 + v[35] * x6 + v[42] * x7; in MatMultAdd_SeqSBAIJ_7()
1146 …z[7 * i + 1] += v[7] * x1 + v[8] * x2 + v[15] * x3 + v[22] * x4 + v[29] * x5 + v[36] * x6 + v[43] … in MatMultAdd_SeqSBAIJ_7()
1147 …z[7 * i + 2] += v[14] * x1 + v[15] * x2 + v[16] * x3 + v[23] * x4 + v[30] * x5 + v[37] * x6 + v[44… in MatMultAdd_SeqSBAIJ_7()
1148 …z[7 * i + 3] += v[21] * x1 + v[22] * x2 + v[23] * x3 + v[24] * x4 + v[31] * x5 + v[38] * x6 + v[45… in MatMultAdd_SeqSBAIJ_7()
1149 …z[7 * i + 4] += v[28] * x1 + v[29] * x2 + v[30] * x3 + v[31] * x4 + v[32] * x5 + v[39] * x6 + v[46… in MatMultAdd_SeqSBAIJ_7()
1150 …z[7 * i + 5] += v[35] * x1 + v[36] * x2 + v[37] * x3 + v[38] * x4 + v[39] * x5 + v[40] * x6 + v[47… in MatMultAdd_SeqSBAIJ_7()
1151 …z[7 * i + 6] += v[42] * x1 + v[43] * x2 + v[44] * x3 + v[45] * x4 + v[46] * x5 + v[47] * x6 + v[48… in MatMultAdd_SeqSBAIJ_7()
1152 v += 49; in MatMultAdd_SeqSBAIJ_7()
1156 … PetscPrefetchBlock(v + 49 * n, 49 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMultAdd_SeqSBAIJ_7()
1160 z[cval] += v[0] * x1 + v[1] * x2 + v[2] * x3 + v[3] * x4 + v[4] * x5 + v[5] * x6 + v[6] * x7; in MatMultAdd_SeqSBAIJ_7()
1161 …z[cval + 1] += v[7] * x1 + v[8] * x2 + v[9] * x3 + v[10] * x4 + v[11] * x5 + v[12] * x6 + v[13] * … in MatMultAdd_SeqSBAIJ_7()
1162 …z[cval + 2] += v[14] * x1 + v[15] * x2 + v[16] * x3 + v[17] * x4 + v[18] * x5 + v[19] * x6 + v[20]… in MatMultAdd_SeqSBAIJ_7()
1163 …z[cval + 3] += v[21] * x1 + v[22] * x2 + v[23] * x3 + v[24] * x4 + v[25] * x5 + v[26] * x6 + v[27]… in MatMultAdd_SeqSBAIJ_7()
1164 …z[cval + 4] += v[28] * x1 + v[29] * x2 + v[30] * x3 + v[31] * x4 + v[32] * x5 + v[33] * x6 + v[34]… in MatMultAdd_SeqSBAIJ_7()
1165 …z[cval + 5] += v[35] * x1 + v[36] * x2 + v[37] * x3 + v[38] * x4 + v[39] * x5 + v[40] * x6 + v[41]… in MatMultAdd_SeqSBAIJ_7()
1166 …z[cval + 6] += v[42] * x1 + v[43] * x2 + v[44] * x3 + v[45] * x4 + v[46] * x5 + v[47] * x6 + v[48]… in MatMultAdd_SeqSBAIJ_7()
1168 …z[7 * i] += v[0] * x[cval] + v[7] * x[cval + 1] + v[14] * x[cval + 2] + v[21] * x[cval + 3] + v[28… in MatMultAdd_SeqSBAIJ_7()
1169 … * i + 1] += v[1] * x[cval] + v[8] * x[cval + 1] + v[15] * x[cval + 2] + v[22] * x[cval + 3] + v[2… in MatMultAdd_SeqSBAIJ_7()
1170 … * i + 2] += v[2] * x[cval] + v[9] * x[cval + 1] + v[16] * x[cval + 2] + v[23] * x[cval + 3] + v[3… in MatMultAdd_SeqSBAIJ_7()
1171 …* i + 3] += v[3] * x[cval] + v[10] * x[cval + 1] + v[17] * x[cval + 2] + v[24] * x[cval + 3] + v[3… in MatMultAdd_SeqSBAIJ_7()
1172 …* i + 4] += v[4] * x[cval] + v[11] * x[cval + 1] + v[18] * x[cval + 2] + v[25] * x[cval + 3] + v[3… in MatMultAdd_SeqSBAIJ_7()
1173 …* i + 5] += v[5] * x[cval] + v[12] * x[cval + 1] + v[19] * x[cval + 2] + v[26] * x[cval + 3] + v[3… in MatMultAdd_SeqSBAIJ_7()
1174 …* i + 6] += v[6] * x[cval] + v[13] * x[cval + 1] + v[20] * x[cval + 2] + v[27] * x[cval + 3] + v[3… in MatMultAdd_SeqSBAIJ_7()
1175 v += 49; in MatMultAdd_SeqSBAIJ_7()
1193 const MatScalar *v; in MatMultAdd_SeqSBAIJ_N() local
1207 v = a->a; in MatMultAdd_SeqSBAIJ_N()
1227 PetscKernel_w_gets_w_plus_Ar_times_v(bs, ncols, work, v, z); in MatMultAdd_SeqSBAIJ_N()
1233 v += bs2; in MatMultAdd_SeqSBAIJ_N()
1240 PetscKernel_w_gets_w_plus_trans_Ar_times_v(bs, ncols, x, v, workt); in MatMultAdd_SeqSBAIJ_N()
1249 v += n * bs2; in MatMultAdd_SeqSBAIJ_N()
1277 const MatScalar *v = a->a; in MatNorm_SeqSBAIJ() local
1294 sum_diag += PetscRealPart(PetscConj(*v) * (*v)); in MatNorm_SeqSBAIJ()
1295 v++; in MatNorm_SeqSBAIJ()
1301 sum_off += PetscRealPart(PetscConj(*v) * (*v)); in MatNorm_SeqSBAIJ()
1302 v++; in MatNorm_SeqSBAIJ()
1324 v = a->a + ik * bs2 + j * bs; in MatNorm_SeqSBAIJ()
1326 sum[j] += PetscAbsScalar(*v); in MatNorm_SeqSBAIJ()
1327 v++; in MatNorm_SeqSBAIJ()
1346 v = a->a + i * bs2 + j; in MatNorm_SeqSBAIJ()
1348 sum[j] += PetscAbsScalar(*v); in MatNorm_SeqSBAIJ()
1349 v += bs; in MatNorm_SeqSBAIJ()
1396 PetscErrorCode MatGetDiagonal_SeqSBAIJ(Mat A, Vec v) in MatGetDiagonal_SeqSBAIJ() argument
1410 PetscCall(VecGetArrayWrite(v, &x)); in MatGetDiagonal_SeqSBAIJ()
1412 PetscCall(VecRestoreArrayWrite(v, &x)); in MatGetDiagonal_SeqSBAIJ()
1416 PetscCall(VecGetLocalSize(v, &n)); in MatGetDiagonal_SeqSBAIJ()
1418 PetscCall(VecGetArrayWrite(v, &x)); in MatGetDiagonal_SeqSBAIJ()
1437 PetscCall(VecRestoreArrayWrite(v, &x)); in MatGetDiagonal_SeqSBAIJ()
1446 MatScalar *aa, *v; in MatDiagonalScale_SeqSBAIJ() local
1466 v = aa + bs2 * ai[i]; in MatDiagonalScale_SeqSBAIJ()
1471 for (tmp = 0; tmp < bs; tmp++) (*v++) *= li[tmp] * x; in MatDiagonalScale_SeqSBAIJ()
1513 PetscErrorCode MatGetRowMaxAbs_SeqSBAIJ(Mat A, Vec v, PetscInt idx[]) in MatGetRowMaxAbs_SeqSBAIJ() argument
1532 PetscCall(VecSet(v, 0.0)); in MatGetRowMaxAbs_SeqSBAIJ()
1533 PetscCall(VecGetArray(v, &x)); in MatGetRowMaxAbs_SeqSBAIJ()
1534 PetscCall(VecGetLocalSize(v, &n)); in MatGetRowMaxAbs_SeqSBAIJ()
1555 PetscCall(VecRestoreArray(v, &x)); in MatGetRowMaxAbs_SeqSBAIJ()
1573 const MatScalar *v = a->a, *vv; in MatMatMult_SeqSBAIJ_1_Private() local
1582 PetscPrefetchBlock(v + n, n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMatMult_SeqSBAIJ_1_Private()
1584 vv = v; in MatMatMult_SeqSBAIJ_1_Private()
1587 v = vv; in MatMatMult_SeqSBAIJ_1_Private()
1591 z[0 + k * cm] += v[0] * x1; in MatMatMult_SeqSBAIJ_1_Private()
1592 if (*idx != i) c[(*idx) + k * cm] += (aconj ? PetscConj(v[0]) : v[0]) * b[i + k * bm]; in MatMatMult_SeqSBAIJ_1_Private()
1593 v += 1; in MatMatMult_SeqSBAIJ_1_Private()
1608 const MatScalar *v = a->a, *vv; in MatMatMult_SeqSBAIJ_2_Private() local
1616 PetscPrefetchBlock(v + 4 * n, 4 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMatMult_SeqSBAIJ_2_Private()
1618 vv = v; in MatMatMult_SeqSBAIJ_2_Private()
1621 v = vv; in MatMatMult_SeqSBAIJ_2_Private()
1626 z[0 + k * cm] += v[0] * x1 + v[2] * x2; in MatMatMult_SeqSBAIJ_2_Private()
1627 z[1 + k * cm] += v[1] * x1 + v[3] * x2; in MatMatMult_SeqSBAIJ_2_Private()
1629 c[2 * (*idx) + 0 + k * cm] += v[0] * b[2 * i + k * bm] + v[1] * b[2 * i + 1 + k * bm]; in MatMatMult_SeqSBAIJ_2_Private()
1630 c[2 * (*idx) + 1 + k * cm] += v[2] * b[2 * i + k * bm] + v[3] * b[2 * i + 1 + k * bm]; in MatMatMult_SeqSBAIJ_2_Private()
1632 v += 4; in MatMatMult_SeqSBAIJ_2_Private()
1647 const MatScalar *v = a->a, *vv; in MatMatMult_SeqSBAIJ_3_Private() local
1655 PetscPrefetchBlock(v + 9 * n, 9 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMatMult_SeqSBAIJ_3_Private()
1657 vv = v; in MatMatMult_SeqSBAIJ_3_Private()
1660 v = vv; in MatMatMult_SeqSBAIJ_3_Private()
1666 z[0 + k * cm] += v[0] * x1 + v[3] * x2 + v[6] * x3; in MatMatMult_SeqSBAIJ_3_Private()
1667 z[1 + k * cm] += v[1] * x1 + v[4] * x2 + v[7] * x3; in MatMatMult_SeqSBAIJ_3_Private()
1668 z[2 + k * cm] += v[2] * x1 + v[5] * x2 + v[8] * x3; in MatMatMult_SeqSBAIJ_3_Private()
1670 …c[3 * (*idx) + 0 + k * cm] += v[0] * b[3 * i + k * bm] + v[3] * b[3 * i + 1 + k * bm] + v[6] * b[3… in MatMatMult_SeqSBAIJ_3_Private()
1671 …c[3 * (*idx) + 1 + k * cm] += v[1] * b[3 * i + k * bm] + v[4] * b[3 * i + 1 + k * bm] + v[7] * b[3… in MatMatMult_SeqSBAIJ_3_Private()
1672 …c[3 * (*idx) + 2 + k * cm] += v[2] * b[3 * i + k * bm] + v[5] * b[3 * i + 1 + k * bm] + v[8] * b[3… in MatMatMult_SeqSBAIJ_3_Private()
1674 v += 9; in MatMatMult_SeqSBAIJ_3_Private()
1689 const MatScalar *v = a->a, *vv; in MatMatMult_SeqSBAIJ_4_Private() local
1697 … PetscPrefetchBlock(v + 16 * n, 16 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMatMult_SeqSBAIJ_4_Private()
1699 vv = v; in MatMatMult_SeqSBAIJ_4_Private()
1702 v = vv; in MatMatMult_SeqSBAIJ_4_Private()
1709 z[0 + k * cm] += v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4; in MatMatMult_SeqSBAIJ_4_Private()
1710 z[1 + k * cm] += v[1] * x1 + v[5] * x2 + v[9] * x3 + v[13] * x4; in MatMatMult_SeqSBAIJ_4_Private()
1711 z[2 + k * cm] += v[2] * x1 + v[6] * x2 + v[10] * x3 + v[14] * x4; in MatMatMult_SeqSBAIJ_4_Private()
1712 z[3 + k * cm] += v[3] * x1 + v[7] * x2 + v[11] * x3 + v[15] * x4; in MatMatMult_SeqSBAIJ_4_Private()
1714 …c[4 * (*idx) + 0 + k * cm] += v[0] * b[4 * i + k * bm] + v[4] * b[4 * i + 1 + k * bm] + v[8] * b[4… in MatMatMult_SeqSBAIJ_4_Private()
1715 …c[4 * (*idx) + 1 + k * cm] += v[1] * b[4 * i + k * bm] + v[5] * b[4 * i + 1 + k * bm] + v[9] * b[4… in MatMatMult_SeqSBAIJ_4_Private()
1716 …c[4 * (*idx) + 2 + k * cm] += v[2] * b[4 * i + k * bm] + v[6] * b[4 * i + 1 + k * bm] + v[10] * b[… in MatMatMult_SeqSBAIJ_4_Private()
1717 …c[4 * (*idx) + 3 + k * cm] += v[3] * b[4 * i + k * bm] + v[7] * b[4 * i + 1 + k * bm] + v[11] * b[… in MatMatMult_SeqSBAIJ_4_Private()
1719 v += 16; in MatMatMult_SeqSBAIJ_4_Private()
1734 const MatScalar *v = a->a, *vv; in MatMatMult_SeqSBAIJ_5_Private() local
1742 … PetscPrefetchBlock(v + 25 * n, 25 * n, 0, PETSC_PREFETCH_HINT_NTA); /* Entries for the next row */ in MatMatMult_SeqSBAIJ_5_Private()
1744 vv = v; in MatMatMult_SeqSBAIJ_5_Private()
1747 v = vv; in MatMatMult_SeqSBAIJ_5_Private()
1755 z[0 + k * cm] += v[0] * x1 + v[5] * x2 + v[10] * x3 + v[15] * x4 + v[20] * x5; in MatMatMult_SeqSBAIJ_5_Private()
1756 z[1 + k * cm] += v[1] * x1 + v[6] * x2 + v[11] * x3 + v[16] * x4 + v[21] * x5; in MatMatMult_SeqSBAIJ_5_Private()
1757 z[2 + k * cm] += v[2] * x1 + v[7] * x2 + v[12] * x3 + v[17] * x4 + v[22] * x5; in MatMatMult_SeqSBAIJ_5_Private()
1758 z[3 + k * cm] += v[3] * x1 + v[8] * x2 + v[13] * x3 + v[18] * x4 + v[23] * x5; in MatMatMult_SeqSBAIJ_5_Private()
1759 z[4 + k * cm] += v[4] * x1 + v[9] * x2 + v[14] * x3 + v[19] * x4 + v[24] * x5; in MatMatMult_SeqSBAIJ_5_Private()
1761 …0 + k * cm] += v[0] * b[5 * i + k * bm] + v[5] * b[5 * i + 1 + k * bm] + v[10] * b[5 * i + 2 + k *… in MatMatMult_SeqSBAIJ_5_Private()
1762 …1 + k * cm] += v[1] * b[5 * i + k * bm] + v[6] * b[5 * i + 1 + k * bm] + v[11] * b[5 * i + 2 + k *… in MatMatMult_SeqSBAIJ_5_Private()
1763 …2 + k * cm] += v[2] * b[5 * i + k * bm] + v[7] * b[5 * i + 1 + k * bm] + v[12] * b[5 * i + 2 + k *… in MatMatMult_SeqSBAIJ_5_Private()
1764 …3 + k * cm] += v[3] * b[5 * i + k * bm] + v[8] * b[5 * i + 1 + k * bm] + v[13] * b[5 * i + 2 + k *… in MatMatMult_SeqSBAIJ_5_Private()
1765 …4 + k * cm] += v[4] * b[5 * i + k * bm] + v[9] * b[5 * i + 1 + k * bm] + v[14] * b[5 * i + 2 + k *… in MatMatMult_SeqSBAIJ_5_Private()
1767 v += 25; in MatMatMult_SeqSBAIJ_5_Private()
1786 const MatScalar *v; in MatMatMultNumeric_SeqSBAIJ_SeqDense() local
1795 b = bd->v; in MatMatMultNumeric_SeqSBAIJ_SeqDense()
1820 v = a->a; in MatMatMultNumeric_SeqSBAIJ_SeqDense()
1828 …if (*idx != i) PetscCallBLAS("BLASgemm", BLASgemm_("T", "N", &bbs, &bcn, &bbs, &_DOne, v, &bbs, b … in MatMatMultNumeric_SeqSBAIJ_SeqDense()
1829 …PetscCallBLAS("BLASgemm", BLASgemm_("N", "N", &bbs, &bcn, &bbs, &_DOne, v, &bbs, b + bs * (*idx++)… in MatMatMultNumeric_SeqSBAIJ_SeqDense()
1830 v += bs2; in MatMatMultNumeric_SeqSBAIJ_SeqDense()