19 static PetscErrorCode private_PetscFECreateDefault_scalar_pk1(DM dm, PetscInt dim, PetscBool isSimp…  in private_PetscFECreateDefault_scalar_pk1()  argument
38   PetscCall(PetscSpaceSetNumVariables(P, dim));  in private_PetscFECreateDefault_scalar_pk1()
46 …PetscCall(DMPlexCreateReferenceCell(PETSC_COMM_SELF, DMPolytopeTypeSimpleShape(dim, isSimplex), &K…  in private_PetscFECreateDefault_scalar_pk1()
70     PetscCall(PetscDTStroudConicalQuadrature(dim, 1, quadPointsPerEdge, -1.0, 1.0, &q));  in private_PetscFECreateDefault_scalar_pk1()
71     PetscCall(PetscDTStroudConicalQuadrature(dim - 1, 1, quadPointsPerEdge, -1.0, 1.0, &fq));  in private_PetscFECreateDefault_scalar_pk1()
73     PetscCall(PetscDTGaussTensorQuadrature(dim, 1, quadPointsPerEdge, -1.0, 1.0, &q));  in private_PetscFECreateDefault_scalar_pk1()
74     PetscCall(PetscDTGaussTensorQuadrature(dim - 1, 1, quadPointsPerEdge, -1.0, 1.0, &fq));  in private_PetscFECreateDefault_scalar_pk1()
85   PetscInt         dim, nfaces, nbasis;  in private_DMSwarmInsertPointsUsingCellDM_PLEX_SubDivide()  local
101   PetscCall(DMGetDimension(dmc, &dim));  in private_DMSwarmInsertPointsUsingCellDM_PLEX_SubDivide()
105   if (nfaces == (dim + 1)) is_simplex = PETSC_TRUE;  in private_DMSwarmInsertPointsUsingCellDM_PLEX_SubDivide()
107   PetscCall(private_PetscFECreateDefault_scalar_pk1(dmc, dim, is_simplex, 0, &fe));  in private_DMSwarmInsertPointsUsingCellDM_PLEX_SubDivide()
141       for (d = 0; d < dim; d++) {  in private_DMSwarmInsertPointsUsingCellDM_PLEX_SubDivide()
142         swarm_coor[dim * pcnt + d] = 0.0;  in private_DMSwarmInsertPointsUsingCellDM_PLEX_SubDivide()
143 …for (k = 0; k < nbasis; k++) swarm_coor[dim * pcnt + d] += T->T[0][q * nbasis + k] * PetscRealPart…  in private_DMSwarmInsertPointsUsingCellDM_PLEX_SubDivide()
159   PetscInt      dim;  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()  local
173   PetscCall(DMGetDimension(dmc, &dim));  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
174   PetscCheck(dim == 2, PetscObjectComm((PetscObject)dm), PETSC_ERR_SUP, "Only 2D is supported");  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
178   if (nfaces == (dim + 1)) is_simplex = PETSC_TRUE;  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
181   PetscCall(PetscMalloc1(dim * npoints * npoints, &xi));  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
187       xi[dim * pcnt + 0] = ii * ds;  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
188       xi[dim * pcnt + 1] = jj * ds;  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
190       xi[dim * pcnt + 0] *= (1.0 - 1.2 * ds2);  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
191       xi[dim * pcnt + 1] *= (1.0 - 1.2 * ds2);  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
193       xi[dim * pcnt + 0] += 0.35 * ds2;  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
194       xi[dim * pcnt + 1] += 0.35 * ds2;  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
205     basis[q][0] = 1.0 - xi[dim * q + 0] - xi[dim * q + 1];  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
206     basis[q][1] = xi[dim * q + 0];  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
207     basis[q][2] = xi[dim * q + 1];  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
231       for (d = 0; d < dim; d++) {  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
232         swarm_coor[dim * pcnt + d] = 0.0;  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
233 …for (k = 0; k < npe; k++) swarm_coor[dim * pcnt + d] += basis[q][k] * PetscRealPart(elcoor[dim * k…  in private_DMSwarmInsertPointsUsingCellDM_PLEX2D_Regular()
251   PetscInt dim;  in private_DMSwarmInsertPointsUsingCellDM_PLEX()  local
254   PetscCall(DMGetDimension(celldm, &dim));  in private_DMSwarmInsertPointsUsingCellDM_PLEX()
257 …PetscCheck(dim != 3, PetscObjectComm((PetscObject)dm), PETSC_ERR_SUP, "No 3D support for REGULAR+P…  in private_DMSwarmInsertPointsUsingCellDM_PLEX()
284   PetscInt      dim, ps, pe, nel, nfaces, Nfc;  in private_DMSwarmSetPointCoordinatesCellwise_PLEX()  local
291   PetscCall(DMGetDimension(dmc, &dim));  in private_DMSwarmSetPointCoordinatesCellwise_PLEX()
298   if (nfaces == (dim + 1)) is_simplex = PETSC_TRUE;  in private_DMSwarmSetPointCoordinatesCellwise_PLEX()
300   switch (dim) {  in private_DMSwarmSetPointCoordinatesCellwise_PLEX()
315 …for (PetscInt d = 0; d < dim; d++) PetscCheck(xi[dim * p + d] >= -1.0, PetscObjectComm((PetscObjec…  in private_DMSwarmSetPointCoordinatesCellwise_PLEX()
317       for (PetscInt d = 0; d < dim; d++) sum += xi[dim * p + d];  in private_DMSwarmSetPointCoordinatesCellwise_PLEX()
322 …for (PetscInt d = 0; d < dim; d++) PetscCheck(PetscAbsReal(xi[dim * p + d]) <= 1.0, PetscObjectCom…  in private_DMSwarmSetPointCoordinatesCellwise_PLEX()
340 … PetscCall(DMPlexReferenceToCoordinates(dmc, e + ps, npoints, xi, &swarm_coor[npoints * dim * e]));  in private_DMSwarmSetPointCoordinatesCellwise_PLEX()