Lines Matching refs:ba
566 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_N() local
650 diag = ba + i * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
651 u = ba + ili * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
660 PetscCall(PetscArraycpy(ba + ili * bs2, uik, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
669 u = ba + j * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
686 diag = ba + k * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
697 u = ba + j * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
739 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering() local
782 diag = ba + i * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
783 u = ba + ili * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
792 PetscCall(PetscArraycpy(ba + ili * bs2, uik, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
801 u = ba + j * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
818 diag = ba + k * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
829 u = ba + j * bs2; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
872 MatScalar *ba = b->a, *aa, *ap; in MatCholeskyFactorNumeric_SeqSBAIJ_2() local
958 diag = ba + i * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
959 u = ba + ili * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
974 PetscCall(PetscArraycpy(ba + ili * 4, uik, 4)); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
983 u = ba + j * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
1003 diag = ba + k * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
1013 u = ba + j * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
1051 MatScalar *ba = b->a, *aa, *ap, dk[8], uik[8]; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering() local
1100 diag = ba + i * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1101 u = ba + ili * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1116 PetscCall(PetscArraycpy(ba + ili * 4, uik, 4)); in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1125 u = ba + j * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1145 diag = ba + k * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1155 u = ba + j * 4; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1196 MatScalar *rtmp, *ba = b->a, *bval, *aa, dk, uikdi; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace() local
1242 bval = ba + bi[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1261 uikdi = -ba[ili] * ba[bi[i]]; /* diagonal(k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1262 dk += uikdi * ba[ili]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1263 ba[ili] = uikdi; /* -U(i,k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1269 for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1301 ba[bi[k]] = 1.0 / dk; /* U(k,k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1307 ba[j] = rtmp[col]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_inplace()
1353 MatScalar *rtmp, *ba = b->a, *bval, *aa = a->a, dk, uikdi; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering() local
1407 bval = ba + bi[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1427 uikdi = -ba[ili] * ba[bdiag[i]]; /* diagonal(k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1428 dk += uikdi * ba[ili]; /* update diag[k] */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1429 ba[ili] = uikdi; /* -U(i,k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1435 for (j = jmin; j < jmax; j++) rtmp[bj[j]] += uikdi * ba[j]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1452 ba[j] = rtmp[col]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1453 rs += PetscAbsScalar(ba[j]); in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1468 ba[bdiag[k]] = 1.0 / dk; /* U(k,k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1505 MatScalar *rtmp, *ba = b->a, *aa = a->a, dk, uikdi, *aval, *bval; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace() local
1538 bval = ba + bi[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1555 uikdi = -ba[ili] * ba[bi[i]]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1556 dk += uikdi * ba[ili]; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1557 ba[ili] = uikdi; /* -U(i,k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1564 bval = ba + jmin; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1597 ba[bi[k]] = 1.0 / dk; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()
1602 bval = ba + jmin; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering_inplace()