| /petsc/src/dm/dt/fe/interface/ ! |
| H A D | fegeom.c | 88 PetscInt dE; in PetscFEGeomGetChunk() local
95 dE = geom->dimEmbed; in PetscFEGeomGetChunk()
102 (*chunkGeom)->v = PetscSafePointerPlusOffset(geom->v, Nq * dE * cStart); in PetscFEGeomGetChunk()
103 (*chunkGeom)->J = PetscSafePointerPlusOffset(geom->J, Nq * dE * dE * cStart); in PetscFEGeomGetChunk()
104 (*chunkGeom)->invJ = PetscSafePointerPlusOffset(geom->invJ, Nq * dE * dE * cStart); in PetscFEGeomGetChunk()
106 (*chunkGeom)->n = PetscSafePointerPlusOffset(geom->n, Nq * dE * cStart); in PetscFEGeomGetChunk()
108 (*chunkGeom)->suppJ[0] = PetscSafePointerPlusOffset(geom->suppJ[0], Nq * dE * dE * cStart); in PetscFEGeomGetChunk()
109 (*chunkGeom)->suppJ[1] = PetscSafePointerPlusOffset(geom->suppJ[1], Nq * dE * dE * cStart); in PetscFEGeomGetChunk()
110 (*chunkGeom)->suppInvJ[0] = PetscSafePointerPlusOffset(geom->suppInvJ[0], Nq * dE * dE * cStart); in PetscFEGeomGetChunk()
111 (*chunkGeom)->suppInvJ[1] = PetscSafePointerPlusOffset(geom->suppInvJ[1], Nq * dE * dE * cStart); in PetscFEGeomGetChunk()
[all …]
|
| H A D | fe.c | 2367 const PetscInt dE = fegeom->dimEmbed; in PetscFEEvaluateFieldJets_Internal() local
2378 for (d = 0; d < dE * Ncf; ++d) u_x[fOffset * dE + d] = 0.0; in PetscFEEvaluateFieldJets_Internal()
2384 …for (d = 0; d < cdim; ++d) u_x[(fOffset + c) * dE + d] += Dq[cidx * cdim + d] * coefficients[dOffs… in PetscFEEvaluateFieldJets_Internal()
2388 for (g = 0; g < Nf; ++g) hOffset += T[g]->Nc * dE; in PetscFEEvaluateFieldJets_Internal()
2389 for (d = 0; d < dE * dE * Ncf; ++d) u_x[hOffset + fOffset * dE * dE + d] = 0.0; in PetscFEEvaluateFieldJets_Internal()
2394 …for (d = 0; d < cdim * cdim; ++d) u_x[hOffset + (fOffset + c) * dE * dE + d] += Hq[cidx * cdim * c… in PetscFEEvaluateFieldJets_Internal()
2397 PetscCall(PetscFEPushforwardHessian(fe, fegeom, 1, &u_x[hOffset + fOffset * dE * dE])); in PetscFEEvaluateFieldJets_Internal()
2400 PetscCall(PetscFEPushforwardGradient(fe, fegeom, 1, &u_x[fOffset * dE])); in PetscFEEvaluateFieldJets_Internal()
2434 const PetscInt dE = fegeom->dimEmbed; in PetscFEEvaluateFieldJets_Hybrid_Internal() local
2449 for (d = 0; d < dE * Ncf; ++d) u_x[fOffset * dE + d] = 0.0; in PetscFEEvaluateFieldJets_Hybrid_Internal()
[all …]
|
| /petsc/src/dm/interface/ ! |
| H A D | dmgeommodel.c | 10 static PetscErrorCode DMSnapToGeomModelBall(DM dm, PetscInt p, PetscInt dE, const PetscScalar mcoor… in DMSnapToGeomModelBall() argument
19 for (PetscInt d = 0; d < dE; ++d) norm += PetscSqr(PetscRealPart(mcoords[d])); in DMSnapToGeomModelBall()
21 for (PetscInt d = 0; d < dE; ++d) gcoords[d] = mcoords[d] / norm; in DMSnapToGeomModelBall()
23 for (PetscInt d = 0; d < dE; ++d) gcoords[d] = mcoords[d]; in DMSnapToGeomModelBall()
28 static PetscErrorCode DMSnapToGeomModelCylinder(DM dm, PetscInt p, PetscInt dE, const PetscScalar m… in DMSnapToGeomModelCylinder() argument
39 for (PetscInt d = 0; d < dE - 1; ++d) norm += PetscSqr(PetscRealPart(mcoords[d])); in DMSnapToGeomModelCylinder()
41 for (PetscInt d = 0; d < dE - 1; ++d) gcoords[d] = mcoords[d] * gmax[0] / norm; in DMSnapToGeomModelCylinder()
42 gcoords[dE - 1] = mcoords[dE - 1]; in DMSnapToGeomModelCylinder()
44 for (PetscInt d = 0; d < dE; ++d) gcoords[d] = mcoords[d]; in DMSnapToGeomModelCylinder()
179 PetscErrorCode DMSnapToGeomModel(DM dm, PetscInt p, PetscInt dE, const PetscScalar mcoords[], Petsc… in DMSnapToGeomModel() argument
[all …]
|
| H A D | dmcoordinates.c | 1013 PetscInt dim, dE, height, cStart, cEnd, gct; in DMCreateAffineCoordinates_Internal() local
1023 PetscCall(DMGetCoordinateDim(dm, &dE)); in DMCreateAffineCoordinates_Internal()
1034 PetscCall(PetscFECreateLagrangeByCell(PETSC_COMM_SELF, dim, dE, ct, 1, -1, &feLinear)); in DMCreateAffineCoordinates_Internal()
|
| H A D | dm.c | 5938 PetscInt dE, Nf = dm->Nf, f, s, Nl, l, Ndef, k; in DMCreateDS() local
5945 PetscCall(DMGetCoordinateDim(dm, &dE)); in DMCreateDS()
6003 if (dsDef) PetscCall(PetscDSSetCoordinateDimension(dsDef, dE)); in DMCreateDS()
6047 PetscCall(PetscDSSetCoordinateDimension(dsDef, dE)); in DMCreateDS()
6073 PetscCall(PetscDSSetCoordinateDimension(ds, dE)); in DMCreateDS()
6096 PetscCall(PetscDSSetCoordinateDimension(dsIn, dE)); in DMCreateDS()
|
| /petsc/src/dm/dt/fe/impls/basic/ ! |
| H A D | febasic.c | 173 …PetscInt dim, dE, Np, numConstants, Nf, NfAux = 0, totDim, totDimAux = 0, cOffset = 0, c… in PetscFEIntegrate_Basic() local
206 dE = cgeom->dimEmbed; in PetscFEIntegrate_Basic()
217 fegeom.J = &cgeom->J[e * Np * dE * dE]; in PetscFEIntegrate_Basic()
218 fegeom.invJ = &cgeom->invJ[e * Np * dE * dE]; in PetscFEIntegrate_Basic()
226 …CoordinatesRefToReal(dE, dim, fegeom.xi, &cgeom->v[e * Np * dE], fegeom.J, &quadPoints[q * dim], x… in PetscFEIntegrate_Basic()
228 fegeom.v = &cgeom->v[(e * Np + q) * dE]; in PetscFEIntegrate_Basic()
229 fegeom.J = &cgeom->J[(e * Np + q) * dE * dE]; in PetscFEIntegrate_Basic()
230 fegeom.invJ = &cgeom->invJ[(e * Np + q) * dE * dE]; in PetscFEIntegrate_Basic()
267 PetscInt qNc, Nq, q, Np, dE; in PetscFEIntegrateBd_Basic() local
301 dE = fgeom->dimEmbed; in PetscFEIntegrateBd_Basic()
[all …]
|
| /petsc/src/dm/field/interface/ ! |
| H A D | dmfield.c | 442 PetscInt dim, dE; in DMFieldCreateFEGeom() local
452 dE = field->numComponents; in DMFieldCreateFEGeom()
453 PetscCall(PetscFEGeomCreate(quad, nPoints, dE, mode, &g)); in DMFieldCreateFEGeom()
456 if (dE > dim) { in DMFieldCreateFEGeom()
463 for (j = 0; j < dE; j++) { in DMFieldCreateFEGeom()
464 for (k = 0; k < dim; k++) J[j * dE + k] = g->J[i * dE * dim + j * dim + k]; in DMFieldCreateFEGeom()
468 for (j = 0; j < dE; j++) { in DMFieldCreateFEGeom()
469 for (k = 0; k < dE; k++) J[j * dE + k] = (j == k) ? 1. : 0.; in DMFieldCreateFEGeom()
473 if (dE == 2) { in DMFieldCreateFEGeom()
519 for (j = 0; j < dE * dE; j++) g->J[i * dE * dE + j] = J[j]; in DMFieldCreateFEGeom()
|
| /petsc/src/dm/impls/plex/ ! |
| H A D | plexproject.c | 208 const PetscInt dE = cgeom->dimEmbed, *cone, *ornt; in DMProjectPoint_Field_Private() local
320 CoordinatesRefToReal(dE, cgeom->dim, fegeom.xi, cgeom->v, fegeom.J, &points[q * dim], x); in DMProjectPoint_Field_Private()
322 fegeom.v = &cgeom->v[tp * dE]; in DMProjectPoint_Field_Private()
323 fegeom.J = &cgeom->J[tp * dE * dE]; in DMProjectPoint_Field_Private()
324 fegeom.invJ = &cgeom->invJ[tp * dE * dE]; in DMProjectPoint_Field_Private()
327 fgeomN[0].J = &cgeom->suppJ[0][tp * dE * dE]; in DMProjectPoint_Field_Private()
328 fgeomN[0].invJ = &cgeom->suppInvJ[0][tp * dE * dE]; in DMProjectPoint_Field_Private()
330 fgeomN[1].J = &cgeom->suppJ[1][tp * dE * dE]; in DMProjectPoint_Field_Private()
331 fgeomN[1].invJ = &cgeom->suppInvJ[1][tp * dE * dE]; in DMProjectPoint_Field_Private()
340 …tscCall(DMPlexBasisTransformApplyReal_Internal(dmIn, fegeom.v, PETSC_TRUE, dE, fegeom.v, fegeom.v,… in DMProjectPoint_Field_Private()
[all …]
|
| H A D | plexegads.c | 90 …ds[], ego obj, const PetscScalar range[], const PetscInt v, const PetscInt dE, PetscScalar paramsV… in DMPlex_EGADS_EDGE_XYZtoUV_Internal() argument
112 dx = coords[v * dE + 0] - eval[0]; in DMPlex_EGADS_EDGE_XYZtoUV_Internal()
113 dy = coords[v * dE + 1] - eval[1]; in DMPlex_EGADS_EDGE_XYZtoUV_Internal()
114 dz = coords[v * dE + 2] - eval[2]; in DMPlex_EGADS_EDGE_XYZtoUV_Internal()
145 …dE + 0] - data[0]) * (coords[v * dE + 0] - data[0]) + (coords[v * dE + 1] - data[1]) * (coords[v *… in DMPlex_EGADS_EDGE_XYZtoUV_Internal()
179 …ds[], ego obj, const PetscScalar range[], const PetscInt v, const PetscInt dE, PetscScalar paramsV… in DMPlex_Geom_EDGE_XYZtoUV_Internal() argument
200 dx = coords[v * dE + 0] - eval[0]; in DMPlex_Geom_EDGE_XYZtoUV_Internal()
201 dy = coords[v * dE + 1] - eval[1]; in DMPlex_Geom_EDGE_XYZtoUV_Internal()
202 dz = coords[v * dE + 2] - eval[2]; in DMPlex_Geom_EDGE_XYZtoUV_Internal()
237 …dE + 0] - data[0]) * (coords[v * dE + 0] - data[0]) + (coords[v * dE + 1] - data[1]) * (coords[v *… in DMPlex_Geom_EDGE_XYZtoUV_Internal()
[all …]
|
| H A D | plexfem.c | 1705 PetscInt dim, depth, dE, Nf, f, Nds, s; in DMComputeL2FieldDiff_Plex() local
1710 PetscCall(DMGetCoordinateDim(dm, &dE)); in DMComputeL2FieldDiff_Plex()
1747 …(totNc, &funcVal, totNc, &interpolant, dE * (Nq + 1), &coords, Nq, &fegeom.detJ, dE * dE * Nq, &fe… in DMComputeL2FieldDiff_Plex()
1803 qgeom.J = &fegeom.J[q * dE * dE]; in DMComputeL2FieldDiff_Plex()
1804 qgeom.invJ = &fegeom.invJ[q * dE * dE]; in DMComputeL2FieldDiff_Plex()
1808 gcoords = &coords[dE * Nq]; in DMComputeL2FieldDiff_Plex()
1809 …tscCall(DMPlexBasisTransformApplyReal_Internal(dm, &coords[dE * q], PETSC_TRUE, dE, &coords[dE * q… in DMComputeL2FieldDiff_Plex()
1811 gcoords = &coords[dE * q]; in DMComputeL2FieldDiff_Plex()
1814 ierr = (*funcs[fields[f]])(dE, time, gcoords, Nc, funcVal, ctx); in DMComputeL2FieldDiff_Plex()
1820 …if (transform) PetscCall(DMPlexBasisTransformApply_Internal(dm, &coords[dE * q], PETSC_FALSE, Nc, … in DMComputeL2FieldDiff_Plex()
[all …]
|
| H A D | plexgeometry.c | 4192 PetscInt dE, d, e; in DMPlexShearGeometry() local
4196 PetscCall(DMGetCoordinateDim(dm, &dE)); in DMPlexShearGeometry()
4197 PetscCall(PetscMalloc1(dE + 1, &moduli)); in DMPlexShearGeometry()
4199 …for (d = 0, e = 0; d < dE; ++d) moduli[d + 1] = d == dir ? 0.0 : (multipliers ? multipliers[e++] :… in DMPlexShearGeometry()
4201 PetscCall(PetscDSSetConstants(cds, dE + 1, moduli)); in DMPlexShearGeometry()
|
| H A D | plexcreate.c | 304 PetscInt dim, cdim, dE, qorder, height; in DMPlexCreateCoordinateSpace() local
309 PetscCall(DMGetCoordinateDim(dm, &dE)); in DMPlexCreateCoordinateSpace()
389 dE = cdim; in DMPlexCreateCoordinateSpace()
404 PetscCall(PetscFECreateLagrangeByCell(PETSC_COMM_SELF, dim, dE, ct, degree, qorder, &fe)); in DMPlexCreateCoordinateSpace()
|
| H A D | plex.c | 9653 PetscInt dim, depth, dE, d, cStart, cEnd, c; in DMPlexCheckGeometry() local
9657 PetscCall(DMGetCoordinateDim(dm, &dE)); in DMPlexCheckGeometry()
9658 if (dim != dE) PetscFunctionReturn(PETSC_SUCCESS); in DMPlexCheckGeometry()
|
| /petsc/src/dm/field/impls/ds/ ! |
| H A D | dmfieldds.c | 789 PetscInt p, dim, dE, numFaces, Nq; in DMFieldComputeFaceData_DS() local
797 dE = geom->dimEmbed; in DMFieldComputeFaceData_DS()
829 PetscInt Np = geom->numPoints, q, dE = geom->dimEmbed, d; in DMFieldComputeFaceData_DS() local
832 for (d = 0; d < dE; ++d) geom->n[(p * Np + q) * dE + d] = -geom->n[(p * Np + q) * dE + d]; in DMFieldComputeFaceData_DS()
892 for (q = 0; q < Nq * dE * dE; q++) { in DMFieldComputeFaceData_DS()
893 geom->suppJ[s][p * Nq * dE * dE + q] = cellGeom->J[offset * Nq * dE * dE + q]; in DMFieldComputeFaceData_DS()
894 geom->suppInvJ[s][p * Nq * dE * dE + q] = cellGeom->invJ[offset * Nq * dE * dE + q]; in DMFieldComputeFaceData_DS()
935 PetscCall(PetscMalloc1(dE * Nq, &cellPoints)); in DMFieldComputeFaceData_DS()
938 PetscCall(PetscQuadratureSetData(cellQuad, dE, 1, Nq, cellPoints, unusedWeights)); in DMFieldComputeFaceData_DS()
1086 … CoordinatesRefToReal(dE, dim, xi0, v0, J, &orientPoints[o][dim * q + 0], &cellPoints[dE * q + 0]); in DMFieldComputeFaceData_DS()
[all …]
|
| /petsc/src/dm/dt/fe/tests/ ! |
| H A D | ex3.c | 206 PetscInt dE, n, i; in DistortMesh() local
209 PetscCall(DMGetCoordinateDim(dm, &dE)); in DistortMesh()
213 for (i = 0; i < (n / dE); ++i) { in DistortMesh()
214 ca[i * dE + 0] += user->shear * ca[i * dE + 0]; in DistortMesh()
215 ca[i * dE + 1] *= user->flatten; in DistortMesh()
|
| /petsc/src/dm/dt/dualspace/interface/ ! |
| H A D | dualspace.c | 1144 PetscInt dim, dE, qNc, c, Nq, q; in PetscDualSpaceApplyDefault() local
1158 dE = cgeom->dimEmbed; in PetscDualSpaceApplyDefault()
1161 CoordinatesRefToReal(dE, cgeom->dim, cgeom->xi, cgeom->v, cgeom->J, &points[q * dim], x); in PetscDualSpaceApplyDefault()
1162 PetscCall((*func)(dE, time, x, Nc, val, ctx)); in PetscDualSpaceApplyDefault()
1164 PetscCall((*func)(dE, time, &cgeom->v[dE * q], Nc, val, ctx)); in PetscDualSpaceApplyDefault()
1942 const PetscInt dim = dsp->dm->dim, dE = fegeom->dimEmbed; in PetscDualSpaceTransformGradient() local
1949 …PetscAssert(dE > 0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid embedding dimension %" Pet… in PetscDualSpaceTransformGradient()
1951 if (dim == dE) { in PetscDualSpaceTransformGradient()
1971 …x_MultTransposeReal_Internal(fegeom->invJ, dim, dE, 1, &vals[(v * Nc + c) * dE], &vals[(v * Nc + c… in PetscDualSpaceTransformGradient()
2075 const PetscInt dim = dsp->dm->dim, dE = fegeom->dimEmbed; in PetscDualSpaceTransformHessian() local
[all …]
|
| /petsc/src/dm/impls/plex/transform/interface/ ! |
| H A D | plextransform.c | 1730 PetscInt v, dE = trdim, d, off = 0; in DMPlexTransformCreateCellVertices_Internal() local
1735 …for (d = 0; d < dE; ++d) PetscCall(PetscPrintf(PETSC_COMM_SELF, "%g ", (double)PetscRealPart(tr->t… in DMPlexTransformCreateCellVertices_Internal()
1807 …Type pct, DMPolytopeType ct, PetscInt p, PetscInt r, PetscInt Nv, PetscInt dE, const PetscScalar i… in DMPlexTransformMapCoordinatesBarycenter_Internal() argument
1813 for (d = 0; d < dE; ++d) out[d] = 0.0; in DMPlexTransformMapCoordinatesBarycenter_Internal()
1815 for (d = 0; d < dE; ++d) out[d] += in[v * dE + d]; in DMPlexTransformMapCoordinatesBarycenter_Internal()
1816 for (d = 0; d < dE; ++d) out[d] /= Nv; in DMPlexTransformMapCoordinatesBarycenter_Internal()
1842 …Type pct, DMPolytopeType ct, PetscInt p, PetscInt r, PetscInt Nv, PetscInt dE, const PetscScalar i… in DMPlexTransformMapCoordinates() argument
1845 if (Nv) PetscUseTypeMethod(tr, mapcoordinates, pct, ct, p, r, Nv, dE, in, out); in DMPlexTransformMapCoordinates()
2170 …PetscInt dE, dEo, d, cStart, cEnd, c, cStartNew, cEndNew, vStartNew, vEndNew, v, pStart,… in DMPlexTransformSetCoordinates() local
2195 PetscCall(DMGetCoordinateDim(rdm, &dE)); in DMPlexTransformSetCoordinates()
[all …]
|
| /petsc/src/dm/impls/plex/transform/impls/extrude/ ! |
| H A D | plextrcohesive.c | 988 …Type pct, DMPolytopeType ct, PetscInt p, PetscInt r, PetscInt Nv, PetscInt dE, const PetscScalar i… in DMPlexTransformMapCoordinates_Cohesive() argument
1002 for (PetscInt d = 0; d < dE; ++d) out[d] = in[d]; in DMPlexTransformMapCoordinates_Cohesive()
1029 for (PetscInt d = 0; d < dE; ++d) avgNormal[d] += normal[d]; in DMPlexTransformMapCoordinates_Cohesive()
1033 for (PetscInt d = 0; d < dE; ++d) norm += PetscSqr(avgNormal[d]); in DMPlexTransformMapCoordinates_Cohesive()
1035 for (PetscInt d = 0; d < dE; ++d) avgNormal[d] /= norm; in DMPlexTransformMapCoordinates_Cohesive()
1037 for (PetscInt d = 0; d < dE; ++d) out[d] = in[d] + width * avgNormal[d] * (r ? -0.5 : 0.5); in DMPlexTransformMapCoordinates_Cohesive()
|
| H A D | plextrextrude.c | 787 …Type pct, DMPolytopeType ct, PetscInt p, PetscInt r, PetscInt Nv, PetscInt dE, const PetscScalar i… in DMPlexTransformMapCoordinates_Extrude() argument
803 …etscCheck(dE == ex->cdim, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Coordinate dim %" PetscInt_FMT "… in DMPlexTransformMapCoordinates_Extrude()
|
| /petsc/src/snes/tutorials/ ! |
| H A D | ex16.c | 354 PetscScalar FtdF[9], dE[9]; in SaintVenantKirchoffJacobian() local
357 TensorTransposeTensor(dF, F, dE); in SaintVenantKirchoffJacobian()
359 for (i = 0; i < 9; i++) dE[i] += FtdF[i]; in SaintVenantKirchoffJacobian()
360 for (i = 0; i < 9; i++) dE[i] *= 0.5; in SaintVenantKirchoffJacobian()
361 for (i = 0; i < 3; i++) dtrE += dE[i + 3 * i]; in SaintVenantKirchoffJacobian()
364 dS[i + 3 * j] = 2. * mu * dE[i + 3 * j]; in SaintVenantKirchoffJacobian()
|
| /petsc/src/dm/impls/plex/transform/impls/refine/bl/ ! |
| H A D | plexrefbl.c | 489 …Type pct, DMPolytopeType ct, PetscInt p, PetscInt r, PetscInt Nv, PetscInt dE, const PetscScalar i… in DMPlexTransformMapCoordinates_BL() argument
500 for (d = 0; d < dE; ++d) out[d] = in[d] + bl->h[r] * (in[d + dE] - in[d]); in DMPlexTransformMapCoordinates_BL()
503 PetscCall(DMPlexTransformMapCoordinatesBarycenter_Internal(tr, pct, ct, p, r, Nv, dE, in, out)); in DMPlexTransformMapCoordinates_BL()
|
| /petsc/src/dm/impls/swarm/ ! |
| H A D | swarmpic.c | 1311 PetscInt dim, dE, Np, n, Nf, Nfc, Nfu, cStart, cEnd, maxC = 0, totbs = 0; in DMProjectFieldLocal_Swarm() local
1320 PetscCall(DMGetCoordinateDim(dmIn, &dE)); in DMProjectFieldLocal_Swarm()
1331 fegeom.dimEmbed = dE; in DMProjectFieldLocal_Swarm()
1384 CoordinatesRealToRef(dE, dim, fegeom.xi, fegeom.v, fegeom.invJ, xr, &xi[p * dim]); in DMProjectFieldLocal_Swarm()
|