| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact2.c | 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()
[all …]
|
| H A D | baijfact13.c | 175 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_3() local
198 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqBAIJ_3()
199 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_3()
221 pv = b->a + bs2 * bdiag[row]; in MatLUFactorNumeric_SeqBAIJ_3()
225 pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */ in MatLUFactorNumeric_SeqBAIJ_3()
226 pv = b->a + bs2 * (bdiag[row + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_3()
227 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries in U(row,:) excluding diag */ in MatLUFactorNumeric_SeqBAIJ_3()
247 pv = b->a + bs2 * bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_3()
248 pj = b->j + bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_3()
254 pj = b->j + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_3()
[all …]
|
| H A D | baijfact11.c | 237 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_4() local
266 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqBAIJ_4()
267 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_4()
289 pv = b->a + bs2 * bdiag[row]; in MatLUFactorNumeric_SeqBAIJ_4()
293 pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */ in MatLUFactorNumeric_SeqBAIJ_4()
294 pv = b->a + bs2 * (bdiag[row + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_4()
295 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries inU(row,:), excluding diag */ in MatLUFactorNumeric_SeqBAIJ_4()
315 pv = b->a + bs2 * bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_4()
316 pj = b->j + bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_4()
322 pv = b->a + bs2 * (bdiag[i + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_4()
[all …]
|
| H A D | baijfact3.c | 159 PetscInt *bdiag, row, nnz, nzi, reallocs = 0, nzbd, *im; in MatLUFactorSymbolic_SeqBAIJ() local
181 PetscCall(PetscMalloc1(n + 1, &bdiag)); in MatLUFactorSymbolic_SeqBAIJ()
182 bi[0] = bdiag[0] = 0; in MatLUFactorSymbolic_SeqBAIJ()
207 nzbd = bdiag[row] + 1; /* num of entries in the row with column index <= row */ in MatLUFactorSymbolic_SeqBAIJ()
224 bdiag[i] = nzbd; /* note : bdaig[i] = nnzL as input for PetscFreeSpaceContiguous_LU() */ in MatLUFactorSymbolic_SeqBAIJ()
247 PetscCall(PetscFreeSpaceContiguous_LU(&free_space, bj, n, bi, bdiag)); in MatLUFactorSymbolic_SeqBAIJ()
256 PetscCall(PetscShmgetAllocateArray((bdiag[0] + 1) * bs2, sizeof(PetscScalar), (void **)&b->a)); in MatLUFactorSymbolic_SeqBAIJ()
260 b->diag = bdiag; in MatLUFactorSymbolic_SeqBAIJ()
272 b->maxnz = b->nz = bdiag[0] + 1; in MatLUFactorSymbolic_SeqBAIJ()
279 B->info.fill_ratio_needed = ((PetscReal)(bdiag[0] + 1)) / ((PetscReal)ai[n]); in MatLUFactorSymbolic_SeqBAIJ()
|
| H A D | baijfact9.c | 315 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_5() local
338 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqBAIJ_5()
339 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_5()
361 pv = b->a + bs2 * bdiag[row]; in MatLUFactorNumeric_SeqBAIJ_5()
365 pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */ in MatLUFactorNumeric_SeqBAIJ_5()
366 pv = b->a + bs2 * (bdiag[row + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_5()
367 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries inU(row,:), excluding diag */ in MatLUFactorNumeric_SeqBAIJ_5()
387 pv = b->a + bs2 * bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_5()
388 pj = b->j + bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_5()
394 pv = b->a + bs2 * (bdiag[i + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_5()
[all …]
|
| H A D | baijfact.c | 14 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag; in MatLUFactorNumeric_SeqBAIJ_2() local
37 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqBAIJ_2()
38 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_2()
62 pv = b->a + bs2 * bdiag[row]; in MatLUFactorNumeric_SeqBAIJ_2()
66 pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */ in MatLUFactorNumeric_SeqBAIJ_2()
67 pv = b->a + bs2 * (bdiag[row + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_2()
68 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries inU(row,:), excluding diag */ in MatLUFactorNumeric_SeqBAIJ_2()
88 pv = b->a + bs2 * bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_2()
89 pj = b->j + bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_2()
99 pv = b->a + bs2 * (bdiag[i + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_2()
[all …]
|
| H A D | baijfact7.c | 373 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_6() local
396 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqBAIJ_6()
397 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_6()
419 pv = b->a + bs2 * bdiag[row]; in MatLUFactorNumeric_SeqBAIJ_6()
423 pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */ in MatLUFactorNumeric_SeqBAIJ_6()
424 pv = b->a + bs2 * (bdiag[row + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_6()
425 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries inU(row,:), excluding diag */ in MatLUFactorNumeric_SeqBAIJ_6()
445 pv = b->a + bs2 * bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_6()
446 pj = b->j + bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_6()
452 pv = b->a + bs2 * (bdiag[i + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_6()
[all …]
|
| H A D | baijfact5.c | 469 const PetscInt *ajtmp, *bjtmp, *bdiag = b->diag, *pj, bs2 = a->bs2; in MatLUFactorNumeric_SeqBAIJ_7() local
492 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqBAIJ_7()
493 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqBAIJ_7()
515 pv = b->a + bs2 * bdiag[row]; in MatLUFactorNumeric_SeqBAIJ_7()
519 pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */ in MatLUFactorNumeric_SeqBAIJ_7()
520 pv = b->a + bs2 * (bdiag[row + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_7()
521 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries inU(row,:), excluding diag */ in MatLUFactorNumeric_SeqBAIJ_7()
541 pv = b->a + bs2 * bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_7()
542 pj = b->j + bdiag[i]; in MatLUFactorNumeric_SeqBAIJ_7()
548 pv = b->a + bs2 * (bdiag[i + 1] + 1); in MatLUFactorNumeric_SeqBAIJ_7()
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | inode.c | 1204 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
1301 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqAIJ_Inode()
1302 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1321 pv = b->a + bdiag[row]; in MatLUFactorNumeric_SeqAIJ_Inode()
1324 pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */ in MatLUFactorNumeric_SeqAIJ_Inode()
1325 pv = b->a + bdiag[row + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1326 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries in U(row,:) excluding diag */ in MatLUFactorNumeric_SeqAIJ_Inode()
1345 pv = b->a + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1346 pj = b->j + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1347 nz = bdiag[i] - bdiag[i + 1] - 1; in MatLUFactorNumeric_SeqAIJ_Inode()
[all …]
|
| H A D | aijfact.c | 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()
[all …]
|
| H A D | inode2.c | 40 PetscCall(PetscFree3(a->inode.ibdiag, a->inode.bdiag, a->inode.ssor_work)); in MatDestroy_SeqAIJ_Inode()
65 b->inode.bdiag = NULL; in MatCreate_SeqAIJ_Inode()
|
| H A D | aij.h | 126 MatScalar *bdiag, *ibdiag, *ssor_work; /* diagonal blocks of matrices */ member
|
| /petsc/src/mat/utils/ |
| H A D | freespace.c | 60 …Contiguous_LU(PetscFreeSpaceList *head, PetscInt *space, PetscInt n, PetscInt *bi, PetscInt *bdiag) in PetscFreeSpaceContiguous_LU() argument
70 nnzL = bdiag[0]; in PetscFreeSpaceContiguous_LU()
87 nnzL = bdiag[row]; in PetscFreeSpaceContiguous_LU()
92 bdiag[row] = bi_temp - 1; in PetscFreeSpaceContiguous_LU()
93 space[bdiag[row]] = row; in PetscFreeSpaceContiguous_LU()
112 bdiag[n] = bdiag[n - 1] - 1; in PetscFreeSpaceContiguous_LU()
|
| /petsc/include/petsc/private/ |
| H A D | matisimpl.h | 22 PetscScalar *bdiag; /* Used by MatInvertBlockDiagonal_IS */ member
|
| /petsc/doc/changes/ |
| H A D | 300.md | 74 - The bdiag matrix formats have been removed
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 1350 PetscInt i, j, mbs = A->rmap->n, *bi = b->i, *bj = b->j, *bdiag = b->diag, *bjtmp; in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering() local
1427 uikdi = -ba[ili] * ba[bdiag[i]]; /* diagonal(k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
1468 ba[bdiag[k]] = 1.0 / dk; /* U(k,k) */ in MatCholeskyFactorNumeric_SeqSBAIJ_1_NaturalOrdering()
|
| /petsc/src/mat/impls/is/ |
| H A D | matis.c | 2108 PetscCall(PetscFree(b->bdiag)); in MatDestroy_IS()
2530 if (!is->bdiag) PetscCall(PetscMalloc1(bs * mat->rmap->n, &is->bdiag)); in MatInvertBlockDiagonal_IS()
2534 PetscCall(PetscSFReduceBegin(is->sf, nodeType, lv, is->bdiag, MPI_REPLACE)); in MatInvertBlockDiagonal_IS()
2535 PetscCall(PetscSFReduceEnd(is->sf, nodeType, lv, is->bdiag, MPI_REPLACE)); in MatInvertBlockDiagonal_IS()
2537 if (values) *values = is->bdiag; in MatInvertBlockDiagonal_IS()
2728 PetscCall(PetscFree(is->bdiag)); in MatSetLocalToGlobalMapping_IS()
|