22   MatScalar       *rtmp, *pc, *mwork, *pv, *vv, work[225];  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()  local
67         pv = b->a + bs2 * bdiag[row];  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
68         PetscKernel_A_gets_A_times_B(bs, pc, pv, mwork);  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
71         pv = b->a + bs2 * (bdiag[row + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
75           PetscKernel_A_gets_A_minus_B_times_C(bs, vv, pc, pv);  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
77           pv += bs2;  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
85     pv = b->a + bs2 * 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()
91     pv = b->a + bs2 * bdiag[i];  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
93     PetscCall(PetscArraycpy(pv, rtmp + bs2 * pj[0], bs2));  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
94 …PetscCall(PetscKernel_A_gets_inverse_A_15(pv, ipvt, work, info->shiftamount, allowzeropivot, &zero…  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
98     pv = b->a + bs2 * (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()
122   MatScalar      *rtmp, *pc, *mwork, *v, *pv, *aa     = a->a;  in MatLUFactorNumeric_SeqBAIJ_N()  local
169         pv = b->a + bs2 * bdiag[row];  in MatLUFactorNumeric_SeqBAIJ_N()
170         PetscKernel_A_gets_A_times_B(bs, pc, pv, mwork); /* *pc = *pc * (*pv); */  in MatLUFactorNumeric_SeqBAIJ_N()
172         pv = b->a + bs2 * (bdiag[row + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_N()
174 …j = 0; j < nz; j++) PetscKernel_A_gets_A_minus_B_times_C(bs, rtmp + bs2 * pj[j], pc, pv + bs2 * j);  in MatLUFactorNumeric_SeqBAIJ_N()
181     pv = b->a + bs2 * 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()
187     pv = b->a + bs2 * bdiag[i];  in MatLUFactorNumeric_SeqBAIJ_N()
189     PetscCall(PetscArraycpy(pv, rtmp + bs2 * pj[0], bs2));  in MatLUFactorNumeric_SeqBAIJ_N()
191 …PetscCall(PetscKernel_A_gets_inverse_A(bs, pv, v_pivots, v_work, allowzeropivot, &zeropivotdetecte…  in MatLUFactorNumeric_SeqBAIJ_N()
195     pv = b->a + bs2 * (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()