175   const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2;  in MatLUFactorNumeric_SeqBAIJ_3()  local
187   PetscCall(PetscMalloc2(bs2 * n, &rtmp, bs2, &mwork));  in MatLUFactorNumeric_SeqBAIJ_3()
188   PetscCall(PetscArrayzero(rtmp, bs2 * n));  in MatLUFactorNumeric_SeqBAIJ_3()
195     for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_3()
200     for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_3()
205     v     = aa + bs2 * ai[r[i]];  in MatLUFactorNumeric_SeqBAIJ_3()
206     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmp[j]], v + bs2 * j, bs2));  in MatLUFactorNumeric_SeqBAIJ_3()
213       pc  = rtmp + bs2 * row;  in MatLUFactorNumeric_SeqBAIJ_3()
214       for (flg = 0, j = 0; j < bs2; j++) {  in MatLUFactorNumeric_SeqBAIJ_3()
221         pv = b->a + bs2 * bdiag[row];  in MatLUFactorNumeric_SeqBAIJ_3()
226         pv = b->a + bs2 * (bdiag[row + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_3()
231           v = rtmp + bs2 * pj[j];  in MatLUFactorNumeric_SeqBAIJ_3()
233           pv += bs2;  in MatLUFactorNumeric_SeqBAIJ_3()
241     pv = b->a + bs2 * bi[i];  in MatLUFactorNumeric_SeqBAIJ_3()
244     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_3()
247     pv = b->a + bs2 * bdiag[i];  in MatLUFactorNumeric_SeqBAIJ_3()
249     PetscCall(PetscArraycpy(pv, rtmp + bs2 * pj[0], bs2));  in MatLUFactorNumeric_SeqBAIJ_3()
255     pv = b->a + bs2 * (bdiag[i + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_3()
257     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_3()
427   const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2;  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()  local
437   PetscCall(PetscMalloc2(bs2 * n, &rtmp, bs2, &mwork));  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
438   PetscCall(PetscArrayzero(rtmp, bs2 * n));  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
445     for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
450     for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * bjtmp[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
455     v     = aa + bs2 * ai[i];  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
456     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ajtmp[j], v + bs2 * j, bs2));  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
463       pc  = rtmp + bs2 * row;  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
464       for (flg = 0, j = 0; j < bs2; j++) {  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
471         pv = b->a + bs2 * bdiag[row];  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
476         pv = b->a + bs2 * (bdiag[row + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
481           v = rtmp + bs2 * pj[j];  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
483           pv += bs2;  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
491     pv = b->a + bs2 * bi[i];  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
494     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
497     pv = b->a + bs2 * bdiag[i];  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
499     PetscCall(PetscArraycpy(pv, rtmp + bs2 * pj[0], bs2));  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
504     pv = b->a + bs2 * (bdiag[i + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()
507     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_3_NaturalOrdering()