237   const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2;  in MatLUFactorNumeric_SeqBAIJ_4()  local
255   PetscCall(PetscMalloc2(bs2 * n, &rtmp, bs2, &mwork));  in MatLUFactorNumeric_SeqBAIJ_4()
256   PetscCall(PetscArrayzero(rtmp, bs2 * n));  in MatLUFactorNumeric_SeqBAIJ_4()
263     for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_4()
268     for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_4()
273     v     = aa + bs2 * ai[r[i]];  in MatLUFactorNumeric_SeqBAIJ_4()
274     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, bs2));  in MatLUFactorNumeric_SeqBAIJ_4()
281       pc  = rtmp + bs2 * row;  in MatLUFactorNumeric_SeqBAIJ_4()
282       for (flg = 0, j = 0; j < bs2; j++) {  in MatLUFactorNumeric_SeqBAIJ_4()
289         pv = b->a + bs2 * bdiag[row];  in MatLUFactorNumeric_SeqBAIJ_4()
294         pv = b->a + bs2 * (bdiag[row + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_4()
299           v = rtmp + bs2 * pj[j];  in MatLUFactorNumeric_SeqBAIJ_4()
301           pv += bs2;  in MatLUFactorNumeric_SeqBAIJ_4()
309     pv = b->a + bs2 * bi[i];  in MatLUFactorNumeric_SeqBAIJ_4()
312     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_4()
315     pv = b->a + bs2 * bdiag[i];  in MatLUFactorNumeric_SeqBAIJ_4()
317     PetscCall(PetscArraycpy(pv, rtmp + bs2 * pj[0], bs2));  in MatLUFactorNumeric_SeqBAIJ_4()
322     pv = b->a + bs2 * (bdiag[i + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_4()
325     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_4()
557   const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2;  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()  local
567   PetscCall(PetscMalloc2(bs2 * n, &rtmp, bs2, &mwork));  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
568   PetscCall(PetscArrayzero(rtmp, bs2 * n));  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
581     for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
586     for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
591     v     = aa + bs2 * ai[i];  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
592     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ajtmp[j], v + bs2 * j, bs2));  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
599       pc  = rtmp + bs2 * row;  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
600       for (flg = 0, j = 0; j < bs2; j++) {  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
607         pv = b->a + bs2 * bdiag[row];  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
612         pv = b->a + bs2 * (bdiag[row + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
617           v = rtmp + bs2 * pj[j];  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
619           pv += bs2;  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
627     pv = b->a + bs2 * bi[i];  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
630     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
633     pv = b->a + bs2 * bdiag[i];  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
635     PetscCall(PetscArraycpy(pv, rtmp + bs2 * pj[0], bs2));  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
640     pv = b->a + bs2 * (bdiag[i + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()
643     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_4_NaturalOrdering()