61   PetscInt          *bdiag, row, nnz, nzi, reallocs = 0, nzbd, *im;  in MatLUFactorSymbolic_SeqAIJ()  local
79   PetscCall(PetscMalloc1(n + 1, &bdiag));  in MatLUFactorSymbolic_SeqAIJ()
80   bi[0] = bdiag[0] = 0;  in MatLUFactorSymbolic_SeqAIJ()
105       nzbd  = bdiag[row] + 1;     /* num of entries in the row with column index <= row */  in MatLUFactorSymbolic_SeqAIJ()
122     bdiag[i] = nzbd; /* note: bdiag[i] = nnzL as input for PetscFreeSpaceContiguous_LU() */  in MatLUFactorSymbolic_SeqAIJ()
146   PetscCall(PetscFreeSpaceContiguous_LU(&free_space, bj, n, bi, bdiag));  in MatLUFactorSymbolic_SeqAIJ()
154   PetscCall(PetscShmgetAllocateArray(bdiag[0] + 1, sizeof(PetscScalar), (void **)&b->a));  in MatLUFactorSymbolic_SeqAIJ()
158   b->diag   = bdiag;  in MatLUFactorSymbolic_SeqAIJ()
169   b->maxnz = b->nz = bdiag[0] + 1;  in MatLUFactorSymbolic_SeqAIJ()
176     B->info.fill_ratio_needed = ((PetscReal)(bdiag[0] + 1)) / ((PetscReal)ai[n]);  in MatLUFactorSymbolic_SeqAIJ()
224   const PetscInt   n = A->rmap->n, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *bdiag = b->diag;  in MatLUFactorNumeric_SeqAIJ()  local
275       nz    = bdiag[i] - bdiag[i + 1];  in MatLUFactorNumeric_SeqAIJ()
276       bjtmp = bj + bdiag[i + 1] + 1;  in MatLUFactorNumeric_SeqAIJ()
294           pv         = ba + bdiag[row];  in MatLUFactorNumeric_SeqAIJ()
298           pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */  in MatLUFactorNumeric_SeqAIJ()
299           pv = ba + bdiag[row + 1] + 1;  in MatLUFactorNumeric_SeqAIJ()
300           nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries in U(row,:) excluding diag */  in MatLUFactorNumeric_SeqAIJ()
320       pv = ba + bdiag[i + 1] + 1;  in MatLUFactorNumeric_SeqAIJ()
321       pj = b->j + bdiag[i + 1] + 1;  in MatLUFactorNumeric_SeqAIJ()
322       nz = bdiag[i] - bdiag[i + 1] - 1;  in MatLUFactorNumeric_SeqAIJ()
335       pv  = ba + bdiag[i];  in MatLUFactorNumeric_SeqAIJ()
1477   PetscInt       i, j, k = 0, nz, *bi, *bj, *bdiag;  in MatILUFactorSymbolic_SeqAIJ_ilu0()  local
1495   bdiag = b->diag;  in MatILUFactorSymbolic_SeqAIJ_ilu0()
1511   bdiag[n] = bi[n] - 1;  in MatILUFactorSymbolic_SeqAIJ_ilu0()
1518     bdiag[i] = bdiag[i + 1] + nz + 1;  in MatILUFactorSymbolic_SeqAIJ_ilu0()
1545   PetscInt          *bdiag, prow, fm, nzbd, reallocs = 0, dcount = 0;  in MatILUFactorSymbolic_SeqAIJ()  local
1577   PetscCall(PetscMalloc1(n + 1, &bdiag));  in MatILUFactorSymbolic_SeqAIJ()
1578   bi[0] = bdiag[0] = 0;  in MatILUFactorSymbolic_SeqAIJ()
1617       nnz      = bdiag[prow];  in MatILUFactorSymbolic_SeqAIJ()
1626     bdiag[i]  = nzbd;  in MatILUFactorSymbolic_SeqAIJ()
1642 …PetscCheck(*(bj_ptr[i] + bdiag[i]) == i, PETSC_COMM_SELF, PETSC_ERR_MAT_LU_ZRPVT, "Row %" PetscInt…  in MatILUFactorSymbolic_SeqAIJ()
1656   PetscCall(PetscFreeSpaceContiguous_LU(&free_space, bj, n, bi, bdiag));  in MatILUFactorSymbolic_SeqAIJ()
1664     PetscReal af = ((PetscReal)(bdiag[0] + 1)) / ((PetscReal)ai[n]);  in MatILUFactorSymbolic_SeqAIJ()
1676   PetscCall(PetscShmgetAllocateArray(bdiag[0] + 1, sizeof(PetscScalar), (void **)&b->a));  in MatILUFactorSymbolic_SeqAIJ()
1680   b->diag   = bdiag;  in MatILUFactorSymbolic_SeqAIJ()
1692   b->maxnz = b->nz = bdiag[0] + 1;  in MatILUFactorSymbolic_SeqAIJ()
1696   fact->info.fill_ratio_needed = ((PetscReal)(bdiag[0] + 1)) / ((PetscReal)ai[n]);  in MatILUFactorSymbolic_SeqAIJ()
1710   PetscInt         i, j, mbs = A->rmap->n, *bi = b->i, *bj = b->j, *bdiag = b->diag, *bjtmp;  in MatCholeskyFactorNumeric_SeqAIJ()  local
1788         uikdi = -ba[ili] * ba[bdiag[i]]; /* diagonal(k) */  in MatCholeskyFactorNumeric_SeqAIJ()
1830       ba[bdiag[k]] = 1.0 / dk; /* U(k,k) */  in MatCholeskyFactorNumeric_SeqAIJ()