Lines Matching refs:z
290 PetscScalar *z, x1, x2, zero = 0.0; in MatMult_SeqSBAIJ_2() local
301 PetscCall(VecGetArray(zz, &z)); 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()
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()
336 PetscCall(VecRestoreArray(zz, &z)); in MatMult_SeqSBAIJ_2()
344 PetscScalar *z, x1, x2, x3, zero = 0.0; in MatMult_SeqSBAIJ_3() local
355 PetscCall(VecGetArray(zz, &z)); 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()
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()
394 PetscCall(VecRestoreArray(zz, &z)); in MatMult_SeqSBAIJ_3()
402 PetscScalar *z, x1, x2, x3, x4, zero = 0.0; in MatMult_SeqSBAIJ_4() local
413 PetscCall(VecGetArray(zz, &z)); 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()
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()
456 PetscCall(VecRestoreArray(zz, &z)); in MatMult_SeqSBAIJ_4()
464 PetscScalar *z, x1, x2, x3, x4, x5, zero = 0.0; in MatMult_SeqSBAIJ_5() local
475 PetscCall(VecGetArray(zz, &z)); 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()
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()
522 PetscCall(VecRestoreArray(zz, &z)); in MatMult_SeqSBAIJ_5()
530 PetscScalar *z, x1, x2, x3, x4, x5, x6, zero = 0.0; in MatMult_SeqSBAIJ_6() local
541 PetscCall(VecGetArray(zz, &z)); 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()
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()
592 PetscCall(VecRestoreArray(zz, &z)); in MatMult_SeqSBAIJ_6()
600 PetscScalar *z, x1, x2, x3, x4, x5, x6, x7, zero = 0.0; in MatMult_SeqSBAIJ_7() local
611 PetscCall(VecGetArray(zz, &z)); 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()
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 …z[7 * i + 1] += v[1] * x[cval] + v[8] * x[cval + 1] + v[15] * x[cval + 2] + v[22] * x[cval + 3] + … in MatMult_SeqSBAIJ_7()
654 …z[7 * i + 2] += v[2] * x[cval] + v[9] * x[cval + 1] + v[16] * x[cval + 2] + v[23] * x[cval + 3] + … in MatMult_SeqSBAIJ_7()
655 …z[7 * i + 3] += v[3] * x[cval] + v[10] * x[cval + 1] + v[17] * x[cval + 2] + v[24] * x[cval + 3] +… in MatMult_SeqSBAIJ_7()
656 …z[7 * i + 4] += v[4] * x[cval] + v[11] * x[cval + 1] + v[18] * x[cval + 2] + v[25] * x[cval + 3] +… in MatMult_SeqSBAIJ_7()
657 …z[7 * i + 5] += v[5] * x[cval] + v[12] * x[cval + 1] + v[19] * x[cval + 2] + v[26] * x[cval + 3] +… in MatMult_SeqSBAIJ_7()
658 …z[7 * i + 6] += v[6] * x[cval] + v[13] * x[cval + 1] + v[20] * x[cval + 2] + v[27] * x[cval + 3] +… in MatMult_SeqSBAIJ_7()
665 PetscCall(VecRestoreArray(zz, &z)); in MatMult_SeqSBAIJ_7()
676 PetscScalar *z, *z_ptr, *zb, *work, *workt, zero = 0.0; in MatMult_SeqSBAIJ_N() local
687 PetscCall(VecGetArray(zz, &z)); in MatMult_SeqSBAIJ_N()
690 z_ptr = z; in MatMult_SeqSBAIJ_N()
713 PetscKernel_w_gets_w_plus_Ar_times_v(bs, ncols, work, v, z); in MatMult_SeqSBAIJ_N()
736 z += bs; in MatMult_SeqSBAIJ_N()
741 PetscCall(VecRestoreArray(zz, &z)); in MatMult_SeqSBAIJ_N()
749 PetscScalar *z, x1; in MatMultAdd_SeqSBAIJ_1() local
761 PetscCall(VecGetArray(zz, &z)); 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()
793 PetscCall(VecRestoreArray(zz, &z)); in MatMultAdd_SeqSBAIJ_1()
802 PetscScalar *z, x1, x2; in MatMultAdd_SeqSBAIJ_2() local
813 PetscCall(VecGetArray(zz, &z)); 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()
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()
847 PetscCall(VecRestoreArray(zz, &z)); in MatMultAdd_SeqSBAIJ_2()
856 PetscScalar *z, x1, x2, x3; in MatMultAdd_SeqSBAIJ_3() local
867 PetscCall(VecGetArray(zz, &z)); 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()
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()
906 PetscCall(VecRestoreArray(zz, &z)); in MatMultAdd_SeqSBAIJ_3()
915 PetscScalar *z, x1, x2, x3, x4; in MatMultAdd_SeqSBAIJ_4() local
926 PetscCall(VecGetArray(zz, &z)); 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()
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()
969 PetscCall(VecRestoreArray(zz, &z)); in MatMultAdd_SeqSBAIJ_4()
978 PetscScalar *z, x1, x2, x3, x4, x5; in MatMultAdd_SeqSBAIJ_5() local
989 PetscCall(VecGetArray(zz, &z)); 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()
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()
1036 PetscCall(VecRestoreArray(zz, &z)); in MatMultAdd_SeqSBAIJ_5()
1045 PetscScalar *z, x1, x2, x3, x4, x5, x6; in MatMultAdd_SeqSBAIJ_6() local
1056 PetscCall(VecGetArray(zz, &z)); 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()
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()
1107 PetscCall(VecRestoreArray(zz, &z)); in MatMultAdd_SeqSBAIJ_6()
1116 PetscScalar *z, x1, x2, x3, x4, x5, x6, x7; in MatMultAdd_SeqSBAIJ_7() local
1127 PetscCall(VecGetArray(zz, &z)); 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()
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 …z[7 * i + 1] += v[1] * x[cval] + v[8] * x[cval + 1] + v[15] * x[cval + 2] + v[22] * x[cval + 3] + … in MatMultAdd_SeqSBAIJ_7()
1170 …z[7 * i + 2] += v[2] * x[cval] + v[9] * x[cval + 1] + v[16] * x[cval + 2] + v[23] * x[cval + 3] + … in MatMultAdd_SeqSBAIJ_7()
1171 …z[7 * i + 3] += v[3] * x[cval] + v[10] * x[cval + 1] + v[17] * x[cval + 2] + v[24] * x[cval + 3] +… in MatMultAdd_SeqSBAIJ_7()
1172 …z[7 * i + 4] += v[4] * x[cval] + v[11] * x[cval + 1] + v[18] * x[cval + 2] + v[25] * x[cval + 3] +… in MatMultAdd_SeqSBAIJ_7()
1173 …z[7 * i + 5] += v[5] * x[cval] + v[12] * x[cval + 1] + v[19] * x[cval + 2] + v[26] * x[cval + 3] +… in MatMultAdd_SeqSBAIJ_7()
1174 …z[7 * i + 6] += v[6] * x[cval] + v[13] * x[cval + 1] + v[20] * x[cval + 2] + v[27] * x[cval + 3] +… in MatMultAdd_SeqSBAIJ_7()
1182 PetscCall(VecRestoreArray(zz, &z)); in MatMultAdd_SeqSBAIJ_7()
1191 PetscScalar *z, *z_ptr = NULL, *zb, *work, *workt; in MatMultAdd_SeqSBAIJ_N() local
1203 PetscCall(VecGetArray(zz, &z)); in MatMultAdd_SeqSBAIJ_N()
1204 z_ptr = z; in MatMultAdd_SeqSBAIJ_N()
1227 PetscKernel_w_gets_w_plus_Ar_times_v(bs, ncols, work, v, z); in MatMultAdd_SeqSBAIJ_N()
1250 z += bs; in MatMultAdd_SeqSBAIJ_N()
1255 PetscCall(VecRestoreArray(zz, &z)); in MatMultAdd_SeqSBAIJ_N()
1570 PetscScalar *z = c; in MatMatMult_SeqSBAIJ_1_Private() local
1591 z[0 + k * cm] += v[0] * x1; in MatMatMult_SeqSBAIJ_1_Private()
1597 z += 1; in MatMatMult_SeqSBAIJ_1_Private()
1605 PetscScalar *z = c; in MatMatMult_SeqSBAIJ_2_Private() local
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()
1636 z += 2; in MatMatMult_SeqSBAIJ_2_Private()
1644 PetscScalar *z = c; in MatMatMult_SeqSBAIJ_3_Private() local
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()
1678 z += 3; in MatMatMult_SeqSBAIJ_3_Private()
1686 PetscScalar *z = c; in MatMatMult_SeqSBAIJ_4_Private() local
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()
1723 z += 4; in MatMatMult_SeqSBAIJ_4_Private()
1731 PetscScalar *z = c; in MatMatMult_SeqSBAIJ_5_Private() local
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()
1771 z += 5; in MatMatMult_SeqSBAIJ_5_Private()
1784 PetscScalar *z = NULL; in MatMatMultNumeric_SeqSBAIJ_SeqDense() local
1823 z = c; in MatMatMultNumeric_SeqSBAIJ_SeqDense()
1829 … BLASgemm_("N", "N", &bbs, &bcn, &bbs, &_DOne, v, &bbs, b + bs * (*idx++), &bbm, &_DOne, z, &bcm)); in MatMatMultNumeric_SeqSBAIJ_SeqDense()
1832 z += bs; in MatMatMultNumeric_SeqSBAIJ_SeqDense()