12   MatScalar    *ba = b->a, *aa, *ap, *dk, *uik;  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()  local
26   PetscCall(PetscMalloc2(16, &dk, 16, &uik));  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
46     PetscCall(PetscArraycpy(dk, rtmp + k * 16, 16));  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
80       dk[0] += uik[0] * u[0] + uik[1] * u[1] + uik[2] * u[2] + uik[3] * u[3];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
81       dk[1] += uik[4] * u[0] + uik[5] * u[1] + uik[6] * u[2] + uik[7] * u[3];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
82       dk[2] += uik[8] * u[0] + uik[9] * u[1] + uik[10] * u[2] + uik[11] * u[3];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
83       dk[3] += uik[12] * u[0] + uik[13] * u[1] + uik[14] * u[2] + uik[15] * u[3];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
85       dk[4] += uik[0] * u[4] + uik[1] * u[5] + uik[2] * u[6] + uik[3] * u[7];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
86       dk[5] += uik[4] * u[4] + uik[5] * u[5] + uik[6] * u[6] + uik[7] * u[7];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
87       dk[6] += uik[8] * u[4] + uik[9] * u[5] + uik[10] * u[6] + uik[11] * u[7];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
88       dk[7] += uik[12] * u[4] + uik[13] * u[5] + uik[14] * u[6] + uik[15] * u[7];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
90       dk[8] += uik[0] * u[8] + uik[1] * u[9] + uik[2] * u[10] + uik[3] * u[11];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
91       dk[9] += uik[4] * u[8] + uik[5] * u[9] + uik[6] * u[10] + uik[7] * u[11];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
92       dk[10] += uik[8] * u[8] + uik[9] * u[9] + uik[10] * u[10] + uik[11] * u[11];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
93       dk[11] += uik[12] * u[8] + uik[13] * u[9] + uik[14] * u[10] + uik[15] * u[11];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
95       dk[12] += uik[0] * u[12] + uik[1] * u[13] + uik[2] * u[14] + uik[3] * u[15];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
96       dk[13] += uik[4] * u[12] + uik[5] * u[13] + uik[6] * u[14] + uik[7] * u[15];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
97       dk[14] += uik[8] * u[12] + uik[9] * u[13] + uik[10] * u[14] + uik[11] * u[15];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
98       dk[15] += uik[12] * u[12] + uik[13] * u[13] + uik[14] * u[14] + uik[15] * u[15];  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
148     PetscCall(PetscArraycpy(diag, dk, 16));  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
179   PetscCall(PetscFree2(dk, uik));  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()