| /petsc/include/petsc/private/ |
| H A D | matimpl.h | 882 …reate(idx_start, lnk_max, nlnk, lnk, bt) ((PetscErrorCode)(PetscMalloc1(nlnk, &lnk) || PetscBTCrea… argument 884 …max, nlnk, lnk, bt, lnk_empty) ((PetscErrorCode)(PetscMalloc1(nlnk, &lnk) || PetscBTCreate(nlnk, &… argument 886 …try, PetscInt *PETSC_RESTRICT nlnk, PetscInt *PETSC_RESTRICT lnkdata, PetscInt *PETSC_RESTRICT lnk) in PetscLLInsertLocation_Private() argument 896 *lnkdata = lnk[location]; in PetscLLInsertLocation_Private() 899 lnk[location] = entry; in PetscLLInsertLocation_Private() 900 lnk[entry] = *lnkdata; in PetscLLInsertLocation_Private() 906 …etscInt idx_start, PetscInt *PETSC_RESTRICT nlnk, PetscInt *PETSC_RESTRICT lnk, PetscBT bt, PetscB… in PetscLLAdd_Private() argument 913 … PetscCall(PetscLLInsertLocation_Private(assume_sorted, k, idx_start, entry, nlnk, &lnkdata, lnk)); in PetscLLAdd_Private() 931 …dices, PetscInt idx_start, PetscInt *PETSC_RESTRICT nlnk, PetscInt *PETSC_RESTRICT lnk, PetscBT bt) in PetscLLAdd() argument 934 PetscCall(PetscLLAdd_Private(nidx, indices, idx_start, nlnk, lnk, bt, PETSC_FALSE)); in PetscLLAdd() [all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact3.c | 161 PetscInt nlnk, *lnk, k, **bi_ptr; in MatLUFactorSymbolic_SeqBAIJ() local 186 PetscCall(PetscLLCreate(n, n, nlnk, lnk, lnkbt)); in MatLUFactorSymbolic_SeqBAIJ() 201 PetscCall(PetscLLAddPerm(nnz, ajtmp, ic, n, &nlnk, lnk, lnkbt)); in MatLUFactorSymbolic_SeqBAIJ() 205 row = lnk[n]; in MatLUFactorSymbolic_SeqBAIJ() 209 PetscCall(PetscLLAddSortedLU(ajtmp, row, &nlnk, lnk, lnkbt, i, nzbd, im)); in MatLUFactorSymbolic_SeqBAIJ() 211 row = lnk[row]; in MatLUFactorSymbolic_SeqBAIJ() 219 k = lnk[n]; in MatLUFactorSymbolic_SeqBAIJ() 222 k = lnk[k]; in MatLUFactorSymbolic_SeqBAIJ() 234 PetscCall(PetscLLClean(n, n, nzi, lnk, current_space->array, lnkbt)); in MatLUFactorSymbolic_SeqBAIJ() 248 PetscCall(PetscLLDestroy(lnk, lnkbt)); in MatLUFactorSymbolic_SeqBAIJ()
|
| H A D | baijfact2.c | 322 PetscInt nlnk, *lnk, *lnk_lvl = NULL; in MatILUFactorSymbolic_SeqBAIJ() local 385 PetscCall(PetscIncompleteLLCreate(n, n, nlnk, lnk, lnk_lvl, lnkbt)); in MatILUFactorSymbolic_SeqBAIJ() 399 lnk[i] = -1; /* marker to indicate if diagonal exists */ in MatILUFactorSymbolic_SeqBAIJ() 400 PetscCall(PetscIncompleteLLInit(nnz, cols, n, ic, &nlnk, lnk, lnk_lvl, lnkbt)); in MatILUFactorSymbolic_SeqBAIJ() 404 if (diagonal_fill && lnk[i] == -1) { in MatILUFactorSymbolic_SeqBAIJ() 406 while (lnk[fm] < i) fm = lnk[fm]; in MatILUFactorSymbolic_SeqBAIJ() 407 lnk[i] = lnk[fm]; /* insert diagonal into linked list */ in MatILUFactorSymbolic_SeqBAIJ() 408 lnk[fm] = i; in MatILUFactorSymbolic_SeqBAIJ() 416 prow = lnk[n]; in MatILUFactorSymbolic_SeqBAIJ() 423 …PetscCall(PetscILULLAddSorted(nnz, cols, levels, cols_lvl, prow, &nlnk, lnk, lnk_lvl, lnkbt, prow)… in MatILUFactorSymbolic_SeqBAIJ() [all …]
|
| H A D | baijfact.c | 1025 …PetscInt nlnk, *lnk, *lnk_lvl = NULL, ncols, ncols_upper, *cols, *cols_lvl, *uj, **uj_pt… in MatICCFactorSymbolic_SeqBAIJ() local 1093 PetscCall(PetscIncompleteLLCreate(am, am, nlnk, lnk, lnk_lvl, lnkbt)); in MatICCFactorSymbolic_SeqBAIJ() 1117 …PetscCall(PetscIncompleteLLAdd(ncols_upper, cols, levels, cols_lvl, am, &nlnk, lnk, lnk_lvl, lnkbt… in MatICCFactorSymbolic_SeqBAIJ() 1133 …PetscCall(PetscIncompleteLLAddSorted(ncols, cols, levels, cols_lvl, am, &nlnk, lnk, lnk_lvl, lnkbt… in MatICCFactorSymbolic_SeqBAIJ() 1157 …PetscCall(PetscIncompleteLLClean(am, am, nzk, lnk, lnk_lvl, current_space->array, current_space_lv… in MatICCFactorSymbolic_SeqBAIJ() 1187 PetscCall(PetscIncompleteLLDestroy(lnk, lnkbt)); in MatICCFactorSymbolic_SeqBAIJ() 1255 PetscInt nlnk, *lnk, ncols, ncols_upper, *cols, *uj, **ui_ptr, *uj_ptr; in MatCholeskyFactorSymbolic_SeqBAIJ() local 1294 PetscCall(PetscLLCreate(mbs, mbs, nlnk, lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqBAIJ() 1313 PetscCall(PetscLLAdd(ncols_upper, cols, mbs, &nlnk, lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqBAIJ() 1326 PetscCall(PetscLLAddSorted(ncols, uj_ptr, mbs, &nlnk, lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqBAIJ() [all …]
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpimatmatmult.c | 202 PetscInt *lnk, i, pnz, row, *api, *apj, *Jptr, apnz, nspacedouble = 0, j, nzi; in MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable() local 242 PetscCall(PetscLLCondensedCreate(pN, pN, &lnk, &lnkbt)); in MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable() 257 PetscCall(PetscLLCondensedAddSorted(pnz, Jptr, lnk, lnkbt)); in MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable() 265 PetscCall(PetscLLCondensedAddSorted(pnz, Jptr, lnk, lnkbt)); in MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable() 270 PetscCall(PetscLLCondensedAddSorted(1, &j, lnk, lnkbt)); in MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable() 273 apnz = lnk[0]; in MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable() 283 PetscCall(PetscLLCondensedClean(pN, apnz, current_space->array, lnk, lnkbt)); in MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable() 296 PetscCall(PetscLLDestroy(lnk, lnkbt)); in MatMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable() 809 …PetscInt i, pnz, row, *api, *apj, *Jptr, apnz, nspacedouble = 0, j, nzi, *lnk, apnz_ma… in MatMatMultSymbolic_MPIAIJ_MPIAIJ() local 847 PetscCall(PetscLLCondensedCreate_Scalable(lsize, &lnk)); in MatMatMultSymbolic_MPIAIJ_MPIAIJ() [all …]
|
| H A D | mpiptap.c | 224 PetscInt *lnk, i, k, pnz, row; in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable() local 286 PetscCall(PetscLLCondensedCreate_Scalable(Crmax, &lnk)); in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable() 310 PetscCall(PetscLLCondensedAddSorted_Scalable(pnz, Jptr, lnk)); in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable() 322 PetscCall(PetscLLCondensedAddSorted_Scalable(pnz, Jptr, lnk)); in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable() 325 apnz = lnk[0]; in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable() 335 PetscCall(PetscLLCondensedClean_Scalable(apnz, current_space->array, lnk)); in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable() 345 PetscCall(PetscLLCondensedDestroy_Scalable(lnk)); in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable() 511 PetscCall(PetscLLCondensedCreate_Scalable(Crmax, &lnk)); in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable() 516 PetscCall(PetscLLCondensedAddSorted_Scalable(nzi, Jptr, lnk)); in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable() 523 PetscCall(PetscLLCondensedAddSorted_Scalable(nzi, Jptr, lnk)); in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable() [all …]
|
| H A D | mpiaij.c | 4775 PetscInt k, anzi, *bi, *bj, *lnk, nlnk, arow, bnzi; in MatCreateMPIAIJSumSeqAIJSymbolic() local 4905 PetscCall(PetscLLCreate(N, N, nlnk, lnk, lnkbt)); in MatCreateMPIAIJSumSeqAIJSymbolic() 4931 PetscCall(PetscLLAddSorted(anzi, aj, N, &nlnk, lnk, lnkbt)); in MatCreateMPIAIJSumSeqAIJSymbolic() 4938 PetscCall(PetscLLAddSorted(anzi, aj, N, &nlnk, lnk, lnkbt)); in MatCreateMPIAIJSumSeqAIJSymbolic() 4949 PetscCall(PetscLLClean(N, N, bnzi, lnk, current_space->array, lnkbt)); in MatCreateMPIAIJSumSeqAIJSymbolic() 4963 PetscCall(PetscLLDestroy(lnk, lnkbt)); in MatCreateMPIAIJSumSeqAIJSymbolic()
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | matptap.c | 71 PetscInt i, j, k, ptnzi, arow, anzj, ptanzi, prow, pnzj, cnzi, nlnk, *lnk; in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy() local 91 PetscCall(PetscLLCreate(pn, pn, nlnk, lnk, lnkbt)); in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy() 121 PetscCall(PetscLLAddSorted(pnzj, pjj, pn, &nlnk, lnk, lnkbt)); in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy() 133 PetscCall(PetscLLClean(pn, pn, cnzi, lnk, current_space->array, lnkbt)); in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy() 151 PetscCall(PetscLLDestroy(lnk, lnkbt)); in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
|
| H A D | aijfact.c | 63 PetscInt nlnk, *lnk, k, **bi_ptr; in MatLUFactorSymbolic_SeqAIJ() local 84 PetscCall(PetscLLCreate(n, n, nlnk, lnk, lnkbt)); in MatLUFactorSymbolic_SeqAIJ() 99 PetscCall(PetscLLAddPerm(nnz, ajtmp, ic, n, &nlnk, lnk, lnkbt)); in MatLUFactorSymbolic_SeqAIJ() 103 row = lnk[n]; in MatLUFactorSymbolic_SeqAIJ() 107 PetscCall(PetscLLAddSortedLU(ajtmp, row, &nlnk, lnk, lnkbt, i, nzbd, im)); in MatLUFactorSymbolic_SeqAIJ() 109 row = lnk[row]; in MatLUFactorSymbolic_SeqAIJ() 117 k = lnk[n]; in MatLUFactorSymbolic_SeqAIJ() 120 k = lnk[k]; in MatLUFactorSymbolic_SeqAIJ() 133 PetscCall(PetscLLClean(n, n, nzi, lnk, current_space->array, lnkbt)); in MatLUFactorSymbolic_SeqAIJ() 147 PetscCall(PetscLLDestroy(lnk, lnkbt)); in MatLUFactorSymbolic_SeqAIJ() [all …]
|
| H A D | matmatmult.c | 150 PetscInt i, j, anzi, brow, bnzj, cnzi, *bj, *aj, *lnk, ndouble = 0, Crmax; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed() local 168 PetscCall(PetscLLCondensedCreate(Crmax, bn, &lnk, &lnkbt)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed() 184 PetscCall(PetscLLCondensedAddSorted(bnzj, bj, lnk, lnkbt)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed() 187 if (C->force_diagonals) PetscCall(PetscLLCondensedAddSorted(1, &i, lnk, lnkbt)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed() 188 cnzi = lnk[0]; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed() 198 PetscCall(PetscLLCondensedClean(bn, cnzi, current_space->array, lnk, lnkbt)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed() 212 PetscCall(PetscLLCondensedDestroy(lnk, lnkbt)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed() 387 PetscInt i, j, anzi, brow, bnzj, cnzi, *bj, *aj, *lnk, ndouble = 0, Crmax; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable_fast() local 403 PetscCall(PetscLLCondensedCreate_fast(Crmax, &lnk)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable_fast() 418 PetscCall(PetscLLCondensedAddSorted_fast(bnzj, bj, lnk)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable_fast() [all …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 229 PetscInt nlnk, *lnk, ncols, *cols, *uj, **ui_ptr, *uj_ptr, *udiag; in MatCholeskyFactorSymbolic_SeqSBAIJ() local 264 PetscCall(PetscLLCreate(mbs, mbs, nlnk, lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqSBAIJ() 279 PetscCall(PetscLLAdd(ncols, cols, mbs, &nlnk, lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqSBAIJ() 292 PetscCall(PetscLLAddSorted(ncols, uj_ptr, mbs, &nlnk, lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqSBAIJ() 314 PetscCall(PetscLLClean(mbs, mbs, nzk, lnk, current_space->array, lnkbt)); in MatCholeskyFactorSymbolic_SeqSBAIJ() 338 PetscCall(PetscLLDestroy(lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqSBAIJ() 394 PetscInt nlnk, *lnk, ncols, *cols, *uj, **ui_ptr, *uj_ptr; in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace() local 439 PetscCall(PetscLLCreate(mbs, mbs, nlnk, lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace() 453 PetscCall(PetscLLAdd(ncols, cols, mbs, &nlnk, lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace() 466 PetscCall(PetscLLAddSorted(ncols, uj_ptr, mbs, &nlnk, lnk, lnkbt)); in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace() [all …]
|
| H A D | sbaijfact2.c | 2478 PetscInt nlnk, *lnk, *lnk_lvl = NULL, ncols, *uj, **uj_ptr, **uj_lvl_ptr; in MatICCFactorSymbolic_SeqSBAIJ() local 2532 PetscCall(PetscIncompleteLLCreate(am, am, nlnk, lnk, lnk_lvl, lnkbt)); in MatICCFactorSymbolic_SeqSBAIJ() 2549 PetscCall(PetscIncompleteLLInit(ncols, cols, am, rip, &nlnk, lnk, lnk_lvl, lnkbt)); in MatICCFactorSymbolic_SeqSBAIJ() 2567 PetscCall(PetscICCLLAddSorted(ncols, cols, levels, uj, am, &nlnk, lnk, lnk_lvl, lnkbt, j)); in MatICCFactorSymbolic_SeqSBAIJ() 2592 …PetscCall(PetscIncompleteLLClean(am, am, nzk, lnk, lnk_lvl, current_space->array, current_space_lv… in MatICCFactorSymbolic_SeqSBAIJ() 2620 PetscCall(PetscIncompleteLLDestroy(lnk, lnkbt)); in MatICCFactorSymbolic_SeqSBAIJ() 2680 PetscInt nlnk, *lnk, *lnk_lvl = NULL, ncols, *cols_lvl, *uj, **uj_ptr, **uj_lvl_ptr; in MatICCFactorSymbolic_SeqSBAIJ_inplace() local 2730 PetscCall(PetscIncompleteLLCreate(am, am, nlnk, lnk, lnk_lvl, lnkbt)); in MatICCFactorSymbolic_SeqSBAIJ_inplace() 2746 PetscCall(PetscIncompleteLLInit(ncols, cols, am, rip, &nlnk, lnk, lnk_lvl, lnkbt)); in MatICCFactorSymbolic_SeqSBAIJ_inplace() 2763 … PetscCall(PetscICCLLAddSorted(ncols, cols, levels, cols_lvl, am, &nlnk, lnk, lnk_lvl, lnkbt, j)); in MatICCFactorSymbolic_SeqSBAIJ_inplace() [all …]
|
| /petsc/src/dm/impls/composite/ |
| H A D | pack.c | 82 struct DMCompositeLink *lnk = com->next; in DMView_Composite() local 88 for (i = 0; lnk; lnk = lnk->next, i++) { in DMView_Composite() 89 …erASCIIPrintf(v, "Link %" PetscInt_FMT ": DM of type %s\n", i, ((PetscObject)lnk->dm)->type_name)); in DMView_Composite() 91 PetscCall(DMView(lnk->dm, v)); in DMView_Composite()
|