15   PetscInt        diag, bs = A->rmap->bs, bs2 = a->bs2, *pj, *v_pivots;  in MatILUFactorNumeric_SeqBAIJ_N_inplace()  local
28   PetscCall(PetscCalloc1(bs2 * (n + 1), &rtmp));  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
30   PetscCall(PetscMalloc3(bs, &v_work, bs2, &multiplier, bs, &v_pivots));  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
35     for (j = 0; j < nz; j++) PetscCall(PetscArrayzero(rtmp + bs2 * ajtmp[j], bs2));  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
39     v        = aa + bs2 * ai[r[i]];  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
40 … for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(rtmp + bs2 * ic[ajtmpold[j]], v + bs2 * j, bs2));  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
43       pc = rtmp + bs2 * row;  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
45       for (flg = 0, k = 0; k < bs2; k++) {  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
52         pv = ba + bs2 * diag_offset[row];  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
56         pv += bs2;  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
57 …j = 0; j < nz; j++) PetscKernel_A_gets_A_minus_B_times_C(bs, rtmp + bs2 * pj[j], pc, pv + bs2 * j);  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
58         PetscCall(PetscLogFlops(2.0 * bs * bs2 * (nz + 1.0) - bs));  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
63     pv = ba + bs2 * bi[i];  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
66     for (j = 0; j < nz; j++) PetscCall(PetscArraycpy(pv + bs2 * j, rtmp + bs2 * pj[j], bs2));  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
69     w = pv + bs2 * diag;  in MatILUFactorNumeric_SeqBAIJ_N_inplace()
84   PetscCall(PetscLogFlops(1.333333333333 * bs * bs2 * b->mbs)); /* from inverting diagonal blocks */  in MatILUFactorNumeric_SeqBAIJ_N_inplace()