Lines Matching refs:uik
12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() local
25 PetscCall(PetscMalloc2(9, &dk, 9, &uik)); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
58 uik[0] = -(diag[0] * u[0] + diag[3] * u[1] + diag[6] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
59 uik[1] = -(diag[1] * u[0] + diag[4] * u[1] + diag[7] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
60 uik[2] = -(diag[2] * u[0] + diag[5] * u[1] + diag[8] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
62 uik[3] = -(diag[0] * u[3] + diag[3] * u[4] + diag[6] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
63 uik[4] = -(diag[1] * u[3] + diag[4] * u[4] + diag[7] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
64 uik[5] = -(diag[2] * u[3] + diag[5] * u[4] + diag[8] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
66 uik[6] = -(diag[0] * u[6] + diag[3] * u[7] + diag[6] * u[8]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
67 uik[7] = -(diag[1] * u[6] + diag[4] * u[7] + diag[7] * u[8]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
68 uik[8] = -(diag[2] * u[6] + diag[5] * u[7] + diag[8] * u[8]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
71 dk[0] += uik[0] * u[0] + uik[1] * u[1] + uik[2] * u[2]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
72 dk[1] += uik[3] * u[0] + uik[4] * u[1] + uik[5] * u[2]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
73 dk[2] += uik[6] * u[0] + uik[7] * u[1] + uik[8] * u[2]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
75 dk[3] += uik[0] * u[3] + uik[1] * u[4] + uik[2] * u[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
76 dk[4] += uik[3] * u[3] + uik[4] * u[4] + uik[5] * u[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
77 dk[5] += uik[6] * u[3] + uik[7] * u[4] + uik[8] * u[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
79 dk[6] += uik[0] * u[6] + uik[1] * u[7] + uik[2] * u[8]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
80 dk[7] += uik[3] * u[6] + uik[4] * u[7] + uik[5] * u[8]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
81 dk[8] += uik[6] * u[6] + uik[7] * u[7] + uik[8] * u[8]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
86 PetscCall(PetscArraycpy(ba + ili * 9, uik, 9)); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
96 rtmp_ptr[0] += uik[0] * u[0] + uik[1] * u[1] + uik[2] * u[2]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
97 rtmp_ptr[1] += uik[3] * u[0] + uik[4] * u[1] + uik[5] * u[2]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
98 rtmp_ptr[2] += uik[6] * u[0] + uik[7] * u[1] + uik[8] * u[2]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
100 rtmp_ptr[3] += uik[0] * u[3] + uik[1] * u[4] + uik[2] * u[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
101 rtmp_ptr[4] += uik[3] * u[3] + uik[4] * u[4] + uik[5] * u[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
102 rtmp_ptr[5] += uik[6] * u[3] + uik[7] * u[4] + uik[8] * u[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
104 rtmp_ptr[6] += uik[0] * u[6] + uik[1] * u[7] + uik[2] * u[8]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
105 rtmp_ptr[7] += uik[3] * u[6] + uik[4] * u[7] + uik[5] * u[8]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
106 rtmp_ptr[8] += uik[6] * u[6] + uik[7] * u[7] + uik[8] * u[8]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
150 PetscCall(PetscFree2(dk, uik)); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()