Lines Matching refs:ibdiag
2423 MatScalar *ibdiag, *bdiag, work[25]; in MatInvertDiagonalForSOR_SeqAIJ_Inode() local
2434 if (!a->inode.ibdiag) { in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2441 …PetscCall(PetscMalloc3(cnt, &a->inode.ibdiag, cnt, &a->inode.bdiag, A->rmap->n, &a->inode.ssor_wor… in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2445 ibdiag = a->inode.ibdiag; in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2453 PetscCall(PetscArraycpy(ibdiag + cnt, bdiag + cnt, nodesz * nodesz)); in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2458 if (PetscAbsScalar(ibdiag[cnt]) < zeropivot) { in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2461 A->factorerror_zeropivot_value = PetscAbsScalar(ibdiag[cnt]); in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2465 ibdiag[cnt] = 1.0 / ibdiag[cnt]; in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2468 …PetscCall(PetscKernel_A_gets_inverse_A_2(ibdiag + cnt, shift, allowzeropivot, &zeropivotdetected)); in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2472 …PetscCall(PetscKernel_A_gets_inverse_A_3(ibdiag + cnt, shift, allowzeropivot, &zeropivotdetected)); in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2476 …PetscCall(PetscKernel_A_gets_inverse_A_4(ibdiag + cnt, shift, allowzeropivot, &zeropivotdetected)); in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2480 …PetscCall(PetscKernel_A_gets_inverse_A_5(ibdiag + cnt, ipvt, work, shift, allowzeropivot, &zeropiv… in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2497 MatScalar *ibdiag, *bdiag, *t; in MatSOR_SeqAIJ_Inode() local
2512 ibdiag = a->inode.ibdiag; in MatSOR_SeqAIJ_Inode()
2547 x[row++] = sum1 * (*ibdiag++); in MatSOR_SeqAIJ_Inode()
2572 x[row++] = sum1 * ibdiag[0] + sum2 * ibdiag[2]; in MatSOR_SeqAIJ_Inode()
2573 x[row++] = sum1 * ibdiag[1] + sum2 * ibdiag[3]; in MatSOR_SeqAIJ_Inode()
2574 ibdiag += 4; in MatSOR_SeqAIJ_Inode()
2605 x[row++] = sum1 * ibdiag[0] + sum2 * ibdiag[3] + sum3 * ibdiag[6]; in MatSOR_SeqAIJ_Inode()
2606 x[row++] = sum1 * ibdiag[1] + sum2 * ibdiag[4] + sum3 * ibdiag[7]; in MatSOR_SeqAIJ_Inode()
2607 x[row++] = sum1 * ibdiag[2] + sum2 * ibdiag[5] + sum3 * ibdiag[8]; in MatSOR_SeqAIJ_Inode()
2608 ibdiag += 9; in MatSOR_SeqAIJ_Inode()
2645 x[row++] = sum1 * ibdiag[0] + sum2 * ibdiag[4] + sum3 * ibdiag[8] + sum4 * ibdiag[12]; in MatSOR_SeqAIJ_Inode()
2646 x[row++] = sum1 * ibdiag[1] + sum2 * ibdiag[5] + sum3 * ibdiag[9] + sum4 * ibdiag[13]; in MatSOR_SeqAIJ_Inode()
2647 x[row++] = sum1 * ibdiag[2] + sum2 * ibdiag[6] + sum3 * ibdiag[10] + sum4 * ibdiag[14]; in MatSOR_SeqAIJ_Inode()
2648 x[row++] = sum1 * ibdiag[3] + sum2 * ibdiag[7] + sum3 * ibdiag[11] + sum4 * ibdiag[15]; in MatSOR_SeqAIJ_Inode()
2649 ibdiag += 16; in MatSOR_SeqAIJ_Inode()
2692 …x[row++] = sum1 * ibdiag[0] + sum2 * ibdiag[5] + sum3 * ibdiag[10] + sum4 * ibdiag[15] + sum5 * … in MatSOR_SeqAIJ_Inode()
2693 …x[row++] = sum1 * ibdiag[1] + sum2 * ibdiag[6] + sum3 * ibdiag[11] + sum4 * ibdiag[16] + sum5 * … in MatSOR_SeqAIJ_Inode()
2694 …x[row++] = sum1 * ibdiag[2] + sum2 * ibdiag[7] + sum3 * ibdiag[12] + sum4 * ibdiag[17] + sum5 * … in MatSOR_SeqAIJ_Inode()
2695 …x[row++] = sum1 * ibdiag[3] + sum2 * ibdiag[8] + sum3 * ibdiag[13] + sum4 * ibdiag[18] + sum5 * … in MatSOR_SeqAIJ_Inode()
2696 …x[row++] = sum1 * ibdiag[4] + sum2 * ibdiag[9] + sum3 * ibdiag[14] + sum4 * ibdiag[19] + sum5 * … in MatSOR_SeqAIJ_Inode()
2697 ibdiag += 25; in MatSOR_SeqAIJ_Inode()
2708 ibdiag = a->inode.ibdiag + a->inode.bdiagsize; in MatSOR_SeqAIJ_Inode()
2711 ibdiag -= nodesz * nodesz; in MatSOR_SeqAIJ_Inode()
2735 x[row--] = sum1 * (*ibdiag); in MatSOR_SeqAIJ_Inode()
2761 x[row--] = sum2 * ibdiag[1] + sum1 * ibdiag[3]; in MatSOR_SeqAIJ_Inode()
2762 x[row--] = sum2 * ibdiag[0] + sum1 * ibdiag[2]; in MatSOR_SeqAIJ_Inode()
2791 x[row--] = sum3 * ibdiag[2] + sum2 * ibdiag[5] + sum1 * ibdiag[8]; in MatSOR_SeqAIJ_Inode()
2792 x[row--] = sum3 * ibdiag[1] + sum2 * ibdiag[4] + sum1 * ibdiag[7]; in MatSOR_SeqAIJ_Inode()
2793 x[row--] = sum3 * ibdiag[0] + sum2 * ibdiag[3] + sum1 * ibdiag[6]; in MatSOR_SeqAIJ_Inode()
2827 x[row--] = sum4 * ibdiag[3] + sum3 * ibdiag[7] + sum2 * ibdiag[11] + sum1 * ibdiag[15]; in MatSOR_SeqAIJ_Inode()
2828 x[row--] = sum4 * ibdiag[2] + sum3 * ibdiag[6] + sum2 * ibdiag[10] + sum1 * ibdiag[14]; in MatSOR_SeqAIJ_Inode()
2829 x[row--] = sum4 * ibdiag[1] + sum3 * ibdiag[5] + sum2 * ibdiag[9] + sum1 * ibdiag[13]; in MatSOR_SeqAIJ_Inode()
2830 x[row--] = sum4 * ibdiag[0] + sum3 * ibdiag[4] + sum2 * ibdiag[8] + sum1 * ibdiag[12]; in MatSOR_SeqAIJ_Inode()
2869 …x[row--] = sum5 * ibdiag[4] + sum4 * ibdiag[9] + sum3 * ibdiag[14] + sum2 * ibdiag[19] + sum1 * ib… in MatSOR_SeqAIJ_Inode()
2870 …x[row--] = sum5 * ibdiag[3] + sum4 * ibdiag[8] + sum3 * ibdiag[13] + sum2 * ibdiag[18] + sum1 * ib… in MatSOR_SeqAIJ_Inode()
2871 …x[row--] = sum5 * ibdiag[2] + sum4 * ibdiag[7] + sum3 * ibdiag[12] + sum2 * ibdiag[17] + sum1 * ib… in MatSOR_SeqAIJ_Inode()
2872 …x[row--] = sum5 * ibdiag[1] + sum4 * ibdiag[6] + sum3 * ibdiag[11] + sum2 * ibdiag[16] + sum1 * ib… in MatSOR_SeqAIJ_Inode()
2873 …x[row--] = sum5 * ibdiag[0] + sum4 * ibdiag[5] + sum3 * ibdiag[10] + sum2 * ibdiag[15] + sum1 * ib… in MatSOR_SeqAIJ_Inode()
2886 …for (i = 0, row = 0, ibdiag = a->inode.ibdiag; i < m; row += nodesz, ibdiag += nodesz * nodesz, i+… in MatSOR_SeqAIJ_Inode()
2935 x[row] = sum1 * (*ibdiag); in MatSOR_SeqAIJ_Inode()
2981 x[row] = sum1 * ibdiag[0] + sum2 * ibdiag[2]; in MatSOR_SeqAIJ_Inode()
2982 x[row + 1] = sum1 * ibdiag[1] + sum2 * ibdiag[3]; in MatSOR_SeqAIJ_Inode()
3039 x[row] = sum1 * ibdiag[0] + sum2 * ibdiag[3] + sum3 * ibdiag[6]; in MatSOR_SeqAIJ_Inode()
3040 x[row + 1] = sum1 * ibdiag[1] + sum2 * ibdiag[4] + sum3 * ibdiag[7]; in MatSOR_SeqAIJ_Inode()
3041 x[row + 2] = sum1 * ibdiag[2] + sum2 * ibdiag[5] + sum3 * ibdiag[8]; in MatSOR_SeqAIJ_Inode()
3109 x[row] = sum1 * ibdiag[0] + sum2 * ibdiag[4] + sum3 * ibdiag[8] + sum4 * ibdiag[12]; in MatSOR_SeqAIJ_Inode()
3110 x[row + 1] = sum1 * ibdiag[1] + sum2 * ibdiag[5] + sum3 * ibdiag[9] + sum4 * ibdiag[13]; in MatSOR_SeqAIJ_Inode()
3111 x[row + 2] = sum1 * ibdiag[2] + sum2 * ibdiag[6] + sum3 * ibdiag[10] + sum4 * ibdiag[14]; in MatSOR_SeqAIJ_Inode()
3112 x[row + 3] = sum1 * ibdiag[3] + sum2 * ibdiag[7] + sum3 * ibdiag[11] + sum4 * ibdiag[15]; in MatSOR_SeqAIJ_Inode()
3191 …x[row] = sum1 * ibdiag[0] + sum2 * ibdiag[5] + sum3 * ibdiag[10] + sum4 * ibdiag[15] + sum5 * … in MatSOR_SeqAIJ_Inode()
3192 …x[row + 1] = sum1 * ibdiag[1] + sum2 * ibdiag[6] + sum3 * ibdiag[11] + sum4 * ibdiag[16] + sum5 * … in MatSOR_SeqAIJ_Inode()
3193 …x[row + 2] = sum1 * ibdiag[2] + sum2 * ibdiag[7] + sum3 * ibdiag[12] + sum4 * ibdiag[17] + sum5 * … in MatSOR_SeqAIJ_Inode()
3194 …x[row + 3] = sum1 * ibdiag[3] + sum2 * ibdiag[8] + sum3 * ibdiag[13] + sum4 * ibdiag[18] + sum5 * … in MatSOR_SeqAIJ_Inode()
3195 …x[row + 4] = sum1 * ibdiag[4] + sum2 * ibdiag[9] + sum3 * ibdiag[14] + sum4 * ibdiag[19] + sum5 * … in MatSOR_SeqAIJ_Inode()
3206 ibdiag = a->inode.ibdiag + a->inode.bdiagsize; in MatSOR_SeqAIJ_Inode()
3209 ibdiag -= nodesz * nodesz; in MatSOR_SeqAIJ_Inode()
3308 …x[row--] += sum5 * ibdiag[4] + sum4 * ibdiag[9] + sum3 * ibdiag[14] + sum2 * ibdiag[19] + sum1 * i… in MatSOR_SeqAIJ_Inode()
3309 …x[row--] += sum5 * ibdiag[3] + sum4 * ibdiag[8] + sum3 * ibdiag[13] + sum2 * ibdiag[18] + sum1 * i… in MatSOR_SeqAIJ_Inode()
3310 …x[row--] += sum5 * ibdiag[2] + sum4 * ibdiag[7] + sum3 * ibdiag[12] + sum2 * ibdiag[17] + sum1 * i… in MatSOR_SeqAIJ_Inode()
3311 …x[row--] += sum5 * ibdiag[1] + sum4 * ibdiag[6] + sum3 * ibdiag[11] + sum2 * ibdiag[16] + sum1 * i… in MatSOR_SeqAIJ_Inode()
3312 …x[row--] += sum5 * ibdiag[0] + sum4 * ibdiag[5] + sum3 * ibdiag[10] + sum2 * ibdiag[15] + sum1 * i… in MatSOR_SeqAIJ_Inode()
3315 x[row--] += sum4 * ibdiag[3] + sum3 * ibdiag[7] + sum2 * ibdiag[11] + sum1 * ibdiag[15]; in MatSOR_SeqAIJ_Inode()
3316 x[row--] += sum4 * ibdiag[2] + sum3 * ibdiag[6] + sum2 * ibdiag[10] + sum1 * ibdiag[14]; in MatSOR_SeqAIJ_Inode()
3317 x[row--] += sum4 * ibdiag[1] + sum3 * ibdiag[5] + sum2 * ibdiag[9] + sum1 * ibdiag[13]; in MatSOR_SeqAIJ_Inode()
3318 x[row--] += sum4 * ibdiag[0] + sum3 * ibdiag[4] + sum2 * ibdiag[8] + sum1 * ibdiag[12]; in MatSOR_SeqAIJ_Inode()
3321 x[row--] += sum3 * ibdiag[2] + sum2 * ibdiag[5] + sum1 * ibdiag[8]; in MatSOR_SeqAIJ_Inode()
3322 x[row--] += sum3 * ibdiag[1] + sum2 * ibdiag[4] + sum1 * ibdiag[7]; in MatSOR_SeqAIJ_Inode()
3323 x[row--] += sum3 * ibdiag[0] + sum2 * ibdiag[3] + sum1 * ibdiag[6]; in MatSOR_SeqAIJ_Inode()
3326 x[row--] += sum2 * ibdiag[1] + sum1 * ibdiag[3]; in MatSOR_SeqAIJ_Inode()
3327 x[row--] += sum2 * ibdiag[0] + sum1 * ibdiag[2]; in MatSOR_SeqAIJ_Inode()
3330 x[row--] += sum1 * (*ibdiag); in MatSOR_SeqAIJ_Inode()
3337 …x[row--] = sum5 * ibdiag[4] + sum4 * ibdiag[9] + sum3 * ibdiag[14] + sum2 * ibdiag[19] + sum1 * ib… in MatSOR_SeqAIJ_Inode()
3338 …x[row--] = sum5 * ibdiag[3] + sum4 * ibdiag[8] + sum3 * ibdiag[13] + sum2 * ibdiag[18] + sum1 * ib… in MatSOR_SeqAIJ_Inode()
3339 …x[row--] = sum5 * ibdiag[2] + sum4 * ibdiag[7] + sum3 * ibdiag[12] + sum2 * ibdiag[17] + sum1 * ib… in MatSOR_SeqAIJ_Inode()
3340 …x[row--] = sum5 * ibdiag[1] + sum4 * ibdiag[6] + sum3 * ibdiag[11] + sum2 * ibdiag[16] + sum1 * ib… in MatSOR_SeqAIJ_Inode()
3341 …x[row--] = sum5 * ibdiag[0] + sum4 * ibdiag[5] + sum3 * ibdiag[10] + sum2 * ibdiag[15] + sum1 * ib… in MatSOR_SeqAIJ_Inode()
3344 x[row--] = sum4 * ibdiag[3] + sum3 * ibdiag[7] + sum2 * ibdiag[11] + sum1 * ibdiag[15]; in MatSOR_SeqAIJ_Inode()
3345 x[row--] = sum4 * ibdiag[2] + sum3 * ibdiag[6] + sum2 * ibdiag[10] + sum1 * ibdiag[14]; in MatSOR_SeqAIJ_Inode()
3346 x[row--] = sum4 * ibdiag[1] + sum3 * ibdiag[5] + sum2 * ibdiag[9] + sum1 * ibdiag[13]; in MatSOR_SeqAIJ_Inode()
3347 x[row--] = sum4 * ibdiag[0] + sum3 * ibdiag[4] + sum2 * ibdiag[8] + sum1 * ibdiag[12]; in MatSOR_SeqAIJ_Inode()
3350 x[row--] = sum3 * ibdiag[2] + sum2 * ibdiag[5] + sum1 * ibdiag[8]; in MatSOR_SeqAIJ_Inode()
3351 x[row--] = sum3 * ibdiag[1] + sum2 * ibdiag[4] + sum1 * ibdiag[7]; in MatSOR_SeqAIJ_Inode()
3352 x[row--] = sum3 * ibdiag[0] + sum2 * ibdiag[3] + sum1 * ibdiag[6]; in MatSOR_SeqAIJ_Inode()
3355 x[row--] = sum2 * ibdiag[1] + sum1 * ibdiag[3]; in MatSOR_SeqAIJ_Inode()
3356 x[row--] = sum2 * ibdiag[0] + sum1 * ibdiag[2]; in MatSOR_SeqAIJ_Inode()
3359 x[row--] = sum1 * (*ibdiag); in MatSOR_SeqAIJ_Inode()
3375 ibdiag = a->inode.ibdiag + a->inode.bdiagsize; in MatSOR_SeqAIJ_Inode()
3378 ibdiag -= nodesz * nodesz; in MatSOR_SeqAIJ_Inode()
3401 x[row] = sum1 * (*ibdiag); in MatSOR_SeqAIJ_Inode()
3428 x[row] = sum2 * ibdiag[1] + sum1 * ibdiag[3]; in MatSOR_SeqAIJ_Inode()
3429 x[row - 1] = sum2 * ibdiag[0] + sum1 * ibdiag[2]; in MatSOR_SeqAIJ_Inode()
3459 x[row] = sum3 * ibdiag[2] + sum2 * ibdiag[5] + sum1 * ibdiag[8]; in MatSOR_SeqAIJ_Inode()
3460 x[row - 1] = sum3 * ibdiag[1] + sum2 * ibdiag[4] + sum1 * ibdiag[7]; in MatSOR_SeqAIJ_Inode()
3461 x[row - 2] = sum3 * ibdiag[0] + sum2 * ibdiag[3] + sum1 * ibdiag[6]; in MatSOR_SeqAIJ_Inode()
3496 x[row] = sum4 * ibdiag[3] + sum3 * ibdiag[7] + sum2 * ibdiag[11] + sum1 * ibdiag[15]; in MatSOR_SeqAIJ_Inode()
3497 x[row - 1] = sum4 * ibdiag[2] + sum3 * ibdiag[6] + sum2 * ibdiag[10] + sum1 * ibdiag[14]; in MatSOR_SeqAIJ_Inode()
3498 x[row - 2] = sum4 * ibdiag[1] + sum3 * ibdiag[5] + sum2 * ibdiag[9] + sum1 * ibdiag[13]; in MatSOR_SeqAIJ_Inode()
3499 x[row - 3] = sum4 * ibdiag[0] + sum3 * ibdiag[4] + sum2 * ibdiag[8] + sum1 * ibdiag[12]; in MatSOR_SeqAIJ_Inode()
3539 …x[row] = sum5 * ibdiag[4] + sum4 * ibdiag[9] + sum3 * ibdiag[14] + sum2 * ibdiag[19] + sum1 * … in MatSOR_SeqAIJ_Inode()
3540 …x[row - 1] = sum5 * ibdiag[3] + sum4 * ibdiag[8] + sum3 * ibdiag[13] + sum2 * ibdiag[18] + sum1 * … in MatSOR_SeqAIJ_Inode()
3541 …x[row - 2] = sum5 * ibdiag[2] + sum4 * ibdiag[7] + sum3 * ibdiag[12] + sum2 * ibdiag[17] + sum1 * … in MatSOR_SeqAIJ_Inode()
3542 …x[row - 3] = sum5 * ibdiag[1] + sum4 * ibdiag[6] + sum3 * ibdiag[11] + sum2 * ibdiag[16] + sum1 * … in MatSOR_SeqAIJ_Inode()
3543 …x[row - 4] = sum5 * ibdiag[0] + sum4 * ibdiag[5] + sum3 * ibdiag[10] + sum2 * ibdiag[15] + sum1 * … in MatSOR_SeqAIJ_Inode()
3654 x[row] += t[row] = sum1 * (*ibdiag++); in MatSOR_SeqAIJ_Inode()
3678 x[row] += t[row] = sum1 * ibdiag[0] + sum2 * ibdiag[2]; in MatSOR_SeqAIJ_Inode()
3679 x[row + 1] += t[row + 1] = sum1 * ibdiag[1] + sum2 * ibdiag[3]; in MatSOR_SeqAIJ_Inode()
3680 ibdiag += 4; in MatSOR_SeqAIJ_Inode()
3709 x[row] += t[row] = sum1 * ibdiag[0] + sum2 * ibdiag[3] + sum3 * ibdiag[6]; in MatSOR_SeqAIJ_Inode()
3710 x[row + 1] += t[row + 1] = sum1 * ibdiag[1] + sum2 * ibdiag[4] + sum3 * ibdiag[7]; in MatSOR_SeqAIJ_Inode()
3711 x[row + 2] += t[row + 2] = sum1 * ibdiag[2] + sum2 * ibdiag[5] + sum3 * ibdiag[8]; in MatSOR_SeqAIJ_Inode()
3712 ibdiag += 9; in MatSOR_SeqAIJ_Inode()
3746 …x[row] += t[row] = sum1 * ibdiag[0] + sum2 * ibdiag[4] + sum3 * ibdiag[8] + sum4 * ibdiag[… in MatSOR_SeqAIJ_Inode()
3747 …x[row + 1] += t[row + 1] = sum1 * ibdiag[1] + sum2 * ibdiag[5] + sum3 * ibdiag[9] + sum4 * ibdiag[… in MatSOR_SeqAIJ_Inode()
3748 …x[row + 2] += t[row + 2] = sum1 * ibdiag[2] + sum2 * ibdiag[6] + sum3 * ibdiag[10] + sum4 * ibdiag… in MatSOR_SeqAIJ_Inode()
3749 …x[row + 3] += t[row + 3] = sum1 * ibdiag[3] + sum2 * ibdiag[7] + sum3 * ibdiag[11] + sum4 * ibdiag… in MatSOR_SeqAIJ_Inode()
3750 ibdiag += 16; in MatSOR_SeqAIJ_Inode()
3789 …x[row] += t[row] = sum1 * ibdiag[0] + sum2 * ibdiag[5] + sum3 * ibdiag[10] + sum4 * ibdiag… in MatSOR_SeqAIJ_Inode()
3790 …x[row + 1] += t[row + 1] = sum1 * ibdiag[1] + sum2 * ibdiag[6] + sum3 * ibdiag[11] + sum4 * ibdiag… in MatSOR_SeqAIJ_Inode()
3791 …x[row + 2] += t[row + 2] = sum1 * ibdiag[2] + sum2 * ibdiag[7] + sum3 * ibdiag[12] + sum4 * ibdiag… in MatSOR_SeqAIJ_Inode()
3792 …x[row + 3] += t[row + 3] = sum1 * ibdiag[3] + sum2 * ibdiag[8] + sum3 * ibdiag[13] + sum4 * ibdiag… in MatSOR_SeqAIJ_Inode()
3793 …x[row + 4] += t[row + 4] = sum1 * ibdiag[4] + sum2 * ibdiag[9] + sum3 * ibdiag[14] + sum4 * ibdiag… in MatSOR_SeqAIJ_Inode()
3794 ibdiag += 25; in MatSOR_SeqAIJ_Inode()
3996 c->inode.ibdiag = NULL; in MatDuplicate_SeqAIJ_Inode()