Home
last modified time | relevance | path

Searched refs:bi (Results 1 – 25 of 49) sorted by relevance

12

/petsc/src/mat/utils/
H A Dfreespace.c60 …Contiguous_LU(PetscFreeSpaceList *head, PetscInt *space, PetscInt n, PetscInt *bi, PetscInt *bdiag) in PetscFreeSpaceContiguous_LU() argument
67 bi_temp = bi[n]; in PetscFreeSpaceContiguous_LU()
76 if (bi[row + 1] > total) break; in PetscFreeSpaceContiguous_LU()
78 nnz = bi[row + 1] - bi[row]; in PetscFreeSpaceContiguous_LU()
81 bi[row] = 0; in PetscFreeSpaceContiguous_LU()
83 bi[row] = bi[row - 1] + nnzL; /* nnzL of previous row */ in PetscFreeSpaceContiguous_LU()
88 bj = space + bi[row]; in PetscFreeSpaceContiguous_LU()
111 bi[n] = bi[n - 1] + nnzL; in PetscFreeSpaceContiguous_LU()
/petsc/src/mat/impls/baij/seq/
H A Dbaijfact2.c20 …const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *ajtmp, *bjtmp, *bdia… in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering() local
39 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
40 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
55 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
56 nzL = bi[i + 1] - bi[i]; 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()
87 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_15_NaturalOrdering()
120 PetscInt i, j, k, n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j; in MatLUFactorNumeric_SeqBAIJ_N() local
141 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_N()
[all …]
H A Dbaijfact13.c15 PetscInt i, j, n = a->mbs, *bi = b->i, *bj = b->j; in MatILUFactorNumeric_SeqBAIJ_3_inplace() local
37 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
38 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
96 nz = bi[row + 1] - diag_offset[row] - 1; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
127 pv = ba + 9 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
128 pj = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
129 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_3_inplace()
174 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j; in MatLUFactorNumeric_SeqBAIJ_3() local
193 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_3()
194 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_3()
[all …]
H A Dbaijfact4.c12 PetscInt i, j, n = a->mbs, *bi = b->i, *bj = b->j; in MatILUFactorNumeric_SeqBAIJ_N_inplace() local
33 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_N_inplace()
34 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_N_inplace()
55 nz = bi[row + 1] - diag_offset[row] - 1; in MatILUFactorNumeric_SeqBAIJ_N_inplace()
63 pv = ba + bs2 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_N_inplace()
64 pj = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_N_inplace()
65 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_N_inplace()
67 diag = diag_offset[i] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_N_inplace()
H A Dbaijfact.c13 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_2() local
32 PetscInt nzL, nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_2()
33 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_2()
48 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_2()
49 nzL = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_2()
82 pv = b->a + bs2 * bi[i]; in MatLUFactorNumeric_SeqBAIJ_2()
83 pj = b->j + bi[i]; in MatLUFactorNumeric_SeqBAIJ_2()
84 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_2()
122 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering() local
140 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_2_NaturalOrdering()
[all …]
H A Dbaijfact11.c15 PetscInt i, j, n = a->mbs, *bi = b->i, *bj = b->j; in MatILUFactorNumeric_SeqBAIJ_4_inplace() local
40 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
41 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
130 nz = bi[row + 1] - diag_offset[row] - 1; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
177 pv = ba + 16 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
178 pj = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
179 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_4_inplace()
236 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j; in MatLUFactorNumeric_SeqBAIJ_4() local
261 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_4()
262 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_4()
[all …]
H A Dbaijfact9.c15 PetscInt *bi = b->i, *bj = b->j, *ajtmpold, *ajtmp; in MatILUFactorNumeric_SeqBAIJ_5_inplace() local
44 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
45 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
180 nz = bi[row + 1] - diag_offset[row] - 1; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
246 pv = ba + 25 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
247 pj = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
248 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_5_inplace()
314 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j; in MatLUFactorNumeric_SeqBAIJ_5() local
333 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_5()
334 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_5()
[all …]
H A Dbaijfact3.c158 PetscInt *bi, *bj, *ajtmp; in MatLUFactorSymbolic_SeqBAIJ() local
180 PetscCall(PetscMalloc1(n + 1, &bi)); in MatLUFactorSymbolic_SeqBAIJ()
182 bi[0] = bdiag[0] = 0; in MatLUFactorSymbolic_SeqBAIJ()
213 bi[i + 1] = bi[i] + nzi; in MatLUFactorSymbolic_SeqBAIJ()
246 PetscCall(PetscMalloc1(bi[n], &bj)); in MatLUFactorSymbolic_SeqBAIJ()
247 PetscCall(PetscFreeSpaceContiguous_LU(&free_space, bj, n, bi, bdiag)); in MatLUFactorSymbolic_SeqBAIJ()
259 b->i = bi; in MatLUFactorSymbolic_SeqBAIJ()
H A Dbaijfact7.c14 …const PetscInt *ajtmpold, *ajtmp, *diag_offset = b->diag, *r, *ic, *bi = b->i, *bj = b->j, *ai = a… in MatILUFactorNumeric_SeqBAIJ_6_inplace() local
40 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
41 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
215 nz = bi[row + 1] - diag_offset[row] - 1; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
304 pv = ba + 36 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
305 pj = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
306 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_6_inplace()
372 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j; in MatLUFactorNumeric_SeqBAIJ_6() local
391 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_6()
392 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_6()
[all …]
H A Dbaijfact5.c13 const PetscInt *r, *ic, *bi = b->i, *bj = b->j, *ajtmp, *ai = a->i, *aj = a->j, *pj, *ajtmpold; in MatILUFactorNumeric_SeqBAIJ_7_inplace() local
43 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
44 ajtmp = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
272 nz = bi[row + 1] - diag_offset[row] - 1; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
387 pv = ba + 49 * bi[i]; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
388 pj = bj + bi[i]; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
389 nz = bi[i + 1] - bi[i]; in MatILUFactorNumeric_SeqBAIJ_7_inplace()
468 const PetscInt n = a->mbs, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j; in MatLUFactorNumeric_SeqBAIJ_7() local
487 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqBAIJ_7()
488 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqBAIJ_7()
[all …]
/petsc/src/tao/util/
H A Dtao_util.c252 PetscReal ai, bi, ci, di, ei; in MatDFischer() local
299 bi = PetscRealPart(u[i]) - PetscRealPart(x[i]); in MatDFischer()
300 ai = fischnorm(bi, PetscRealPart(f[i])); in MatDFischer()
303 da[i] = bi / ai - 1.0; in MatDFischer()
313 bi = PetscRealPart(x[i]) - PetscRealPart(l[i]); in MatDFischer()
314 ai = fischnorm(bi, PetscRealPart(f[i])); in MatDFischer()
317 da[i] = bi / ai - 1.0; in MatDFischer()
330 bi = PetscRealPart(x[i]) - PetscRealPart(u[i]); in MatDFischer()
331 ai = fischnorm(bi, PetscRealPart(f[i])); in MatDFischer()
334 ci = bi / ai + 1.0; in MatDFischer()
[all …]
/petsc/src/mat/impls/sbaij/seq/
H A Daijsbaij.c10 …PetscInt *ai = a->i, *aj = a->j, m = A->rmap->N, n = A->cmap->n, i, j, k, *bi, *bj, *rowlength… in MatConvert_SeqSBAIJ_SeqAIJ() local
49 bi = b->i; in MatConvert_SeqSBAIJ_SeqAIJ()
54 bi[0] = 0; in MatConvert_SeqSBAIJ_SeqAIJ()
61 bi[i + 1] = bi[i] + rowlengths[i * bs] / bs; in MatConvert_SeqSBAIJ_SeqAIJ()
63 …etscCheck(bi[mbs] == 2 * a->nz - diagcnt, PETSC_COMM_SELF, PETSC_ERR_PLIB, "bi[mbs]: %" PetscInt_F… in MatConvert_SeqSBAIJ_SeqAIJ()
168 …PetscInt *ai = a->i, *aj, m = A->rmap->N, n = A->cmap->N, i, j, *bi, *bj, *rowlengths, bs = … in MatConvert_SeqAIJ_SeqSBAIJ() local
201 bi = b->i; in MatConvert_SeqAIJ_SeqSBAIJ()
205 bi[0] = 0; in MatConvert_SeqAIJ_SeqSBAIJ()
217 bi[i + 1] = bi[i] + rowlengths[i]; in MatConvert_SeqAIJ_SeqSBAIJ()
240 …PetscInt *ai = a->i, *aj = a->j, m = A->rmap->N, n = A->cmap->n, *bi, *bj, *browlengths, nz, *… in MatConvert_SeqSBAIJ_SeqBAIJ() local
[all …]
H A Dsbaijfact.c562 const PetscInt *ai, *aj, *perm_ptr, mbs = a->mbs, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqSBAIJ_N() local
664 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
692 jmin = bi[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
693 jmax = bi[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_N()
736 PetscInt i, j, mbs = a->mbs, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering() local
796 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
824 jmin = bi[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
825 jmax = bi[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
870 PetscInt i, j, mbs = a->mbs, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqSBAIJ_2() local
978 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
[all …]
H A Dsbaijfact4.c10 PetscInt i, j, mbs = a->mbs, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() local
90 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
127 jmin = bi[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
128 jmax = bi[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
H A Dsbaijfact3.c9 const PetscInt *ai, *aj, *perm_ptr, mbs = a->mbs, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqSBAIJ_3() local
128 jmax = bi[i + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
165 jmin = bi[k]; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
166 jmax = bi[k + 1]; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
/petsc/src/mat/tests/
H A Dex303.c11 PetscInt *bi, *bj, nd; in main() local
55 …PetscCall(MatGetRowIJ(B, 0, PETSC_FALSE, PETSC_FALSE, &nd, (const PetscInt **)&bi, (const PetscInt… in main()
62 …for (int j = 0; j < bi[i + 1] - bi[i]; j++) PetscCall(MatSetValue(B2, i, garray[bj[bi[i] + j]], ba… in main()
67 …PetscCall(MatRestoreRowIJ(B, 0, PETSC_FALSE, PETSC_FALSE, &nd, (const PetscInt **)&bi, (const Pets… in main()
/petsc/src/mat/impls/aij/seq/
H A Dmatmatmult.c147 PetscInt *ai = a->i, *bi = b->i, *ci, *cj; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed() local
171 …PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, PetscIntSumTruncate(ai[am], bi[bm])), &… in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
181 bnzj = bi[brow + 1] - bi[brow]; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
182 bj = b->j + bi[brow]; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
229 afill = (PetscReal)ci[am] / (ai[am] + bi[bm]) + 1.e-5; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
252 PetscInt *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *bjj, *ci = c->i, *cj = c->j; in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted() local
288 bnzi = bi[brow + 1] - bi[brow]; in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted()
289 bjj = PetscSafePointerPlusOffset(bj, bi[brow]); in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted()
290 baj = PetscSafePointerPlusOffset(ba, bi[brow]); in MatMatMultNumeric_SeqAIJ_SeqAIJ_Sorted()
325 PetscInt *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *bjj, *ci = c->i, *cj = c->j; in MatMatMultNumeric_SeqAIJ_SeqAIJ_Scalable() local
[all …]
H A Daijfact.c60 PetscInt *bi, *bj; in MatLUFactorSymbolic_SeqAIJ() local
78 PetscCall(PetscShmgetAllocateArray(n + 1, sizeof(PetscInt), (void **)&bi)); in MatLUFactorSymbolic_SeqAIJ()
80 bi[0] = bdiag[0] = 0; in MatLUFactorSymbolic_SeqAIJ()
111 bi[i + 1] = bi[i] + nzi; in MatLUFactorSymbolic_SeqAIJ()
145 PetscCall(PetscShmgetAllocateArray(bi[n], sizeof(PetscInt), (void **)&bj)); in MatLUFactorSymbolic_SeqAIJ()
146 PetscCall(PetscFreeSpaceContiguous_LU(&free_space, bj, n, bi, bdiag)); in MatLUFactorSymbolic_SeqAIJ()
157 b->i = bi; 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
270 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqAIJ()
271 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqAIJ()
[all …]
H A Dinode.c1204 const PetscInt n = A->rmap->n, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *bdiag = b->diag; in MatLUFactorNumeric_SeqAIJ_Inode() local
1296 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
1297 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
1315 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
1317 nzL = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
1336 pv = b->a + bi[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
1337 pj = b->j + bi[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
1338 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
1367 nz = bi[i + 1] - bi[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
1368 bjtmp = bj + bi[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
[all …]
/petsc/src/mat/impls/aij/seq/bas/
H A Dbasfactor.c100 PetscInt mbs = A->rmap->n, *bi = b->i, *bj = b->j; in MatCholeskyFactorNumeric_SeqAIJ_Bas() local
112 …PetscCall(spbas_compress_pattern(bi, bj, mbs, mbs, SPBAS_DIAGONAL_OFFSETS, &Pattern, &mem_reductio… in MatCholeskyFactorNumeric_SeqAIJ_Bas()
113 PetscCall(PetscFree(bi)); in MatCholeskyFactorNumeric_SeqAIJ_Bas()
144 PetscCall(spbas_matrix_to_crs(matrix_L, &ba, &bi, &bj)); in MatCholeskyFactorNumeric_SeqAIJ_Bas()
145 b->i = bi; in MatCholeskyFactorNumeric_SeqAIJ_Bas()
/petsc/src/mat/impls/aij/mpi/mkl_cpardiso/
H A Dmkl_cpardiso.c145 const PetscInt *ai, *aj, *bi, *bj, *garray, m = A->rmap->n, *ajj, *bjj; in MatConvertToTriples_mpiaij_mpiaij_MKL_CPARDISO() local
158 bi = bb->i; in MatConvertToTriples_mpiaij_mpiaij_MKL_CPARDISO()
183 countB = bi[i + 1] - bi[i]; in MatConvertToTriples_mpiaij_mpiaij_MKL_CPARDISO()
185 bjj = bj + bi[i]; in MatConvertToTriples_mpiaij_mpiaij_MKL_CPARDISO()
216 …const PetscInt *ai, *aj, *bi, *bj, *garray, bs = A->rmap->bs, bs2 = bs * bs, m = A->rmap->n / b… in MatConvertToTriples_mpibaij_mpibaij_MKL_CPARDISO() local
229 bi = bb->i; in MatConvertToTriples_mpibaij_mpibaij_MKL_CPARDISO()
254 countB = bi[i + 1] - bi[i]; in MatConvertToTriples_mpibaij_mpibaij_MKL_CPARDISO()
256 bjj = bj + bi[i]; in MatConvertToTriples_mpibaij_mpibaij_MKL_CPARDISO()
289 …const PetscInt *ai, *aj, *bi, *bj, *garray, bs = A->rmap->bs, bs2 = bs * bs, m = A->rmap->n / b… in MatConvertToTriples_mpisbaij_mpisbaij_MKL_CPARDISO() local
301 bi = bb->i; in MatConvertToTriples_mpisbaij_mpisbaij_MKL_CPARDISO()
[all …]
/petsc/src/ksp/ksp/impls/cr/
H A Dcr.c14 PetscScalar ai, bi; in KSPSolve_CR() local
130 bi = btop / bbot; in KSPSolve_CR()
131 PetscCall(VecAYPX(P, bi, RT)); /* P <- RT + Bi P */ in KSPSolve_CR()
132 PetscCall(VecAYPX(AP, bi, ART)); /* AP <- ART + Bi AP */ in KSPSolve_CR()
/petsc/src/sys/tests/
H A Dex15.cxx53 PetscBLASInt bi; in main() local
136 bi = 0; in main()
183 (void)bi; in main()
/petsc/src/sys/utils/
H A Dsorti.c1163 PetscInt ai, bi, k; in PetscMergeMPIIntArray() local
1167 for (ai = 0, bi = 0, k = 0; ai < an || bi < bn;) { in PetscMergeMPIIntArray()
1169 for (; ai < an && (!bn || aI[ai] <= bI[bi]); ai++) (*L)[k++] = t = aI[ai]; in PetscMergeMPIIntArray()
1170 for (; bi < bn && bI[bi] == t; bi++); in PetscMergeMPIIntArray()
1171 for (; bi < bn && (!an || bI[bi] <= aI[ai]); bi++) (*L)[k++] = t = bI[bi]; in PetscMergeMPIIntArray()
/petsc/src/mat/impls/aij/mpi/
H A Dmpiov.c286 const PetscInt *gcols, *ai, *aj, *bi, *bj; in MatIncreaseOverlap_MPIAIJ_Send_Scalable() local
299 PetscCall(MatGetRowIJ(bmat, 0, PETSC_FALSE, PETSC_FALSE, &bn, &bi, &bj, &done)); in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
302 tnz = ai[an] + bi[bn]; in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
334 start = bi[row]; in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
335 end = bi[row + 1]; in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
374 start = bi[row]; in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
375 end = bi[row + 1]; in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
398 PetscCall(MatRestoreRowIJ(bmat, 0, PETSC_FALSE, PETSC_FALSE, &bn, &bi, &bj, &done)); in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
404 const PetscInt *gcols, *ai, *aj, *bi, *bj, *indices; in MatIncreaseOverlap_MPIAIJ_Local_Scalable() local
419 PetscCall(MatGetRowIJ(bmat, 0, PETSC_FALSE, PETSC_FALSE, &bn, &bi, &bj, &done)); in MatIncreaseOverlap_MPIAIJ_Local_Scalable()
[all …]

12