| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex2.c | 450 …PetscInt dim, coordDim, numFields, numComponents = 0, qNc, Nq, cStart, cEnd, vStart, vEnd,… in InterpolateGradient() local
456 PetscCall(DMGetCoordinateDim(dm, &coordDim)); in InterpolateGradient()
457 fegeom.dimEmbed = coordDim; in InterpolateGradient()
482 …coordDim * numComponents * 2, &gradsum, coordDim * numComponents, &interpolant, coordDim * Nq, &co… in InterpolateGradient()
489 PetscCall(PetscArrayzero(gradsum, coordDim * numComponents)); in InterpolateGradient()
493 PetscScalar *grad = &gradsum[coordDim * numComponents]; in InterpolateGradient()
504 PetscCall(PetscArrayzero(grad, coordDim * numComponents)); in InterpolateGradient()
521 …for (d = 0; d < coordDim; ++d) grad[fc * coordDim + d] += interpolant[fc * dim + d] * wt * fegeom.… in InterpolateGradient()
528 for (d = 0; d < coordDim; ++d) gradsum[fc * coordDim + d] += grad[fc * coordDim + d]; in InterpolateGradient()
533 for (d = 0; d < coordDim; ++d) { in InterpolateGradient()
[all …]
|
| /petsc/src/dm/impls/plex/cgns/ |
| H A D | plexcgns2.c | 476 int dim = 0, physDim = 0, coordDim = 0, numVertices = 0, numCells = 0; in DMPlexCreateCGNS_Internal_Serial() local
505 coordDim = dim; in DMPlexCreateCGNS_Internal_Serial()
509 PetscCallMPI(MPI_Bcast(&coordDim, 1, MPI_INT, 0, comm)); in DMPlexCreateCGNS_Internal_Serial()
620 PetscCall(DMSetCoordinateDim(*dm, coordDim)); in DMPlexCreateCGNS_Internal_Serial()
624 PetscCall(PetscSectionSetFieldComponents(coordSection, 0, coordDim)); in DMPlexCreateCGNS_Internal_Serial()
628 PetscCall(PetscSectionSetFieldDof(coordSection, v, 0, coordDim)); in DMPlexCreateCGNS_Internal_Serial()
652 …PetscCheck(ncoords == coordDim, PETSC_COMM_SELF, PETSC_ERR_LIB, "CGNS file must have a coordinate … in DMPlexCreateCGNS_Internal_Serial()
653 for (d = 0; d < coordDim; ++d) { in DMPlexCreateCGNS_Internal_Serial()
657 if (coordDim >= 1) { in DMPlexCreateCGNS_Internal_Serial()
658 for (v = 0; v < sizes[0]; ++v) coords[(v + off) * coordDim + 0] = x[0][v]; in DMPlexCreateCGNS_Internal_Serial()
[all …]
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexgmsh.c | 1580 PetscInt dim = 0, coordDim = -1, order = 0, maxHeight = 0; in DMPlexCreateGmsh() local
1603 …x_gmsh_spacedim", "Embedding space dimension", "DMPlexCreateGmsh", coordDim, &coordDim, NULL, PETS… in DMPlexCreateGmsh()
1968 if (coordDim < 0) coordDim = dim; in DMPlexCreateGmsh()
1969 PetscCall(DMSetCoordinateDim(*dm, coordDim)); in DMPlexCreateGmsh()
1978 PetscCall(GmshCreateFE(comm, NULL, isSimplex, continuity, nodeType, dim, coordDim, order, &fe)); in DMPlexCreateGmsh()
1991 PetscInt maxDof = GmshNumNodes_HEX(order) * coordDim; in DMPlexCreateGmsh()
2010 …for (d = 0; d < coordDim; ++d) cellCoords[(n + s) * coordDim + d] = (PetscReal)coords[node * 3 + d… in DMPlexCreateGmsh()
2031 for (d = 0; d < coordDim; ++d) cellCoords[n * coordDim + d] = 0.0; in DMPlexCreateGmsh()
2034 const PetscReal *weights = sdWeights[coordDim][n]; in DMPlexCreateGmsh()
2036 …for (d = 0; d < coordDim; ++d) cellCoords[n * coordDim + d] += weights[bn] * (PetscReal)coords[bno… in DMPlexCreateGmsh()
[all …]
|
| H A D | plexfem.c | 1469 …PetscInt dim, coordDim, numFields, numComponents = 0, qNc, Nq, cellHeight, cStart, cEnd, c… in DMPlexComputeL2DiffLocal() local
1474 PetscCall(DMGetCoordinateDim(dm, &coordDim)); in DMPlexComputeL2DiffLocal()
1475 fegeom.dimEmbed = coordDim; in DMPlexComputeL2DiffLocal()
1504 …mponents, &interpolant, coordDim * (Nq + 1), &coords, Nq, &fegeom.detJ, coordDim * coordDim * Nq, … in DMPlexComputeL2DiffLocal()
1540 qgeom.J = &fegeom.J[q * coordDim * coordDim]; in DMPlexComputeL2DiffLocal()
1541 qgeom.invJ = &fegeom.invJ[q * coordDim * coordDim]; in DMPlexComputeL2DiffLocal()
1545 gcoords = &coords[coordDim * Nq]; in DMPlexComputeL2DiffLocal()
1546 …MPlexBasisTransformApplyReal_Internal(dm, &coords[coordDim * q], PETSC_TRUE, coordDim, &coords[coo… in DMPlexComputeL2DiffLocal()
1548 gcoords = &coords[coordDim * q]; in DMPlexComputeL2DiffLocal()
1551 ierr = (*funcs[field])(coordDim, time, gcoords, Nc, funcVal, ctx); in DMPlexComputeL2DiffLocal()
[all …]
|
| H A D | plexproject.c | 73 PetscInt coordDim, Nf, *Nc, f, spDim, d, v, tp; in DMProjectPoint_Func_Private() local
77 PetscCall(DMGetCoordinateDim(dmIn, &coordDim)); in DMProjectPoint_Func_Private()
104 PetscCall(DMGetWorkArray(rdm, coordDim, MPIU_REAL, &x)); in DMProjectPoint_Func_Private()
119 … CoordinatesRefToReal(coordDim, fegeom->dim, fegeom->xi, fegeom->v, fegeom->J, refpoint, x); in DMProjectPoint_Func_Private()
122 v0 = &fegeom->v[tp * coordDim]; in DMProjectPoint_Func_Private()
125 …PetscCall(DMPlexBasisTransformApplyReal_Internal(dmIn, v0, PETSC_TRUE, coordDim, v0, x, dm->transf… in DMProjectPoint_Func_Private()
143 PetscCall((*funcs[f])(coordDim, time, v0, Nc[f], &pointEval[Nc[f] * q], ctx)); in DMProjectPoint_Func_Private()
148 PetscCall(DMRestoreWorkArray(rdm, coordDim, MPIU_REAL, &x)); in DMProjectPoint_Func_Private()
|
| H A D | plexmetric.c | 904 PetscInt coordDim, Nd; in DMPlexMetricCreate() local
907 PetscCall(DMGetCoordinateDim(dm, &coordDim)); in DMPlexMetricCreate()
908 Nd = coordDim * coordDim; in DMPlexMetricCreate()
|
| H A D | plexgeometry.c | 2368 PetscInt depth, dim, coordDim, coneSize, i; in DMPlexComputeCellGeometryFEM_Implicit() local
2381 PetscCall(DMGetCoordinateDim(dm, &coordDim)); in DMPlexComputeCellGeometryFEM_Implicit()
2382 …PetscCheck(coordDim <= 3, PetscObjectComm((PetscObject)dm), PETSC_ERR_SUP, "Unsupported coordinate… in DMPlexComputeCellGeometryFEM_Implicit()
2424 …for (i = 0; i < Nq; i++) CoordinatesRefToReal(coordDim, dim, xi0, v0, J0, &points[dim * i], &v[coo… in DMPlexComputeCellGeometryFEM_Implicit()
2435 for (j = 0; j < coordDim * coordDim; j++, k++) J[k] = J0[j]; in DMPlexComputeCellGeometryFEM_Implicit()
2440 switch (coordDim) { in DMPlexComputeCellGeometryFEM_Implicit()
2453 for (i = 1, k = coordDim * coordDim; i < Nq; i++) { in DMPlexComputeCellGeometryFEM_Implicit()
2456 for (j = 0; j < coordDim * coordDim; j++, k++) invJ[k] = invJ[j]; in DMPlexComputeCellGeometryFEM_Implicit()
|
| /petsc/src/dm/impls/forest/p4est/ |
| H A D | pforest.h | 201 PetscInt coordDim; member
211 PetscInt i, d = PetscMin(3, geom_pforest->coordDim); in GeometryMapping_pforest()
218 …ierr = (geom_pforest->map)(geom_pforest->base, (PetscInt)which_tree, geom_pforest->coordDim, Petsc… in GeometryMapping_pforest()
788 PetscCall(DMGetCoordinateDim(dm, &geom_pforest->coordDim)); in DMSetUp_pforest()
1764 PetscInt coordDim; in DMPlexCreateConnectivity_pforest() local
1766 PetscCall(DMGetCoordinateDim(dm, &coordDim)); in DMPlexCreateConnectivity_pforest()
1775 …coordDim, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Need coordinates at the corners: (dof) %" PetscIn… in DMPlexCreateConnectivity_pforest()
1777 PetscInt i, lim = PetscMin(3, coordDim); in DMPlexCreateConnectivity_pforest()
1783 …3 * (P4EST_CHILDREN * (c - cStart) + p4estVert) + i] = PetscRealPart(cellCoords[v * coordDim + i]); in DMPlexCreateConnectivity_pforest()
3592 PetscInt closureSize, c, coordStart, coordEnd, coordDim; in DMPforestMapCoordinates_Cell() local
[all …]
|
| /petsc/src/snes/utils/dm/ |
| H A D | dmadapt.c | 1316 PetscInt dim, coordDim, Nf; in DMAdaptorAdapt_Sequence_Private() local
1325 PetscCall(DMGetCoordinateDim(adaptor->idm, &coordDim)); in DMAdaptorAdapt_Sequence_Private()
1443 PetscInt Nd = coordDim * coordDim, f, vStart, vEnd; in DMAdaptorAdapt_Sequence_Private()
1474 …PetscCall(PetscFECreateDefault(PetscObjectComm((PetscObject)dmGrad), dim, Nc * coordDim, PETSC_TRU… in DMAdaptorAdapt_Sequence_Private()
|