11   MatScalar      *ba = b->a, *aa, *ap, *dk, *uik;  in MatCholeskyFactorNumeric_SeqSBAIJ_4()  local
25   PetscCall(PetscMalloc2(16, &dk, 16, &uik));  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
50             dk[k1]          = aa[k * 16 + k1];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
52             aa[j * 16 + k1] = dk[k1];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
59           for (k = 0; k < 16; k++) dk[k] = ap[k]; /* dk <- j-th block of aa */  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
61             for (k1 = 0; k1 < 4; k1++) *ap++ = dk[k + 4 * k1];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
84     PetscCall(PetscArraycpy(dk, rtmp + k * 16, 16));  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
118       dk[0] += uik[0] * u[0] + uik[1] * u[1] + uik[2] * u[2] + uik[3] * u[3];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
119       dk[1] += uik[4] * u[0] + uik[5] * u[1] + uik[6] * u[2] + uik[7] * u[3];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
120       dk[2] += uik[8] * u[0] + uik[9] * u[1] + uik[10] * u[2] + uik[11] * u[3];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
121       dk[3] += uik[12] * u[0] + uik[13] * u[1] + uik[14] * u[2] + uik[15] * u[3];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
123       dk[4] += uik[0] * u[4] + uik[1] * u[5] + uik[2] * u[6] + uik[3] * u[7];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
124       dk[5] += uik[4] * u[4] + uik[5] * u[5] + uik[6] * u[6] + uik[7] * u[7];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
125       dk[6] += uik[8] * u[4] + uik[9] * u[5] + uik[10] * u[6] + uik[11] * u[7];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
126       dk[7] += uik[12] * u[4] + uik[13] * u[5] + uik[14] * u[6] + uik[15] * u[7];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
128       dk[8] += uik[0] * u[8] + uik[1] * u[9] + uik[2] * u[10] + uik[3] * u[11];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
129       dk[9] += uik[4] * u[8] + uik[5] * u[9] + uik[6] * u[10] + uik[7] * u[11];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
130       dk[10] += uik[8] * u[8] + uik[9] * u[9] + uik[10] * u[10] + uik[11] * u[11];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
131       dk[11] += uik[12] * u[8] + uik[13] * u[9] + uik[14] * u[10] + uik[15] * u[11];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
133       dk[12] += uik[0] * u[12] + uik[1] * u[13] + uik[2] * u[14] + uik[3] * u[15];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
134       dk[13] += uik[4] * u[12] + uik[5] * u[13] + uik[6] * u[14] + uik[7] * u[15];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
135       dk[14] += uik[8] * u[12] + uik[9] * u[13] + uik[10] * u[14] + uik[11] * u[15];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
136       dk[15] += uik[12] * u[12] + uik[13] * u[13] + uik[14] * u[14] + uik[15] * u[15];  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
186     PetscCall(PetscArraycpy(diag, dk, 16));  in MatCholeskyFactorNumeric_SeqSBAIJ_4()
218   PetscCall(PetscFree2(dk, uik));  in MatCholeskyFactorNumeric_SeqSBAIJ_4()