| /petsc/src/mat/utils/ |
| H A D | zerodiag.c | 72 PetscInt prow, k, nz, n, repl, *j, *col, *row, m, *icol, nnz, *jj, kk; in MatReorderForNonzeroDiagonal_SeqAIJ() local
85 for (prow = 0; prow < n; prow++) { in MatReorderForNonzeroDiagonal_SeqAIJ()
86 PetscCall(MatGetRow_SeqAIJ(mat, row[prow], &nz, &j, &v)); in MatReorderForNonzeroDiagonal_SeqAIJ()
88 if (icol[j[k]] == prow) break; in MatReorderForNonzeroDiagonal_SeqAIJ()
97 if (icol[j[k]] > prow && PetscAbsScalar(v[k]) > repla) { in MatReorderForNonzeroDiagonal_SeqAIJ()
100 SWAP(icol[col[prow]], icol[col[repl]]); in MatReorderForNonzeroDiagonal_SeqAIJ()
101 SWAP(col[prow], col[repl]); in MatReorderForNonzeroDiagonal_SeqAIJ()
111 if (icol[j[k]] < prow && PetscAbsScalar(v[k]) > repla) { in MatReorderForNonzeroDiagonal_SeqAIJ()
116 if (icol[jj[kk]] == prow && PetscAbsScalar(vv[kk]) > abstol) { in MatReorderForNonzeroDiagonal_SeqAIJ()
118 SWAP(icol[col[prow]], icol[col[repl]]); in MatReorderForNonzeroDiagonal_SeqAIJ()
[all …]
|
| /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
117 prow = *ptaj++; in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
118 pnzj = pi[prow + 1] - pi[prow]; in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
119 pjj = pj + pi[prow]; in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
196 PetscInt i, j, k, anzi, pnzi, apnzj, nextap, pnzj, prow, crow; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy() local
227 prow = *aj++; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy()
228 pnzj = pi[prow + 1] - pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy()
229 pjj = pj + pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy()
230 paj = pa + pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy()
|
| H A D | aijfact.c | 1545 PetscInt *bdiag, prow, fm, nzbd, reallocs = 0, dcount = 0; in MatILUFactorSymbolic_SeqAIJ() local
1615 prow = lnk[n]; in MatILUFactorSymbolic_SeqAIJ()
1616 while (prow < i) { in MatILUFactorSymbolic_SeqAIJ()
1617 nnz = bdiag[prow]; in MatILUFactorSymbolic_SeqAIJ()
1618 cols = bj_ptr[prow] + nnz + 1; in MatILUFactorSymbolic_SeqAIJ()
1619 cols_lvl = bjlvl_ptr[prow] + nnz + 1; in MatILUFactorSymbolic_SeqAIJ()
1620 nnz = bi[prow + 1] - bi[prow] - nnz - 1; in MatILUFactorSymbolic_SeqAIJ()
1621 …etscCall(PetscILULLAddSorted(nnz, cols, levels, cols_lvl, prow, &nlnk, lnk, lnk_lvl, lnkbt, prow)); in MatILUFactorSymbolic_SeqAIJ()
1623 prow = lnk[prow]; in MatILUFactorSymbolic_SeqAIJ()
2058 PetscInt jmin, jmax, nzk, k, j, *jl, prow, *il, nextprow; in MatICCFactorSymbolic_SeqAIJ() local
[all …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 44 PetscInt m, reallocs = 0, prow; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR() local
107 prow = k; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
108 prow = jl[prow]; /* next pivot row (== mbs for symbolic factorization) */ in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
110 while (prow < k) { in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
112 jmin = iu[prow] + 1; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
113 jmax = iu[prow + 1]; in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
128 prow = jl[prow]; /* next pivot row */ in MatCholeskyFactorSymbolic_SeqSBAIJ_MSR()
227 PetscInt i, mbs = a->mbs, bs = A->rmap->bs, reallocs = 0, prow; in MatCholeskyFactorSymbolic_SeqSBAIJ() local
283 prow = jl[k]; /* 1st pivot row */ in MatCholeskyFactorSymbolic_SeqSBAIJ()
285 while (prow < k) { in MatCholeskyFactorSymbolic_SeqSBAIJ()
[all …]
|
| H A D | sbaijfact2.c | 2286 PetscInt incrlev, *lev, shift, prow, nz; in MatICCFactorSymbolic_SeqSBAIJ_MSR() local
2343 prow = k; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2344 prow = prowl[prow]; /* next pivot row (== 0 for symbolic factorization) */ in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2346 while (prow < k) { in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2348 jmin = iu[prow] + 1; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2349 jmax = iu[prow + 1]; in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2367 prow = prowl[prow]; /* next pivot row */ in MatICCFactorSymbolic_SeqSBAIJ_MSR()
2477 PetscInt jmin, jmax, nzk, k, j, *jl, prow, *il, nextprow; in MatICCFactorSymbolic_SeqSBAIJ() local
2553 prow = jl[k]; /* 1st pivot row */ in MatICCFactorSymbolic_SeqSBAIJ()
2555 while (prow < k) { in MatICCFactorSymbolic_SeqSBAIJ()
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact2.c | 318 PetscInt *bdiag, prow, fm, nzbd, reallocs = 0, dcount = 0; in MatILUFactorSymbolic_SeqBAIJ() local
416 prow = lnk[n]; in MatILUFactorSymbolic_SeqBAIJ()
417 while (prow < i) { in MatILUFactorSymbolic_SeqBAIJ()
418 nnz = bdiag[prow]; in MatILUFactorSymbolic_SeqBAIJ()
419 cols = bj_ptr[prow] + nnz + 1; in MatILUFactorSymbolic_SeqBAIJ()
420 cols_lvl = bjlvl_ptr[prow] + nnz + 1; in MatILUFactorSymbolic_SeqBAIJ()
421 nnz = bi[prow + 1] - bi[prow] - nnz - 1; in MatILUFactorSymbolic_SeqBAIJ()
423 …etscCall(PetscILULLAddSorted(nnz, cols, levels, cols_lvl, prow, &nlnk, lnk, lnk_lvl, lnkbt, prow)); in MatILUFactorSymbolic_SeqBAIJ()
425 prow = lnk[prow]; in MatILUFactorSymbolic_SeqBAIJ()
|
| H A D | baijfact.c | 1024 PetscInt jmin, jmax, nzk, k, j, *jl, prow, *il, nextprow; in MatICCFactorSymbolic_SeqBAIJ() local
1121 prow = jl[k]; /* 1st pivot row */ in MatICCFactorSymbolic_SeqBAIJ()
1123 while (prow < k) { in MatICCFactorSymbolic_SeqBAIJ()
1124 nextprow = jl[prow]; in MatICCFactorSymbolic_SeqBAIJ()
1127 jmin = il[prow] + 1; /* index of the 2nd nzero entry in U(prow,k:am-1) */ in MatICCFactorSymbolic_SeqBAIJ()
1128 jmax = ui[prow + 1]; in MatICCFactorSymbolic_SeqBAIJ()
1130 i = jmin - ui[prow]; in MatICCFactorSymbolic_SeqBAIJ()
1131 cols = uj_ptr[prow] + i; /* points to the 2nd nzero entry in U(prow,k:am-1) */ in MatICCFactorSymbolic_SeqBAIJ()
1132 for (j = 0; j < ncols; j++) cols_lvl[j] = *(uj_lvl_ptr[prow] + i + j); in MatICCFactorSymbolic_SeqBAIJ()
1138 il[prow] = jmin; in MatICCFactorSymbolic_SeqBAIJ()
[all …]
|
| /petsc/src/mat/impls/maij/ |
| H A D | maij.c | 594 …PetscInt i, j, k, pshift, poffset, anzi, pnzi, apnzj, nextap, pnzj, prow, crow, *apj, *apj… in MatPtAPNumeric_SeqAIJ_SeqMAIJ() local
613 prow = *aj++ / ppdof; /* integer division */ in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
614 pnzj = pi[prow + 1] - pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
615 pjj = pj + pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
616 paj = pa + pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
634 prow = i / ppdof; /* integer division */ in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
636 poffset = pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
637 pnzi = pi[prow + 1] - poffset; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
677 PetscInt i, j, k, dof, pshift, ptnzi, arow, anzj, ptanzi, prow, pnzj, cnzi, cn; in MatPtAPSymbolic_SeqAIJ_SeqMAIJ() local
729 prow = (*ptaj++) / ppdof; /* integer division */ in MatPtAPSymbolic_SeqAIJ_SeqMAIJ()
[all …]
|
| /petsc/include/petsc/private/ |
| H A D | petscscalapack.h | 53 BLAS_EXTERN PetscBLASInt Cblacs_pnum(PetscBLASInt context, PetscBLASInt prow, PetscBLASInt pcol);
|
| H A D | matimpl.h | 1205 …ICT nlnk, PetscInt *PETSC_RESTRICT lnk, PetscInt *PETSC_RESTRICT lnklvl, PetscBT bt, PetscInt prow) in PetscILULLAddSorted() argument
1208 …Add_Private(nidx, idx, level, idxlvl, idx_start, nlnk, lnk, lnklvl, bt, lnklvl[prow], PETSC_TRUE)); in PetscILULLAddSorted()
|