Lines Matching refs:pdim
1888 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()
2529 …PetscAssert(pdim == T->Nb, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Nb %" PetscInt_FMT " != %" PetscIn… in DMPlexComputeCellGeometryFEM_FE()
2537 for (k = 0; k < pdim; ++k) { in DMPlexComputeCellGeometryFEM_FE()
2539 …for (i = 0; i < cdim; ++i) v[q * cdim + i] += basis[(q * pdim + k) * cdim + i] * PetscRealPart(coo… in DMPlexComputeCellGeometryFEM_FE()
2541 PetscCall(PetscLogFlops(2.0 * pdim * cdim)); in DMPlexComputeCellGeometryFEM_FE()
2550 for (k = 0; k < pdim; ++k) { in DMPlexComputeCellGeometryFEM_FE()
2553 …for (i = 0; i < cdim; ++i) J[(q * cdim + i) * cdim + j] += basisDer[((q * pdim + k) * cdim + i) * … in DMPlexComputeCellGeometryFEM_FE()
2556 PetscCall(PetscLogFlops(2.0 * pdim * dim * cdim)); in DMPlexComputeCellGeometryFEM_FE()
3401 PetscInt dim, pdim, cStart, cEnd, cEndInterior, c; in DMPlexComputeGradientFVM() local
3405 PetscCall(PetscFVGetNumComponents(fvm, &pdim)); in DMPlexComputeGradientFVM()
3425 for (c = cStart; c < cEnd; ++c) PetscCall(PetscSectionSetDof(sectionGrad, c, pdim * dim)); in DMPlexComputeGradientFVM()
3729 PetscInt numComp, pdim, i, j, k, l, m, coordSize; in DMPlexCoordinatesToReference_FE() local
3737 PetscCall(PetscFEGetDimension(fe, &pdim)); in DMPlexCoordinatesToReference_FE()
3743 PetscCall(DMGetWorkArray(dm, pdim, MPIU_REAL, &modes)); in DMPlexCoordinatesToReference_FE()
3745 for (i = 0; i < pdim; ++i) { in DMPlexCoordinatesToReference_FE()
3747 for (j = 0; j < pdim; ++j) modes[i] += invV[i * pdim + j] * PetscRealPart(nodes[j]); in DMPlexCoordinatesToReference_FE()
3749 PetscCall(DMGetWorkArray(dm, pdim * Nc + pdim * Nc * dimR + Nc, MPIU_REAL, &B)); in DMPlexCoordinatesToReference_FE()
3750 D = &B[pdim * Nc]; in DMPlexCoordinatesToReference_FE()
3751 resNeg = &D[pdim * Nc * dimR]; in DMPlexCoordinatesToReference_FE()
3764 for (k = 0; k < pdim; k++) { in DMPlexCoordinatesToReference_FE()
3785 PetscCall(DMRestoreWorkArray(dm, pdim * Nc + pdim * Nc * dimR + Nc, MPIU_REAL, &B)); in DMPlexCoordinatesToReference_FE()
3786 PetscCall(DMRestoreWorkArray(dm, pdim, MPIU_REAL, &modes)); in DMPlexCoordinatesToReference_FE()
3794 PetscInt numComp, pdim, i, j, k, l, coordSize; in DMPlexReferenceToCoordinates_FE() local
3800 PetscCall(PetscFEGetDimension(fe, &pdim)); in DMPlexReferenceToCoordinates_FE()
3806 PetscCall(DMGetWorkArray(dm, pdim, MPIU_REAL, &modes)); in DMPlexReferenceToCoordinates_FE()
3808 for (i = 0; i < pdim; ++i) { in DMPlexReferenceToCoordinates_FE()
3810 for (j = 0; j < pdim; ++j) modes[i] += invV[i * pdim + j] * PetscRealPart(nodes[j]); in DMPlexReferenceToCoordinates_FE()
3812 PetscCall(DMGetWorkArray(dm, numPoints * pdim * Nc, MPIU_REAL, &B)); in DMPlexReferenceToCoordinates_FE()
3818 for (k = 0; k < pdim; k++) { in DMPlexReferenceToCoordinates_FE()
3819 for (l = 0; l < Nc; l++) mapped[l] += modes[k] * B[(j * pdim + k) * Nc + l]; in DMPlexReferenceToCoordinates_FE()
3822 PetscCall(DMRestoreWorkArray(dm, numPoints * pdim * Nc, MPIU_REAL, &B)); in DMPlexReferenceToCoordinates_FE()
3823 PetscCall(DMRestoreWorkArray(dm, pdim, MPIU_REAL, &modes)); in DMPlexReferenceToCoordinates_FE()