13   MatScalar    *u, *diag, *rtmp, *rtmp_ptr;  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()  local
56       diag = ba + i * 16;  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
59       uik[0] = -(diag[0] * u[0] + diag[4] * u[1] + diag[8] * u[2] + diag[12] * u[3]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
60       uik[1] = -(diag[1] * u[0] + diag[5] * u[1] + diag[9] * u[2] + diag[13] * u[3]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
61       uik[2] = -(diag[2] * u[0] + diag[6] * u[1] + diag[10] * u[2] + diag[14] * u[3]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
62       uik[3] = -(diag[3] * u[0] + diag[7] * u[1] + diag[11] * u[2] + diag[15] * u[3]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
64       uik[4] = -(diag[0] * u[4] + diag[4] * u[5] + diag[8] * u[6] + diag[12] * u[7]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
65       uik[5] = -(diag[1] * u[4] + diag[5] * u[5] + diag[9] * u[6] + diag[13] * u[7]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
66       uik[6] = -(diag[2] * u[4] + diag[6] * u[5] + diag[10] * u[6] + diag[14] * u[7]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
67       uik[7] = -(diag[3] * u[4] + diag[7] * u[5] + diag[11] * u[6] + diag[15] * u[7]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
69       uik[8]  = -(diag[0] * u[8] + diag[4] * u[9] + diag[8] * u[10] + diag[12] * u[11]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
70       uik[9]  = -(diag[1] * u[8] + diag[5] * u[9] + diag[9] * u[10] + diag[13] * u[11]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
71       uik[10] = -(diag[2] * u[8] + diag[6] * u[9] + diag[10] * u[10] + diag[14] * u[11]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
72       uik[11] = -(diag[3] * u[8] + diag[7] * u[9] + diag[11] * u[10] + diag[15] * u[11]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
74       uik[12] = -(diag[0] * u[12] + diag[4] * u[13] + diag[8] * u[14] + diag[12] * u[15]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
75       uik[13] = -(diag[1] * u[12] + diag[5] * u[13] + diag[9] * u[14] + diag[13] * u[15]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
76       uik[14] = -(diag[2] * u[12] + diag[6] * u[13] + diag[10] * u[14] + diag[14] * u[15]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
77       uik[15] = -(diag[3] * u[12] + diag[7] * u[13] + diag[11] * u[14] + diag[15] * u[15]);  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
147     diag = ba + k * 16;  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
148     PetscCall(PetscArraycpy(diag, dk, 16));  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
150       PetscCall(PetscKernel_A_gets_inverse_A_4(diag, shift, allowzeropivot, &zeropivotdetected));  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
153       PetscCall(PetscKernel_A_gets_inverse_A_4_nopivot(diag));  in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()