Home
last modified time | relevance | path

Searched refs:bm (Results 1 – 22 of 22) sorted by relevance

/petsc/src/sys/classes/bm/interfaces/
H A Dbm.c79 PetscErrorCode PetscBenchReset(PetscBench bm) in PetscBenchReset() argument
82 PetscValidHeaderSpecific(bm, BM_CLASSID, 1); in PetscBenchReset()
83 …PetscCall(PetscLogHandlerDestroy(&bm->lhdlr)); // Temporarily here until PetscLogHandlerReset() ex… in PetscBenchReset()
84 PetscTryTypeMethod(bm, reset); in PetscBenchReset()
85 bm->setupcalled = PETSC_FALSE; in PetscBenchReset()
101 PetscErrorCode PetscBenchDestroy(PetscBench *bm) in PetscBenchDestroy() argument
104 PetscAssertPointer(bm, 1); in PetscBenchDestroy()
105 if (!*bm) PetscFunctionReturn(PETSC_SUCCESS); in PetscBenchDestroy()
106 PetscValidHeaderSpecific(*bm, BM_CLASSID, 1); in PetscBenchDestroy()
107 if (--((PetscObject)*bm)->refct > 0) { in PetscBenchDestroy()
[all …]
/petsc/src/mat/utils/hpl/
H A Dhplcreate.c45 static PetscErrorCode PetscBenchSetUp_HPL(PetscBench bm) in PetscBenchSetUp_HPL() argument
47 PetscBench_HPL *hp = (PetscBench_HPL *)bm->data; in PetscBenchSetUp_HPL()
51 PetscCallMPI(MPI_Comm_rank(PetscObjectComm((PetscObject)bm), &rank)); in PetscBenchSetUp_HPL()
53 if (bm->size == PETSC_DECIDE) bm->size = 2500; in PetscBenchSetUp_HPL()
54 PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, bm->size, bm->size, NULL, &hp->A)); in PetscBenchSetUp_HPL()
61 static PetscErrorCode PetscBenchRun_HPL(PetscBench bm) in PetscBenchRun_HPL() argument
63 PetscBench_HPL *hp = (PetscBench_HPL *)bm->data; in PetscBenchRun_HPL()
67 PetscCallMPI(MPI_Comm_rank(PetscObjectComm((PetscObject)bm), &rank)); in PetscBenchRun_HPL()
73 static PetscErrorCode PetscBenchView_HPL(PetscBench bm, PetscViewer viewer) in PetscBenchView_HPL() argument
79 …PetscCall(PetscLogHandlerGetEventPerfInfo(bm->lhdlr, 0, MAT_LUFactor, &info)); // because symbolic… in PetscBenchView_HPL()
[all …]
/petsc/src/mat/tests/
H A Dtmpbench.c8 PetscBench bm; in main() local
13 PetscCall(PetscBenchCreate(PETSC_COMM_SELF, &bm)); in main()
14 PetscCall(PetscBenchSetType(bm, PETSCBMHPL)); in main()
15 PetscCall(PetscBenchSetFromOptions(bm)); in main()
16 PetscCall(PetscBenchSetUp(bm)); in main()
17 PetscCall(PetscBenchRun(bm)); in main()
18 PetscCall(PetscBenchView(bm, PETSC_VIEWER_STDOUT_SELF)); in main()
19 PetscCall(PetscBenchSetSize(bm, 5000)); in main()
20 PetscCall(PetscBenchRun(bm)); in main()
21 PetscCall(PetscBenchView(bm, PETSC_VIEWER_STDOUT_SELF)); in main()
[all …]
/petsc/doc/tutorials/
H A Dguide_to_examples_by_physics.md37 -\nabla\cdot \bm \sigma = \bm f
40 where $\bm\sigma$ is the stress tensor. Linear, isotropic elasticity governing infinitesimal strain…
43 …\nabla\cdot \left( \lambda I \operatorname{trace}(\bm\varepsilon) + 2\mu \bm\varepsilon \right) = …
46 where the strain tensor $\bm \varepsilon$ is given by
49 \bm \varepsilon = \frac{1}{2} \left(\nabla \bm u + (\nabla \bm u)^T \right)
52 where $\bm u$ is the infinitesimal displacement of the body. The resulting discretizations use PETS…
/petsc/src/mat/utils/
H A Dmultequal.c14 PetscInt am, an, bm, bn, k; in MatMultEqual_Private() local
30 PetscCall(MatGetLocalSize(B, &bm, &bn)); in MatMultEqual_Private()
31 …bm && an == bn, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Mat A,Mat B: local dim %" PetscInt_FMT " %" P… in MatMultEqual_Private()
121 PetscInt am, an, bm, bn, cm, cn, k; in MatMatMultEqual_Private() local
139 PetscCall(MatGetLocalSize(B, &bm, &bn)); in MatMatMultEqual_Private()
148 bn = bm; in MatMatMultEqual_Private()
149 bm = tt; in MatMatMultEqual_Private()
151 …bm && am == cm && bn == cn, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Mat A, B, C local dim %" PetscInt… in MatMatMultEqual_Private()
442 PetscInt i, am, an, bm, bn, cm, cn; in MatProjMultEqual_Private() local
447 PetscCall(MatGetLocalSize(B, &bm, &bn)); in MatProjMultEqual_Private()
[all …]
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaij2.c1567 static PetscErrorCode MatMatMult_SeqSBAIJ_1_Private(Mat A, PetscScalar *b, PetscInt bm, PetscScalar… in MatMatMult_SeqSBAIJ_1_Private() argument
1590 x1 = xb[0 + k * bm]; in MatMatMult_SeqSBAIJ_1_Private()
1592 if (*idx != i) c[(*idx) + k * cm] += (aconj ? PetscConj(v[0]) : v[0]) * b[i + k * bm]; in MatMatMult_SeqSBAIJ_1_Private()
1602 static PetscErrorCode MatMatMult_SeqSBAIJ_2_Private(Mat A, PetscScalar *b, PetscInt bm, PetscScalar… in MatMatMult_SeqSBAIJ_2_Private() argument
1624 x1 = xb[0 + k * bm]; in MatMatMult_SeqSBAIJ_2_Private()
1625 x2 = xb[1 + k * bm]; in MatMatMult_SeqSBAIJ_2_Private()
1629 c[2 * (*idx) + 0 + k * cm] += v[0] * b[2 * i + k * bm] + v[1] * b[2 * i + 1 + k * bm]; in MatMatMult_SeqSBAIJ_2_Private()
1630 c[2 * (*idx) + 1 + k * cm] += v[2] * b[2 * i + k * bm] + v[3] * b[2 * i + 1 + k * bm]; in MatMatMult_SeqSBAIJ_2_Private()
1641 static PetscErrorCode MatMatMult_SeqSBAIJ_3_Private(Mat A, PetscScalar *b, PetscInt bm, PetscScalar… in MatMatMult_SeqSBAIJ_3_Private() argument
1663 x1 = xb[0 + k * bm]; in MatMatMult_SeqSBAIJ_3_Private()
[all …]
/petsc/src/ksp/pc/impls/tfs/
H A Dbitmask.c106 PetscErrorCode PCTFS_set_bit_mask(PetscInt *bm, PetscInt len, PetscInt val) in PCTFS_set_bit_mask() argument
115 cptr = (char *)bm; in PCTFS_set_bit_mask()
119 *bm = 0; in PCTFS_set_bit_mask()
120 bm++; in PCTFS_set_bit_mask()
/petsc/src/mat/impls/aij/seq/
H A Dmatrart.c109 …PetscInt cn = B->cmap->n, bm = B->rmap->n, col, i, j, n, *ai = a->i, *aj, am = A->rmap->… in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense() local
110 PetscInt am2 = 2 * am, am3 = 3 * am, bm4 = 4 * bm; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
117 …eck(bm == A->cmap->n, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Number columns in A %" PetscInt_FMT " n… in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
142 b2 = b1 + bm; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
143 b3 = b2 + bm; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
144 b4 = b3 + bm; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
199 b1 += bm; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
H A Dmattransposematmult.c110 PetscInt bm = B->rmap->n; in MatTMatTMultNumeric_SeqAIJ_SeqDense() local
114 for (i = 0; i < bm; i++) ctarray[j * bm + i] = Barray[j * blda + i]; in MatTMatTMultNumeric_SeqAIJ_SeqDense()
H A Dmatmatmult.c148 PetscInt am = A->rmap->N, bn = B->cmap->N, bm = B->rmap->N; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed() local
164 MatRowMergeMax_SeqAIJ(b, bm, ta); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
171 …PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, PetscIntSumTruncate(ai[am], bi[bm])), &… in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
229 afill = (PetscReal)ci[am] / (ai[am] + bi[bm]) + 1.e-5; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_LLCondensed()
384 PetscInt am = A->rmap->N, bn = B->cmap->N, bm = B->rmap->N; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable_fast() local
399 MatRowMergeMax_SeqAIJ(b, bm, ta); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable_fast()
406 …PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, PetscIntSumTruncate(ai[am], bi[bm])), &… in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable_fast()
466 afill = (PetscReal)ci[am] / (ai[am] + bi[bm]) + 1.e-5; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable_fast()
487 PetscInt am = A->rmap->N, bn = B->cmap->N, bm = B->rmap->N; in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable() local
502 MatRowMergeMax_SeqAIJ(b, bm, ta); in MatMatMultSymbolic_SeqAIJ_SeqAIJ_Scalable()
[all …]
H A Daij.c4116 …PetscInt *i, *j, m, n, p, q, nnz = 0, am = A->rmap->n, bm = B->rmap->n, an = A->cmap->n, … in MatSeqAIJKron_SeqAIJ() local
4128 PetscCall(PetscMalloc2(am * bm + 1, &i, a->i[am] * b->i[bm], &j)); in MatSeqAIJKron_SeqAIJ()
4130 PetscCall(MatSetSizes(newmat, am * bm, an * bn, am * bm, an * bn)); in MatSeqAIJKron_SeqAIJ()
4134 for (p = 0; p < bm; ++p) { in MatSeqAIJKron_SeqAIJ()
4135 i[m * bm + p + 1] = i[m * bm + p] + (a->i[m + 1] - a->i[m]) * (b->i[p + 1] - b->i[p]); in MatSeqAIJKron_SeqAIJ()
4150 for (p = 0; p < bm; ++p) { in MatSeqAIJKron_SeqAIJ()
/petsc/src/ksp/ksp/utils/lmvm/blas_cyclic/
H A Dblas_cyclic.c250 PetscBLASInt b_one = 1, blda, bm; in TRSVCyclic_Private() local
256 PetscCall(PetscBLASIntCast(m, &bm)); in TRSVCyclic_Private()
266 PetscBLASInt bn = bm - i_oldest; in TRSVCyclic_Private()
345 PetscBLASInt b_one = 1, blda, bm; in HEMVCyclic_Private() local
351 PetscCall(PetscBLASIntCast(m, &bm)); in HEMVCyclic_Private()
359 PetscBLASInt bn = bm - i_oldest; in HEMVCyclic_Private()
422 PetscBLASInt b_one = 1, blda, bm; in GEMVCyclic_Private() local
429 PetscCall(PetscBLASIntCast(m, &bm)); in GEMVCyclic_Private()
433 …PetscCallBLAS("BLASgemv", BLASgemv_(trans, &bm, &bm, &alpha, A, &blda, x, &b_one, &beta, y, &b_one… in GEMVCyclic_Private()
439 PetscBLASInt bn = bm - i_oldest; in GEMVCyclic_Private()
/petsc/src/ts/impls/rosw/
H A Drosw.c966 PetscReal A[4][4], Gamma[4][4], b[4], bm[4]; in TSRosWRegisterRos4() local
1061 …bm[3] = b[3] - e4 * gamma; /* using definition of E4 … in TSRosWRegisterRos4()
1062 bm[2] = (p32 - beta4jbetajp * bm[3]) / (beta32 * beta2p); /* fourth row of 7.18 */ in TSRosWRegisterRos4()
1063 bm[1] = (one / two - gamma - beta3p * bm[2] - beta4p * bm[3]) / beta2p; /* second row */ in TSRosWRegisterRos4()
1064 bm[0] = one - bm[1] - bm[2] - bm[3]; /* first row */ in TSRosWRegisterRos4()
1067 const PetscReal misfit = a2 * a2 * bm[1] + a3 * a3 * bm[2] + a4 * a4 * bm[3] - one / three; in TSRosWRegisterRos4()
1070 PetscCall(TSRosWRegister(name, 4, 4, &A[0][0], &Gamma[0][0], b, bm, 0, NULL)); in TSRosWRegisterRos4()
/petsc/src/mat/impls/baij/seq/
H A Dbaij2.c3362 static PetscErrorCode MatMatMult_SeqBAIJ_1_Private(Mat A, PetscScalar *b, PetscInt bm, PetscScalar … in MatMatMult_SeqBAIJ_1_Private() argument
3399 x1 = xb[0 + k * bm]; in MatMatMult_SeqBAIJ_1_Private()
3410 static PetscErrorCode MatMatMult_SeqBAIJ_2_Private(Mat A, PetscScalar *b, PetscInt bm, PetscScalar … in MatMatMult_SeqBAIJ_2_Private() argument
3448 x1 = xb[0 + k * bm]; in MatMatMult_SeqBAIJ_2_Private()
3449 x2 = xb[1 + k * bm]; in MatMatMult_SeqBAIJ_2_Private()
3462 static PetscErrorCode MatMatMult_SeqBAIJ_3_Private(Mat A, PetscScalar *b, PetscInt bm, PetscScalar … in MatMatMult_SeqBAIJ_3_Private() argument
3501 x1 = xb[0 + k * bm]; in MatMatMult_SeqBAIJ_3_Private()
3502 x2 = xb[1 + k * bm]; in MatMatMult_SeqBAIJ_3_Private()
3503 x3 = xb[2 + k * bm]; in MatMatMult_SeqBAIJ_3_Private()
3518 static PetscErrorCode MatMatMult_SeqBAIJ_4_Private(Mat A, PetscScalar *b, PetscInt bm, PetscScalar … in MatMatMult_SeqBAIJ_4_Private() argument
[all …]
/petsc/src/dm/dt/interface/
H A Ddt.c376 PetscBLASInt bm, bn, info; in PetscDTJacobianInverse_Internal() local
380 PetscCall(PetscBLASIntCast(m, &bm)); in PetscDTJacobianInverse_Internal()
396 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&bm, &bm, Jinvs, &bm, pivots, &info)); in PetscDTJacobianInverse_Internal()
398 PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&bm, Jinvs, &bm, pivots, W, &bm, &info)); in PetscDTJacobianInverse_Internal()
417 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&bm, &bm, JJT, &bm, pivots, &info)); in PetscDTJacobianInverse_Internal()
419 PetscCallBLAS("LAPACKgetri", LAPACKgetri_(&bm, JJT, &bm, pivots, W, &bm, &info)); in PetscDTJacobianInverse_Internal()
1427 PetscBLASInt bn, bm, ldz; /* bm will equal bn on exit */ in PetscDTSymmetricTridiagonalEigensolve() local
1446 …ACKstegr_(&jobz, &range, &bn, diag, subdiag, &VL, &VU, &IL, &IU, &abstol, &bm, eigs, V, &ldz, isup… in PetscDTSymmetricTridiagonalEigensolve()
1452 …ACKstegr_(&jobz, &range, &bn, diag, subdiag, &VL, &VU, &IL, &IU, &abstol, &bm, eigs, V, &ldz, isup… in PetscDTSymmetricTridiagonalEigensolve()
/petsc/src/ksp/pc/tutorials/
H A Dex4.c166 …PetscInt boundary_indices_size, am, an, bm, bn, condensed_am, astart, aend, Dstart, Dend,… in main() local
237 PetscCall(MatGetLocalSize(Bcondensed, &bm, &bn)); in main()
/petsc/src/mat/impls/nest/
H A Dmatnest.c1964 PetscInt bm, br; in MatAXPY_Dense_Nest() local
1980 PetscCall(ISGetLocalSize(nest->isglobal.row[i], &bm)); in MatAXPY_Dense_Nest()
1983 for (br = 0; br < bm; ++br) { in MatAXPY_Dense_Nest()
2123 PetscInt bm, *sub_dnnz, *sub_onnz, br; in MatConvert_Nest_AIJ() local
2126 PetscCall(ISGetLocalSize(nest->isglobal.row[i], &bm)); in MatConvert_Nest_AIJ()
2129 PetscCall(PetscMalloc1(bm, &iremote)); in MatConvert_Nest_AIJ()
2130 PetscCall(PetscMalloc1(bm, &sub_dnnz)); in MatConvert_Nest_AIJ()
2131 PetscCall(PetscMalloc1(bm, &sub_onnz)); in MatConvert_Nest_AIJ()
2132 for (k = 0; k < bm; ++k) { in MatConvert_Nest_AIJ()
2155 for (br = 0; br < bm; ++br) { in MatConvert_Nest_AIJ()
[all …]
/petsc/src/ksp/ksp/impls/gmres/dgmres/
H A Ddgmres.c848 PetscBLASInt NbrEig = 0, nr, bm; in KSPDGMRESImproveEig_DGMRES() local
960 PetscCall(PetscBLASIntCast(bmax, &bm)); in KSPDGMRESImproveEig_DGMRES()
963 PetscCallBLAS("LAPACKgetrf", LAPACKgetrf_(&nr, &nr, TTF, &bm, INVP, &info)); in KSPDGMRESImproveEig_DGMRES()
/petsc/src/dm/dt/dualspace/impls/lagrange/
H A Ddspacelagrange.c2575 PetscBLASInt bm, bn, bnrhs, blda, bldb, blwork, info; in PetscDualSpaceCreateInteriorSymmetryMatrix_Lagrange() local
2577 PetscCall(PetscBLASIntCast(nodeVecDim, &bm)); in PetscDualSpaceCreateInteriorSymmetryMatrix_Lagrange()
2580 PetscCall(PetscBLASIntCast(bm, &blda)); in PetscDualSpaceCreateInteriorSymmetryMatrix_Lagrange()
2581 PetscCall(PetscBLASIntCast(bm, &bldb)); in PetscDualSpaceCreateInteriorSymmetryMatrix_Lagrange()
2583 …PetscCallBLAS("LAPACKgels", LAPACKgels_(&transpose, &bm, &bn, &bnrhs, V, &blda, W, &bldb, work, &b… in PetscDualSpaceCreateInteriorSymmetryMatrix_Lagrange()
/petsc/src/mat/impls/dense/seq/
H A Ddense.c1010 PetscBLASInt o = 1, bm = 0; in MatSOR_SeqDense() local
1017 PetscCall(PetscBLASIntCast(m, &bm)); in MatSOR_SeqDense()
1029 PetscCallBLAS("BLASdotu", xt = b[i] - BLASdotu_(&bm, v + i, &bm, x, &o)); in MatSOR_SeqDense()
1035 PetscCallBLAS("BLASdotu", xt = b[i] - BLASdotu_(&bm, v + i, &bm, x, &o)); in MatSOR_SeqDense()
/petsc/src/mat/impls/sell/mpi/
H A Dmpisell.c136 …MatSeqXSELLReallocateSELL(B, bm, 1, nrow2, b->sliidx, b->sliceheight, row / sliceheight, row, col,…
/petsc/src/mat/impls/aij/mpi/
H A Dmpiaij.c478 …MatSeqXAIJReallocateAIJ(B, bm, 1, nrow2, row, col, rmax2, ba, bi, bj, rp2, ap2, bimax, nonew, MatS…
544 …PetscInt *bimax = b->imax, *bi = b->i, *bilen = b->ilen, *bj = b->j, bm = aij->B->rmap->n, am = … in MatSetValues_MPIAIJ() local
609 bm = aij->B->rmap->n; in MatSetValues_MPIAIJ()
8081 …PetscInt *bimax = b->imax, *bi = b->i, *bilen = b->ilen, *bj = b->j, bm = aij->B->rmap->n, am = … in matsetvaluesmpiaij_() local
8149 bm = aij->B->rmap->n; in matsetvaluesmpiaij_()