Home
last modified time | relevance | path

Searched refs:pdim (Results 1 – 14 of 14) sorted by relevance

/petsc/src/dm/dt/space/impls/point/
H A Dspacepoint.c66 PetscInt dim = sp->Nv, pdim = pt->quad->numPoints, d, p, i, c; in PetscSpaceEvaluate_Point() local
70 PetscCall(PetscArrayzero(B, npoints * pdim)); in PetscSpaceEvaluate_Point()
72 for (i = 0; i < pdim; ++i) { in PetscSpaceEvaluate_Point()
77 B[p * pdim + i] = 1.0; in PetscSpaceEvaluate_Point()
85 for (i = 0; i < pdim; ++i) B[(c * npoints + p) * pdim + i] = B[p * pdim + i]; in PetscSpaceEvaluate_Point()
88 if (D) PetscCall(PetscArrayzero(D, npoints * pdim * dim)); in PetscSpaceEvaluate_Point()
89 if (H) PetscCall(PetscArrayzero(H, npoints * pdim * dim * dim)); in PetscSpaceEvaluate_Point()
/petsc/src/dm/dt/space/impls/tensor/
H A Dspacetensor.c306 PetscInt pdim; in PetscSpaceEvaluate_Tensor() local
315 PetscCall(PetscSpaceGetDimension(sp, &pdim)); in PetscSpaceEvaluate_Tensor()
317 if (B || D || H) PetscCall(DMGetWorkArray(dm, npoints * pdim * Nc, MPIU_REAL, &sB)); in PetscSpaceEvaluate_Tensor()
318 if (D || H) PetscCall(DMGetWorkArray(dm, npoints * pdim * Nc * Nv, MPIU_REAL, &sD)); in PetscSpaceEvaluate_Tensor()
319 if (H) PetscCall(DMGetWorkArray(dm, npoints * pdim * Nc * Nv * Nv, MPIU_REAL, &sH)); in PetscSpaceEvaluate_Tensor()
321 for (PetscInt i = 0; i < npoints * pdim * Nc; i++) B[i] = 1.; in PetscSpaceEvaluate_Tensor()
324 for (PetscInt i = 0; i < npoints * pdim * Nc * Nv; i++) D[i] = 1.; in PetscSpaceEvaluate_Tensor()
327 for (PetscInt i = 0; i < npoints * pdim * Nc * Nv * Nv; i++) H[i] = 1.; in PetscSpaceEvaluate_Tensor()
336 …pdim % vstep) && !(pdim % spdim), PETSC_COMM_SELF, PETSC_ERR_PLIB, "Bad tensor loop: Nv %" PetscIn… in PetscSpaceEvaluate_Tensor()
338 vskip = pdim / (vstep * spdim); in PetscSpaceEvaluate_Tensor()
[all …]
/petsc/src/dm/dt/fe/impls/composite/
H A Dfecomposite.c26 PetscInt dim, pdim, spdim, j, s; in PetscFESetUp_Composite() local
39 PetscCall(PetscDualSpaceGetDimension(fem->dualSpace, &pdim)); in PetscFESetUp_Composite()
105 PetscInt pdim; /* Dimension of FE space P */ in PetscFEComputeTabulation_Composite() local
121 PetscCall(PetscDualSpaceGetDimension(fem->dualSpace, &pdim)); in PetscFEComputeTabulation_Composite()
150 if (K >= 0) PetscCall(PetscArrayzero(B, npoints * pdim * comp)); in PetscFEComputeTabulation_Composite()
151 if (K >= 1) PetscCall(PetscArrayzero(D, npoints * pdim * comp * dim)); in PetscFEComputeTabulation_Composite()
152 if (K >= 2) PetscCall(PetscArrayzero(H, npoints * pdim * comp * dim * dim)); in PetscFEComputeTabulation_Composite()
159 const PetscInt i = (p * pdim + cmp->embedding[s * spdim + j]) * comp; in PetscFEComputeTabulation_Composite()
169 const PetscInt i = ((p * pdim + cmp->embedding[s * spdim + j]) * comp + 0) * dim + d; in PetscFEComputeTabulation_Composite()
180 … const PetscInt i = ((p * pdim + cmp->embedding[s * spdim + j]) * comp + 0) * dim * dim + d; in PetscFEComputeTabulation_Composite()
/petsc/src/dm/impls/plex/
H A Dplexextrude.c105 PetscInt pNew, pdim = DMPolytopeTypeGetDim(rct[n]); in DMPlexExtrude() local
107 if (ct == rct[n] || pdim > dim) { in DMPlexExtrude()
109 PetscCall(DMLabelSetValue(cutLabel, pNew, !pdim ? 1 : 2)); in DMPlexExtrude()
H A Dplexgeometry.c1888 const PetscInt pdim = 2; in DMPlexComputeTriangleGeometry_Internal() local
1890 for (d = 0; d < pdim; d++) { in DMPlexComputeTriangleGeometry_Internal()
1891 … (PetscInt f = 0; f < pdim; f++) J0[d * dim + f] = 0.5 * (PetscRealPart(coords[(f + 1) * pdim + d]… in DMPlexComputeTriangleGeometry_Internal()
1950 const PetscInt pdim = 2; in DMPlexComputeRectangleGeometry_Internal() local
1952 for (d = 0; d < pdim; d++) { in DMPlexComputeRectangleGeometry_Internal()
1953 …dim + 0] = 0.5 * (PetscRealPart(coords[vorder[1] * pdim + d]) - PetscRealPart(coords[vorder[0] * p… in DMPlexComputeRectangleGeometry_Internal()
1954 …dim + 1] = 0.5 * (PetscRealPart(coords[vorder[2] * pdim + d]) - PetscRealPart(coords[vorder[1] * p… in DMPlexComputeRectangleGeometry_Internal()
2498 PetscInt dim, cdim, pdim, qdim, Nq, q; in DMPlexComputeCellGeometryFEM_FE() local
2514 PetscCall(PetscFEGetDimension(fe, &pdim)); in DMPlexComputeCellGeometryFEM_FE()
2518 …pdim * cdim, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "There are %" PetscInt_FMT " coordinates for poin… in DMPlexComputeCellGeometryFEM_FE()
[all …]
H A Dplexfem.c1047 PetscInt dim, nleaves, loc, fStart, fEnd, pdim, i; in DMPlexInsertBoundaryValuesRiemann() local
1069 PetscCall(PetscFVGetNumComponents(fv, &pdim)); in DMPlexInsertBoundaryValuesRiemann()
1070 PetscCall(DMGetWorkArray(dm, pdim, MPIU_SCALAR, &fx)); in DMPlexInsertBoundaryValuesRiemann()
1103 for (d = 0; d < pdim; ++d) fx[d] = cx[d] + DMPlex_DotD_Internal(dim, &cgrad[d * dim], dx); in DMPlexInsertBoundaryValuesRiemann()
1125 PetscCall(DMRestoreWorkArray(dm, pdim, MPIU_SCALAR, &fx)); in DMPlexInsertBoundaryValuesRiemann()
4563 PetscInt pdim, d; in DMPlexComputeResidual_Patch_Internal() local
4569 PetscCall(PetscFVGetNumComponents(fv, &pdim)); in DMPlexComputeResidual_Patch_Internal()
4582 for (d = 0; d < pdim; ++d) r[d] += u_t[d]; in DMPlexComputeResidual_Patch_Internal()
5373 PetscInt cdim, foff, pdim; in DMPlexComputeResidualByKey() local
5381 PetscCall(PetscFVGetNumComponents(fv, &pdim)); in DMPlexComputeResidualByKey()
[all …]
H A Dplexsubmesh.c3426 PetscInt closureSize, c, pdim; in DMPlexCreateSubmeshGeneric_Interpolated() local
3431 PetscCall(DMLabelGetValue(depth, closure[c], &pdim)); in DMPlexCreateSubmeshGeneric_Interpolated()
3432 PetscCall(DMLabelSetValue(subpointMap, closure[c], pdim)); in DMPlexCreateSubmeshGeneric_Interpolated()
/petsc/src/dm/dt/fv/interface/
H A Dfv.c1640 PetscInt pdim; // Dimension of approximation space P in PetscFVCreateTabulation() local
1655 pdim = Nc; in PetscFVCreateTabulation()
1660 (*T)->Nb = pdim; in PetscFVCreateTabulation()
1664 …for (k = 0; k <= (*T)->K; ++k) PetscCall(PetscMalloc1(nrepl * npoints * pdim * Nc * PetscPowInt(cd… in PetscFVCreateTabulation()
1667 for (d = 0; d < pdim; ++d) in PetscFVCreateTabulation()
1668 for (c = 0; c < Nc; ++c) (*T)->T[0][(p * pdim + d) * Nc + c] = 1.; in PetscFVCreateTabulation()
1672 for (d = 0; d < pdim; ++d) in PetscFVCreateTabulation()
1674 for (e = 0; e < cdim; ++e) (*T)->T[1][((p * pdim + d) * Nc + c) * cdim + e] = 0.0; in PetscFVCreateTabulation()
1678 for (d = 0; d < pdim; ++d) in PetscFVCreateTabulation()
1680 … for (e = 0; e < cdim * cdim; ++e) (*T)->T[2][((p * pdim + d) * Nc + c) * cdim * cdim + e] = 0.0; in PetscFVCreateTabulation()
[all …]
/petsc/src/dm/dt/fe/impls/basic/
H A Dfebasic.c51 PetscInt pdim, j; in PetscFESetUp_Basic() local
54 PetscCall(PetscDualSpaceGetDimension(fem->dualSpace, &pdim)); in PetscFESetUp_Basic()
55 PetscCall(PetscMalloc1(pdim * pdim, &fem->invV)); in PetscFESetUp_Basic()
56 for (j = 0; j < pdim; ++j) { in PetscFESetUp_Basic()
64 PetscCall(PetscMalloc1(Nc * Nq * pdim, &Bf)); in PetscFESetUp_Basic()
66 for (k = 0; k < pdim; ++k) { in PetscFESetUp_Basic()
68 fem->invV[j * pdim + k] = 0.0; in PetscFESetUp_Basic()
71 …for (c = 0; c < Nc; ++c) fem->invV[j * pdim + k] += Bf[(q * pdim + k) * Nc + c] * weights[q * Nc +… in PetscFESetUp_Basic()
77 PetscCall(PetscMalloc2(pdim, &pivots, pdim, &work)); in PetscFESetUp_Basic()
78 PetscCall(PetscBLASIntCast(pdim, &n)); in PetscFESetUp_Basic()
[all …]
/petsc/src/dm/dt/space/interface/
H A Dspace.c149 PetscInt pdim; in PetscSpaceView() local
156 PetscCall(PetscSpaceGetDimension(sp, &pdim)); in PetscSpaceView()
160 …T " variables with %" PetscInt_FMT " components, size %" PetscInt_FMT "\n", sp->Nv, sp->Nc, pdim)); in PetscSpaceView()
/petsc/src/mat/impls/fft/fftw/
H A Dfftw.c1233 ptrdiff_t *pdim; in MatCreate_FFTW() local
1246 pdim = (ptrdiff_t *)calloc(ndim, sizeof(ptrdiff_t)); in MatCreate_FFTW()
1247 pdim[0] = dim[0]; in MatCreate_FFTW()
1253 pdim[ctr] = dim[ctr]; in MatCreate_FFTW()
1316 fftw_mpi_local_size(ndim, pdim, comm, &local_n0, &local_0_start); in MatCreate_FFTW()
1321 temp = pdim[ndim - 1]; in MatCreate_FFTW()
1323 pdim[ndim - 1] = temp / 2 + 1; in MatCreate_FFTW()
1325 …fftw_mpi_local_size_transposed(ndim, pdim, comm, &local_n0, &local_0_start, &local_n1, &local_1_st… in MatCreate_FFTW()
1327 … PetscCall(PetscIntCast(2 * (PetscInt)local_n0 * partial_dim * pdim[ndim - 1] / temp, &fft->n)); in MatCreate_FFTW()
1328 N1 = 2 * fft->N * (PetscInt)pdim[ndim - 1] / ((PetscInt)temp); in MatCreate_FFTW()
[all …]
/petsc/src/dm/dt/dualspace/impls/lagrange/
H A Ddspacelagrange.c547 PetscInt dim, pdim, c, i, j, k, n, v, vStart, vEnd; in PetscLagNodeIndicesPushForward() local
568 pdim = pdepth != depth ? pdepth != 0 ? pdepth : 0 : dim; in PetscLagNodeIndicesPushForward()
628 for (j = 0; j < pdim; j++) J[i * pdim + j] = J[i * dim + j]; in PetscLagNodeIndicesPushForward()
633 PetscCall(PetscDTBinomialInt(pdim, PetscAbsInt(formDegree), &pNk)); in PetscLagNodeIndicesPushForward()
635 PetscCall(PetscDTAltVPullbackMatrix(pdim, dim, J, formDegree, Jstar)); in PetscLagNodeIndicesPushForward()
1514 PetscInt pdim, pNk; in PetscDualSpaceCreateAllDataFromInteriorData() local
1521 PetscCall(DMGetDimension(pdm, &pdim)); in PetscDualSpaceCreateAllDataFromInteriorData()
1522 if (pdim < PetscAbsInt(k)) continue; in PetscDualSpaceCreateAllDataFromInteriorData()
1523 PetscCall(PetscDTBinomialInt(pdim, PetscAbsInt(k), &pNk)); in PetscDualSpaceCreateAllDataFromInteriorData()
1551 PetscInt pdim, pNk; in PetscDualSpaceCreateAllDataFromInteriorData() local
[all …]
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A DDMDA.pyx112 cdef PetscInt pdim = PETSC_DECIDE
116 pdim = asDims(psizes, &m, &n, &p)
117 if gdim>=0 and pdim>=0:
118 assert gdim == pdim
350 cdef PetscInt pdim = PETSC_DECIDE
354 pdim = asDims(psizes, &m, &n, &p)
358 CHKERR(DMSetDimension(self.dm, pdim))
/petsc/src/dm/dt/dualspace/interface/
H A Ddualspace.c179 PetscInt pdim, f; in PetscDualSpaceView_ASCII() local
182 PetscCall(PetscDualSpaceGetDimension(sp, &pdim)); in PetscDualSpaceView_ASCII()
186 …" PetscInt_FMT "\n", PetscAbsInt(sp->k), sp->k < 0 ? "(stored in dual form) " : "", sp->Nc, pdim)); in PetscDualSpaceView_ASCII()
188 …Printf(v, "Dual space with %" PetscInt_FMT " components, size %" PetscInt_FMT "\n", sp->Nc, pdim)); in PetscDualSpaceView_ASCII()
194 for (f = 0; f < pdim; ++f) { in PetscDualSpaceView_ASCII()