Lines Matching refs:b
18 Mat_SeqBAIJ *a = (Mat_SeqBAIJ *)A->data, *b = (Mat_SeqBAIJ *)C->data; in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering() local
20 … PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *ajtmp, *bjtmp, *bdiag = b… 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()
85 pv = b->a + bs2 * bi[i]; in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
86 pj = b->j + bi[i]; 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()
110 PetscCall(PetscLogFlops(1.333333333333 * bs * bs2 * b->mbs)); /* from inverting diagonal blocks */ in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
117 Mat_SeqBAIJ *a = (Mat_SeqBAIJ *)A->data, *b = (Mat_SeqBAIJ *)C->data; in MatLUFactorNumeric_SeqBAIJ_N() local
118 IS isrow = b->row, isicol = b->icol; in MatLUFactorNumeric_SeqBAIJ_N()
120 PetscInt i, j, k, n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j; in MatLUFactorNumeric_SeqBAIJ_N()
121 PetscInt *ajtmp, *bjtmp, nz, nzL, row, *bdiag = b->diag, *pj; 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()
181 pv = b->a + bs2 * bi[i]; in MatLUFactorNumeric_SeqBAIJ_N()
182 pj = b->j + bi[i]; 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()
242 PetscCall(PetscLogFlops(1.333333333333 * bs * bs2 * b->mbs)); /* from inverting diagonal blocks */ in MatLUFactorNumeric_SeqBAIJ_N()
254 Mat_SeqBAIJ *a = (Mat_SeqBAIJ *)A->data, *b; in MatILUFactorSymbolic_SeqBAIJ_ilu0() local
261 b = (Mat_SeqBAIJ *)fact->data; in MatILUFactorSymbolic_SeqBAIJ_ilu0()
264 PetscCall(PetscShmgetAllocateArray(bs2 * ai[n], sizeof(PetscScalar), (void **)&b->a)); in MatILUFactorSymbolic_SeqBAIJ_ilu0()
265 PetscCall(PetscShmgetAllocateArray(ai[n], sizeof(PetscInt), (void **)&b->j)); in MatILUFactorSymbolic_SeqBAIJ_ilu0()
266 PetscCall(PetscShmgetAllocateArray(n + 1, sizeof(PetscInt), (void **)&b->i)); in MatILUFactorSymbolic_SeqBAIJ_ilu0()
267 b->free_a = PETSC_TRUE; in MatILUFactorSymbolic_SeqBAIJ_ilu0()
268 b->free_ij = PETSC_TRUE; in MatILUFactorSymbolic_SeqBAIJ_ilu0()
271 if (!b->diag) PetscCall(PetscMalloc1(n + 1, &b->diag)); in MatILUFactorSymbolic_SeqBAIJ_ilu0()
272 bdiag = b->diag; in MatILUFactorSymbolic_SeqBAIJ_ilu0()
274 if (n > 0) PetscCall(PetscArrayzero(b->a, bs2 * ai[n])); in MatILUFactorSymbolic_SeqBAIJ_ilu0()
277 bi = b->i; in MatILUFactorSymbolic_SeqBAIJ_ilu0()
278 bj = b->j; in MatILUFactorSymbolic_SeqBAIJ_ilu0()
313 Mat_SeqBAIJ *a = (Mat_SeqBAIJ *)A->data, *b; in MatILUFactorSymbolic_SeqBAIJ() local
359 b = (Mat_SeqBAIJ *)fact->data; in MatILUFactorSymbolic_SeqBAIJ()
360 b->row = isrow; in MatILUFactorSymbolic_SeqBAIJ()
361 b->col = iscol; in MatILUFactorSymbolic_SeqBAIJ()
362 b->icol = isicol; in MatILUFactorSymbolic_SeqBAIJ()
365 b->pivotinblocks = (info->pivotinblocks) ? PETSC_TRUE : PETSC_FALSE; in MatILUFactorSymbolic_SeqBAIJ()
367 PetscCall(PetscMalloc1((n + 1) * bs, &b->solve_work)); in MatILUFactorSymbolic_SeqBAIJ()
482 b = (Mat_SeqBAIJ *)fact->data; in MatILUFactorSymbolic_SeqBAIJ()
483 b->free_ij = PETSC_TRUE; in MatILUFactorSymbolic_SeqBAIJ()
484 PetscCall(PetscShmgetAllocateArray(bs2 * (bdiag[0] + 1), sizeof(PetscScalar), (void **)&b->a)); in MatILUFactorSymbolic_SeqBAIJ()
485 b->free_a = PETSC_TRUE; in MatILUFactorSymbolic_SeqBAIJ()
487 b->j = bj; in MatILUFactorSymbolic_SeqBAIJ()
488 b->i = bi; in MatILUFactorSymbolic_SeqBAIJ()
489 b->diag = bdiag; in MatILUFactorSymbolic_SeqBAIJ()
490 b->ilen = NULL; in MatILUFactorSymbolic_SeqBAIJ()
491 b->imax = NULL; in MatILUFactorSymbolic_SeqBAIJ()
492 b->row = isrow; in MatILUFactorSymbolic_SeqBAIJ()
493 b->col = iscol; in MatILUFactorSymbolic_SeqBAIJ()
496 b->icol = isicol; in MatILUFactorSymbolic_SeqBAIJ()
498 PetscCall(PetscMalloc1(bs * n + bs, &b->solve_work)); in MatILUFactorSymbolic_SeqBAIJ()
501 b->maxnz = b->nz = bdiag[0] + 1; in MatILUFactorSymbolic_SeqBAIJ()