11   MatScalar      *ba = b->a, *aa, *ap, *dk, *uik;  in MatCholeskyFactorNumeric_SeqSBAIJ_5()  local
24   PetscCall(PetscMalloc2(25, &dk, 25, &uik));  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
49             dk[k1]          = aa[k * 25 + k1];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
51             aa[j * 25 + k1] = dk[k1];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
58           for (k = 0; k < 25; k++) dk[k] = ap[k]; /* dk <- j-th block of aa */  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
60             for (k1 = 0; k1 < 5; k1++) *ap++ = dk[k + 5 * k1];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
83     PetscCall(PetscArraycpy(dk, rtmp + k * 25, 25));  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
127       dk[0] += uik[0] * u[0] + uik[1] * u[1] + uik[2] * u[2] + uik[3] * u[3] + uik[4] * u[4];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
128       dk[1] += uik[5] * u[0] + uik[6] * u[1] + uik[7] * u[2] + uik[8] * u[3] + uik[9] * u[4];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
129       dk[2] += uik[10] * u[0] + uik[11] * u[1] + uik[12] * u[2] + uik[13] * u[3] + uik[14] * u[4];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
130       dk[3] += uik[15] * u[0] + uik[16] * u[1] + uik[17] * u[2] + uik[18] * u[3] + uik[19] * u[4];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
131       dk[4] += uik[20] * u[0] + uik[21] * u[1] + uik[22] * u[2] + uik[23] * u[3] + uik[24] * u[4];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
133       dk[5] += uik[0] * u[5] + uik[1] * u[6] + uik[2] * u[7] + uik[3] * u[8] + uik[4] * u[9];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
134       dk[6] += uik[5] * u[5] + uik[6] * u[6] + uik[7] * u[7] + uik[8] * u[8] + uik[9] * u[9];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
135       dk[7] += uik[10] * u[5] + uik[11] * u[6] + uik[12] * u[7] + uik[13] * u[8] + uik[14] * u[9];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
136       dk[8] += uik[15] * u[5] + uik[16] * u[6] + uik[17] * u[7] + uik[18] * u[8] + uik[19] * u[9];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
137       dk[9] += uik[20] * u[5] + uik[21] * u[6] + uik[22] * u[7] + uik[23] * u[8] + uik[24] * u[9];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
139       dk[10] += uik[0] * u[10] + uik[1] * u[11] + uik[2] * u[12] + uik[3] * u[13] + uik[4] * u[14];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
140       dk[11] += uik[5] * u[10] + uik[6] * u[11] + uik[7] * u[12] + uik[8] * u[13] + uik[9] * u[14];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
141 … dk[12] += uik[10] * u[10] + uik[11] * u[11] + uik[12] * u[12] + uik[13] * u[13] + uik[14] * u[14];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
142 … dk[13] += uik[15] * u[10] + uik[16] * u[11] + uik[17] * u[12] + uik[18] * u[13] + uik[19] * u[14];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
143 … dk[14] += uik[20] * u[10] + uik[21] * u[11] + uik[22] * u[12] + uik[23] * u[13] + uik[24] * u[14];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
145       dk[15] += uik[0] * u[15] + uik[1] * u[16] + uik[2] * u[17] + uik[3] * u[18] + uik[4] * u[19];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
146       dk[16] += uik[5] * u[15] + uik[6] * u[16] + uik[7] * u[17] + uik[8] * u[18] + uik[9] * u[19];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
147 … dk[17] += uik[10] * u[15] + uik[11] * u[16] + uik[12] * u[17] + uik[13] * u[18] + uik[14] * u[19];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
148 … dk[18] += uik[15] * u[15] + uik[16] * u[16] + uik[17] * u[17] + uik[18] * u[18] + uik[19] * u[19];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
149 … dk[19] += uik[20] * u[15] + uik[21] * u[16] + uik[22] * u[17] + uik[23] * u[18] + uik[24] * u[19];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
151       dk[20] += uik[0] * u[20] + uik[1] * u[21] + uik[2] * u[22] + uik[3] * u[23] + uik[4] * u[24];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
152       dk[21] += uik[5] * u[20] + uik[6] * u[21] + uik[7] * u[22] + uik[8] * u[23] + uik[9] * u[24];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
153 … dk[22] += uik[10] * u[20] + uik[11] * u[21] + uik[12] * u[22] + uik[13] * u[23] + uik[14] * u[24];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
154 … dk[23] += uik[15] * u[20] + uik[16] * u[21] + uik[17] * u[22] + uik[18] * u[23] + uik[19] * u[24];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
155 … dk[24] += uik[20] * u[20] + uik[21] * u[21] + uik[22] * u[22] + uik[23] * u[23] + uik[24] * u[24];  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
215     PetscCall(PetscArraycpy(d, dk, 25));  in MatCholeskyFactorNumeric_SeqSBAIJ_5()
242   PetscCall(PetscFree2(dk, uik));  in MatCholeskyFactorNumeric_SeqSBAIJ_5()