Home
last modified time | relevance | path

Searched refs:rmax (Results 1 – 25 of 28) sorted by relevance

12

/petsc/src/mat/impls/aij/seq/crl/
H A Dcrl.c37 PetscInt i, j, rmax = a->rmax, *icols, *ilen = a->ilen; in MatSeqAIJCRL_create_aijcrl() local
44 aijcrl->rmax = rmax; in MatSeqAIJCRL_create_aijcrl()
47 PetscCall(PetscMalloc2(rmax * m, &aijcrl->acols, rmax * m, &aijcrl->icols)); in MatSeqAIJCRL_create_aijcrl()
55 for (; j < rmax; j++) { /* empty column entries */ in MatSeqAIJCRL_create_aijcrl()
60 …iply %g. Rmax= %" PetscInt_FMT "\n", 1.0 - ((double)a->nz) / PetscMax((double)rmax * m, 1), rmax)); in MatSeqAIJCRL_create_aijcrl()
90 PetscInt rmax = aijcrl->rmax, *icols = aijcrl->icols; in MatMult_AIJCRL() local
115 fortranmultcrl_(&m, &rmax, x, y, icols, acols); in MatMult_AIJCRL()
125 for (i = 1; i < rmax; i++) { in MatMult_AIJCRL()
H A Dcrl.h8 PetscInt rmax; /* maximum number of columns in a row */ member
/petsc/src/mat/impls/aij/seq/crl/ftn-kernels/
H A Dfmultcrl.F9010 PetscInt, intent(in) :: m, rmax, icols(m, rmax)
11 PetscScalar, intent(in) :: x(0:m - 1), acols(m, rmax)
17 do i = 2, rmax
/petsc/src/mat/impls/aij/mpi/crl/
H A Dmcrl.c43 PetscInt i, j, rmax = 0, *icols, *ailen = Aij->ilen, *bilen = Bij->ilen; in MatMPIAIJCRL_create_aijcrl() local
48 for (i = 0; i < m; i++) rmax = PetscMax(rmax, ailen[i] + bilen[i]); in MatMPIAIJCRL_create_aijcrl()
51 aijcrl->rmax = rmax; in MatMPIAIJCRL_create_aijcrl()
54 PetscCall(PetscMalloc2(rmax * m, &aijcrl->acols, rmax * m, &aijcrl->icols)); in MatMPIAIJCRL_create_aijcrl()
66 for (; j < rmax; j++) { /* empty column entries */ in MatMPIAIJCRL_create_aijcrl()
71 …oduced for vectorized multiply %g\n", 1.0 - ((double)aijcrl->nz) / PetscMax((double)rmax * m, 1))); in MatMPIAIJCRL_create_aijcrl()
/petsc/src/vec/vec/utils/
H A Dvecstash.c63 stash->rmax = 0; in VecStashCreate_Private()
115 stash->rmax = 0; in VecStashScatterEnd_Private()
309 stash->rmax = nmax; in VecStashScatterBegin_Private()
363 *rows = stash->rindices + i2 * stash->rmax; in VecStashScatterGetMesg_Private()
364 *vals = stash->rvalues + i1 * bs * stash->rmax; in VecStashScatterGetMesg_Private()
/petsc/src/mat/impls/sell/mpi/
H A Dmpisell.h17 PetscInt rmax; /* maximum message length */ member
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaij.c578 PetscInt *rp, k, low, high, t, ii, jj, row, nrow, i, col, l, rmax, N, lastcol = -1; in MatSetValuesBlocked_SeqSBAIJ() local
594 rmax = imax[row]; in MatSetValuesBlocked_SeqSBAIJ()
648 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, row, col, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValuesBlocked_SeqSBAIJ()
679 PetscInt mbs = a->mbs, bs2 = a->bs2, rmax = 0; in MatAssemblyEnd_SeqSBAIJ() local
685 if (m) rmax = ailen[0]; in MatAssemblyEnd_SeqSBAIJ()
689 rmax = PetscMax(rmax, ailen[i]); in MatAssemblyEnd_SeqSBAIJ()
711 PetscCall(PetscInfo(A, "Most nonzeros blocks in any row is %" PetscInt_FMT "\n", rmax)); in MatAssemblyEnd_SeqSBAIJ()
717 a->rmax = rmax; in MatAssemblyEnd_SeqSBAIJ()
740 PetscInt *rp, k, low, high, t, ii, row, nrow, i, col, l, rmax, N, lastcol = -1; in MatSetValues_SeqSBAIJ() local
754 rmax = imax[brow]; /* maximum space allocated for this row */ in MatSetValues_SeqSBAIJ()
[all …]
/petsc/lib/petsc/bin/maint/
H A Dtoclapack.sh626 static real rmin, rmax, rmach;
1653 /* = 'O' or 'o', DLAMCH := rmax */
1666 /* rmax = overflow threshold - (base**emax)*(1-eps) */
1685 dlamc2_(&beta, &it, &lrnd, &eps, &imin, &rmin, &imax, &rmax);
1701 small = 1. / rmax;
1730 rmach = rmax;
1953 doublereal *rmax)
2253 *rmax = lrmax;
2405 logical *ieee, integer *emax, doublereal *rmax)
2569 *rmax = y;
[all …]
/petsc/src/mat/impls/blockmat/seq/
H A Dblockmat.c262 …MatSeqXAIJReallocateAIJ(A, a->mbs, 1, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValues_BlockMat()
600 PetscInt m = a->mbs, *ip, N, *ailen = a->ilen, rmax = 0; in MatAssemblyEnd_BlockMat() local
606 if (m) rmax = ailen[0]; /* determine row with most nonzeros */ in MatAssemblyEnd_BlockMat()
610 rmax = PetscMax(rmax, ailen[i]); in MatAssemblyEnd_BlockMat()
636 PetscCall(PetscInfo(A, "Maximum nonzeros in any row is %" PetscInt_FMT "\n", rmax)); in MatAssemblyEnd_BlockMat()
641 a->rmax = rmax; in MatAssemblyEnd_BlockMat()
/petsc/src/snes/interface/
H A Dsnesut.c480 PetscReal rmax, pwork; in SNESMonitorRange_Private() local
486 PetscCall(VecNorm(resid, NORM_INFINITY, &rmax)); in SNESMonitorRange_Private()
491 for (i = 0; i < n; i++) pwork += (PetscAbsScalar(r[i]) > .20 * rmax); in SNESMonitorRange_Private()
/petsc/src/mat/impls/baij/seq/
H A Dbaij.c2108 PetscInt *rp, k, low, high, t, ii, jj, row, nrow, i, col, l, rmax, N, lastcol = -1; in MatSetValuesBlocked_SeqBAIJ() local
2127 rmax = imax[row]; in MatSetValuesBlocked_SeqBAIJ()
2184 …MatSeqXAIJReallocateAIJ_structure_only(A, a->mbs, bs2, nrow, row, col, rmax, ai, aj, rp, imax, non… in MatSetValuesBlocked_SeqBAIJ()
2186 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, row, col, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValuesBlocked_SeqBAIJ()
2219 PetscInt mbs = a->mbs, bs2 = a->bs2, rmax = 0; in MatAssemblyEnd_SeqBAIJ() local
2226 if (m) rmax = ailen[0]; in MatAssemblyEnd_SeqBAIJ()
2230 rmax = PetscMax(rmax, ailen[i]); in MatAssemblyEnd_SeqBAIJ()
2262 PetscCall(PetscInfo(A, "Most nonzeros blocks in any row is %" PetscInt_FMT "\n", rmax)); in MatAssemblyEnd_SeqBAIJ()
2267 a->rmax = rmax; in MatAssemblyEnd_SeqBAIJ()
2446 PetscInt *rp, k, low, high, t, ii, row, nrow, i, col, l, rmax, N, lastcol = -1; in MatSetValues_SeqBAIJ() local
[all …]
/petsc/src/ksp/pc/impls/gamg/
H A Dclassical.c90 PetscReal rmax; in PCGAMGCreateGraph_Classical() local
110 rmax = 0.; in PCGAMGCreateGraph_Classical()
113 if (PetscRealPart(-rval[c]) > rmax && rcol[c] != r) rmax = PetscRealPart(-rval[c]); in PCGAMGCreateGraph_Classical()
115 Amax[r - s] = rmax; in PCGAMGCreateGraph_Classical()
/petsc/src/mat/impls/sbaij/mpi/
H A Dmpisbaij.c206 rmax = aimax[brow]; \
229 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, aimax, a->none…
246 rmax = bimax[brow]; \
269 …MatSeqXAIJReallocateAIJ(B, b->mbs, bs2, nrow, brow, bcol, rmax, ba, bi, bj, rp, ap, bimax, b->none…
305 PetscInt *rp, ii, nrow, _i, rmax, N, brow, bcol; in MatSetValues_MPISBAIJ() local
403 PetscInt *rp, low, high, t, ii, jj, nrow, i, rmax, N; in MatSetValuesBlocked_SeqSBAIJ_Inlined() local
417 rmax = imax[row]; in MatSetValuesBlocked_SeqSBAIJ_Inlined()
458 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, row, col, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValuesBlocked_SeqSBAIJ_Inlined()
486 PetscInt *rp, low, high, t, ii, jj, nrow, i, rmax, N; in MatSetValuesBlocked_SeqBAIJ_Inlined() local
496 rmax = imax[row]; in MatSetValuesBlocked_SeqBAIJ_Inlined()
[all …]
/petsc/src/mat/impls/aij/seq/
H A Daij.c405 PetscInt *rp, k, low, high, t, ii, row, nrow, i, col, l, rmax, N; in MatSetValues_SeqAIJ() local
420 rmax = imax[row]; in MatSetValues_SeqAIJ()
456 …MatSeqXAIJReallocateAIJ_structure_only(A, A->rmap->n, 1, nrow, row, col, rmax, ai, aj, rp, imax, n… in MatSetValues_SeqAIJ()
458 …MatSeqXAIJReallocateAIJ(A, A->rmap->n, 1, nrow, row, col, rmax, aa, ai, aj, rp, ap, imax, nonew, M… in MatSetValues_SeqAIJ()
1101 PetscInt m = A->rmap->n, *ip, N, *ailen = a->ilen, rmax = 0; in MatAssemblyEnd_SeqAIJ() local
1113 if (m) rmax = ailen[0]; /* determine row with most nonzeros */ in MatAssemblyEnd_SeqAIJ()
1117 rmax = PetscMax(rmax, ailen[i]); in MatAssemblyEnd_SeqAIJ()
1146 PetscCall(PetscInfo(A, "Maximum nonzeros in any row is %" PetscInt_FMT "\n", rmax)); in MatAssemblyEnd_SeqAIJ()
1151 a->rmax = rmax; in MatAssemblyEnd_SeqAIJ()
4538 *nz = aij->rmax; in MatSeqAIJGetMaxRowNonzeros()
[all …]
H A Dmatmatmult.c48 for (ii = 0, aij->nonzerorowcnt = 0, aij->rmax = 0; ii < m; ii++) { in MatSetSeqAIJWithArrays_private()
51 aij->rmax = PetscMax(aij->rmax, rnz); in MatSetSeqAIJWithArrays_private()
607 PetscCall(PetscHeapCreate(a->rmax, &h)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Heap()
608 PetscCall(PetscMalloc1(a->rmax, &bb)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Heap()
711 PetscCall(PetscHeapCreate(a->rmax, &h)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_BTHeap()
712 PetscCall(PetscMalloc1(a->rmax, &bb)); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_BTHeap()
H A Daij.h58 PetscInt rmax; /* max nonzeros in any row */ \
/petsc/src/mat/utils/
H A Dmatstash.c338 …, PetscInt n, const PetscInt idxn[], const PetscScalar values[], PetscInt rmax, PetscInt cmax, Pet… in MatStashValuesRowBlocked_Private() argument
389 …, PetscInt n, const PetscInt idxn[], const PetscScalar values[], PetscInt rmax, PetscInt cmax, Pet… in MatStashValuesColBlocked_Private() argument
412 vals += rmax * bs; in MatStashValuesColBlocked_Private()
/petsc/src/mat/impls/baij/mpi/
H A Dmpibaij.c165 rmax = aimax[brow]; \
188 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, aimax, a->none…
205 rmax = bimax[brow]; \
228 …MatSeqXAIJReallocateAIJ(B, b->mbs, bs2, nrow, brow, bcol, rmax, ba, bi, bj, rp, ap, bimax, b->none…
261 PetscInt *rp, ii, nrow, _i, rmax, N, brow, bcol; in MatSetValues_MPIBAIJ() local
329 PetscInt *rp, low, high, t, ii, jj, nrow, i, rmax, N; in MatSetValuesBlocked_SeqBAIJ_Inlined() local
339 rmax = imax[row]; in MatSetValuesBlocked_SeqBAIJ_Inlined()
381 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, row, col, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValuesBlocked_SeqBAIJ_Inlined()
3588 PetscInt *bindx, rmax = a->rmax, j; in MatCreateMPIMatConcatenateSeqMat_MPIBAIJ() local
3597 PetscCall(PetscMalloc1(rmax, &bindx)); in MatCreateMPIMatConcatenateSeqMat_MPIBAIJ()
/petsc/src/ksp/ksp/interface/
H A Diterativ.c298 PetscReal rmax, pwork; in KSPMonitorRange_Private() local
303 PetscCall(VecNorm(resid, NORM_INFINITY, &rmax)); in KSPMonitorRange_Private()
308 for (i = 0; i < n; ++i) pwork += (PetscAbsScalar(r[i]) > .20 * rmax); in KSPMonitorRange_Private()
/petsc/src/mat/impls/aij/mpi/
H A Dmpiaij.h51 PetscInt rmax; /* maximum message length */ \
H A Dmpiov.c1226 PetscInt *lens, rmax, ncols, *cols, Crow; in MatCreateSubMatrices_MPIAIJ_SingleIS_Local() local
1646 rmax = 0; in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
1647 for (PetscMPIInt i = 0; i < nrow; i++) rmax = PetscMax(rmax, lens[i]); in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
1654 rmax = subc->rmax; in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
1684 PetscCall(PetscMalloc3(nrqs, &rbuf4, rmax, &subcols, rmax, &subvals)); in MatCreateSubMatrices_MPIAIJ_SingleIS_Local()
H A Dmpimatmatmult.c1778 PetscInt rstart = P->cmap->rstart, rmax, Armax; in MatTransposeMatMultSymbolic_MPIAIJ_MPIAIJ() local
1988 rmax = 0; in MatTransposeMatMultSymbolic_MPIAIJ_MPIAIJ()
2034 if (nnz > rmax) rmax = nnz; in MatTransposeMatMultSymbolic_MPIAIJ_MPIAIJ()
/petsc/include/petsc/private/
H A Dvecimpl.h145 PetscInt rmax; /* maximum message length */ member
/petsc/src/mat/impls/aij/seq/seqhipsparse/
H A Daijhipsparse.hip.cxx2513 PetscCall(PetscInfo(C, "Maximum nonzeros in any row is %" PetscInt_FMT "\n", c->rmax)); in MatProductNumeric_SeqAIJHIPSPARSE_SeqAIJHIPSPARSE()
2840 c->rmax = 0; in MatProductSymbolic_SeqAIJHIPSPARSE_SeqAIJHIPSPARSE()
2845 c->rmax = PetscMax(c->rmax, nn); in MatProductSymbolic_SeqAIJHIPSPARSE_SeqAIJHIPSPARSE()
4288 c->rmax = 0; in MatSeqAIJHIPSPARSEMergeMats()
4293 c->rmax = PetscMax(c->rmax, nn); in MatSeqAIJHIPSPARSEMergeMats()
/petsc/src/mat/impls/aij/seq/seqcusparse/
H A Daijcusparse.cu3015 PetscCall(PetscInfo(C, "Maximum nonzeros in any row is %" PetscInt_FMT "\n", c->rmax)); in MatProductNumeric_SeqAIJCUSPARSE_SeqAIJCUSPARSE()
3367 c->rmax = 0; in MatProductSymbolic_SeqAIJCUSPARSE_SeqAIJCUSPARSE()
3372 c->rmax = PetscMax(c->rmax, nn); in MatProductSymbolic_SeqAIJCUSPARSE_SeqAIJCUSPARSE()
4965 c->rmax = 0; in MatSeqAIJCUSPARSEMergeMats()
4970 c->rmax = PetscMax(c->rmax, nn); in MatSeqAIJCUSPARSEMergeMats()

12