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
44     nz    = bdiag[i] - bdiag[i + 1];  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
45     bjtmp = bj + bdiag[i + 1] + 1;  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
67         pv = b->a + bs2 * bdiag[row];  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
70         pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
71         pv = b->a + bs2 * (bdiag[row + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
72         nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries inU(row,:), excluding diag */  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
91     pv = b->a + bs2 * bdiag[i];  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
92     pj = b->j + bdiag[i];  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
98     pv = b->a + bs2 * (bdiag[i + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
99     pj = b->j + bdiag[i + 1] + 1;  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
100     nz = bdiag[i] - bdiag[i + 1] - 1;  in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
121   PetscInt       *ajtmp, *bjtmp, nz, nzL, row, *bdiag = b->diag, *pj;  in MatLUFactorNumeric_SeqBAIJ_N()  local
146     nz    = bdiag[i] - bdiag[i + 1];  in MatLUFactorNumeric_SeqBAIJ_N()
147     bjtmp = bj + bdiag[i + 1] + 1;  in MatLUFactorNumeric_SeqBAIJ_N()
169         pv = b->a + bs2 * bdiag[row];  in MatLUFactorNumeric_SeqBAIJ_N()
171         pj = b->j + bdiag[row + 1] + 1;                  /* beginning of U(row,:) */  in MatLUFactorNumeric_SeqBAIJ_N()
172         pv = b->a + bs2 * (bdiag[row + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_N()
173         nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries inU(row,:), excluding diag */  in MatLUFactorNumeric_SeqBAIJ_N()
187     pv = b->a + bs2 * bdiag[i];  in MatLUFactorNumeric_SeqBAIJ_N()
188     pj = b->j + bdiag[i];  in MatLUFactorNumeric_SeqBAIJ_N()
195     pv = b->a + bs2 * (bdiag[i + 1] + 1);  in MatLUFactorNumeric_SeqBAIJ_N()
196     pj = b->j + bdiag[i + 1] + 1;  in MatLUFactorNumeric_SeqBAIJ_N()
197     nz = bdiag[i] - bdiag[i + 1] - 1;  in MatLUFactorNumeric_SeqBAIJ_N()
256   PetscInt       i, j, nz, *bi, *bj, *bdiag, bi_temp;  in MatILUFactorSymbolic_SeqBAIJ_ilu0()  local
272   bdiag = b->diag;  in MatILUFactorSymbolic_SeqBAIJ_ilu0()
294   bdiag[n] = bi[n] - 1;  in MatILUFactorSymbolic_SeqBAIJ_ilu0()
306     bdiag[i] = bi_temp - 1;  in MatILUFactorSymbolic_SeqBAIJ_ilu0()
318   PetscInt          *bdiag, prow, fm, nzbd, reallocs = 0, dcount = 0;  in MatILUFactorSymbolic_SeqBAIJ()  local
378   PetscCall(PetscMalloc1(n + 1, &bdiag));  in MatILUFactorSymbolic_SeqBAIJ()
379   bdiag[0] = 0;  in MatILUFactorSymbolic_SeqBAIJ()
418       nnz      = bdiag[prow];  in MatILUFactorSymbolic_SeqBAIJ()
428     bdiag[i]  = nzbd;  in MatILUFactorSymbolic_SeqBAIJ()
446 …PetscCheck(*(bj_ptr[i] + bdiag[i]) == i, PETSC_COMM_SELF, PETSC_ERR_MAT_LU_ZRPVT, "Row %" PetscInt…  in MatILUFactorSymbolic_SeqBAIJ()
462   PetscCall(PetscFreeSpaceContiguous_LU(&free_space, bj, n, bi, bdiag));  in MatILUFactorSymbolic_SeqBAIJ()
470     PetscReal af = ((PetscReal)(bdiag[0] + 1)) / ((PetscReal)ai[n]);  in MatILUFactorSymbolic_SeqBAIJ()
484   PetscCall(PetscShmgetAllocateArray(bs2 * (bdiag[0] + 1), sizeof(PetscScalar), (void **)&b->a));  in MatILUFactorSymbolic_SeqBAIJ()
489   b->diag = bdiag;  in MatILUFactorSymbolic_SeqBAIJ()
501   b->maxnz = b->nz = bdiag[0] + 1;  in MatILUFactorSymbolic_SeqBAIJ()
505   fact->info.fill_ratio_needed = ((PetscReal)(bdiag[0] + 1)) / ((PetscReal)ai[n]);  in MatILUFactorSymbolic_SeqBAIJ()