20 … n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *ajtmp, *bjtmp, *bdiag = b->diag, *pj;  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()  local
70         pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
74           vv = rtmp + bs2 * pj[j];  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
86     pj = b->j + bi[i];  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
88     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
92     pj = b->j + bdiag[i];  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
93     PetscCall(PetscArraycpy(pv, rtmp + bs2 * pj[0], bs2));  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
99     pj = b->j + bdiag[i + 1] + 1;  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
101     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
121   PetscInt       *ajtmp, *bjtmp, nz, nzL, row, *bdiag = b->diag, *pj;  in MatLUFactorNumeric_SeqBAIJ_N()  local
171         pj = b->j + bdiag[row + 1] + 1;                  /* beginning of U(row,:) */  in MatLUFactorNumeric_SeqBAIJ_N()
174 …for (j = 0; j < nz; j++) PetscKernel_A_gets_A_minus_B_times_C(bs, rtmp + bs2 * pj[j], pc, pv + bs2…  in MatLUFactorNumeric_SeqBAIJ_N()
182     pj = b->j + bi[i];  in MatLUFactorNumeric_SeqBAIJ_N()
184     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_N()
188     pj = b->j + bdiag[i];  in MatLUFactorNumeric_SeqBAIJ_N()
189     PetscCall(PetscArraycpy(pv, rtmp + bs2 * pj[0], bs2));  in MatLUFactorNumeric_SeqBAIJ_N()
196     pj = b->j + bdiag[i + 1] + 1;  in MatLUFactorNumeric_SeqBAIJ_N()
198     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatLUFactorNumeric_SeqBAIJ_N()