Home
last modified time | relevance | path

Searched refs:invJ (Results 1 – 25 of 35) sorted by relevance

12

/petsc/src/dm/dt/fe/interface/ !
H A Dfegeom.c40 PetscCall(PetscCalloc1(N * dimEmbed * dimEmbed, &g->invJ)); in PetscFEGeomCreate()
60 PetscCall(PetscFree((*geom)->invJ)); in PetscFEGeomDestroy()
104 (*chunkGeom)->invJ = PetscSafePointerPlusOffset(geom->invJ, Nq * dE * dE * cStart); in PetscFEGeomGetChunk()
175 pgeom->invJ = &geom->invJ[c * Np * dE * dE]; in PetscFEGeomGetPoint()
183 pgeom->invJ = &geom->invJ[(c * Np + p) * dE * dE]; in PetscFEGeomGetPoint()
228 pgeom->invJ = &geom->suppInvJ[0][c * Np * dE * dE]; in PetscFEGeomGetCellPoint()
232 pgeom->invJ = &geom->invJ[c * Np * dE * dE]; in PetscFEGeomGetCellPoint()
239 pgeom->invJ = &geom->suppInvJ[0][(c * Np + p) * dE * dE]; in PetscFEGeomGetCellPoint()
243 pgeom->invJ = &geom->invJ[(c * Np + p) * dE * dE]; in PetscFEGeomGetCellPoint()
271 …if (geom->invJ) DMPlex_Invert3D_Internal(&geom->invJ[dE * dE * i], &geom->J[dE * dE * i], geom->de… in PetscFEGeomComplete()
[all …]
/petsc/src/dm/impls/plex/ !
H A Dplexgeometry.c425 PetscReal v0, J, invJ, detJ; in DMPlexLocatePoint_Simplex_1D_Internal() local
429 PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, &v0, &J, &invJ, &detJ)); in DMPlexLocatePoint_Simplex_1D_Internal()
430 xi = invJ * (x - v0); in DMPlexLocatePoint_Simplex_1D_Internal()
441 PetscReal x[3], v0[3], J[9], invJ[9], detJ; in DMPlexLocatePoint_Simplex_2D_Internal() local
446 PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, v0, J, invJ, &detJ)); in DMPlexLocatePoint_Simplex_2D_Internal()
449 for (PetscInt j = 0; j < embedDim; ++j) xi[i] += invJ[i * embedDim + j] * (x[j] - v0[j]); in DMPlexLocatePoint_Simplex_2D_Internal()
461 PetscReal v0[2], J[4], invJ[4], detJ; in DMPlexClosestPoint_Simplex_2D_Internal() local
465 PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, v0, J, invJ, &detJ)); in DMPlexClosestPoint_Simplex_2D_Internal()
466 xi = invJ[0 * embedDim + 0] * (x - v0[0]) + invJ[0 * embedDim + 1] * (y - v0[1]); in DMPlexClosestPoint_Simplex_2D_Internal()
467 eta = invJ[1 * embedDim + 0] * (x - v0[0]) + invJ[1 * embedDim + 1] * (y - v0[1]); in DMPlexClosestPoint_Simplex_2D_Internal()
[all …]
H A Dplexproject.c277 fegeom.invJ = cgeom->invJ; in DMProjectPoint_Field_Private()
281 fgeomN[0].invJ = cgeom->suppInvJ[0]; in DMProjectPoint_Field_Private()
284 fgeomN[1].invJ = cgeom->suppInvJ[1]; in DMProjectPoint_Field_Private()
324 fegeom.invJ = &cgeom->invJ[tp * dE * dE]; in DMProjectPoint_Field_Private()
328 fgeomN[0].invJ = &cgeom->suppInvJ[0][tp * dE * dE]; in DMProjectPoint_Field_Private()
331 fgeomN[1].invJ = &cgeom->suppInvJ[1][tp * dE * dE]; in DMProjectPoint_Field_Private()
442 fegeom.invJ = fgeom->invJ; in DMProjectPoint_BdField_Private()
447 cgeom.invJ = fgeom->suppInvJ[0]; in DMProjectPoint_BdField_Private()
452 fgeomN[0].invJ = fgeom->suppInvJ[0]; in DMProjectPoint_BdField_Private()
455 fgeomN[1].invJ = fgeom->suppInvJ[1]; in DMProjectPoint_BdField_Private()
[all …]
H A Dplexfem.c1504 …s, Nq, &fegeom.detJ, coordDim * coordDim * Nq, &fegeom.J, coordDim * coordDim * Nq, &fegeom.invJ)); in DMPlexComputeL2DiffLocal()
1512 … PetscCall(DMPlexComputeCellGeometryFEM(dm, c, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ)); in DMPlexComputeL2DiffLocal()
1541 qgeom.invJ = &fegeom.invJ[q * coordDim * coordDim]; in DMPlexComputeL2DiffLocal()
1555 PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMPlexComputeL2DiffLocal()
1576 PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ)); in DMPlexComputeL2DiffLocal()
1621 …s, coordDim * coordDim * Nq, &fegeom.J, coordDim * coordDim * Nq, &fegeom.invJ, numComponents * co… in DMComputeL2GradientDiff_Plex()
1628 … PetscCall(DMPlexComputeCellGeometryFEM(dm, c, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ)); in DMComputeL2GradientDiff_Plex()
1650 qgeom.invJ = &fegeom.invJ[q * coordDim * coordDim]; in DMComputeL2GradientDiff_Plex()
1664 PetscCall(PetscFree6(funcVal, coords, fegeom.J, fegeom.invJ, interpolant, fegeom.detJ)); in DMComputeL2GradientDiff_Plex()
1690 PetscCall(PetscFree6(funcVal, coords, fegeom.J, fegeom.invJ, interpolant, fegeom.detJ)); in DMComputeL2GradientDiff_Plex()
[all …]
/petsc/include/ !
H A Dpetscfe.h49 …PetscReal *invJ; // invJ[Nc*Np*dE*dE]: The inverse of the Jacobian of the map from refer… member
188 const CeedScalar invJ[2][2] = { \
192 …const CeedScalar u_x[2] = {invJ[0][0] * du[i + Q * 0] + invJ[1][0] * du[i + Q * 1], invJ[0][1] * d…
201 dv[i + Q * 0] = qdata[i + Q * 0] * (invJ[0][0] * f1[0] + invJ[0][1] * f1[1]); \
202 dv[i + Q * 1] = qdata[i + Q * 0] * (invJ[1][0] * f1[0] + invJ[1][1] * f1[1]); \
H A Dpetsclandau.h65 void *invJ; // nip*dim*dim member
/petsc/src/dm/impls/plex/transform/impls/refine/regular/ !
H A Dplexrefregular.c45 …olytopeType ct, PetscInt *Nf, PetscReal *v0[], PetscReal *J[], PetscReal *invJ[], PetscReal *detJ[… in DMPlexRefineRegularGetAffineFaceTransforms() argument
97 if (invJ) *invJ = tri_invJ; in DMPlexRefineRegularGetAffineFaceTransforms()
104 if (invJ) *invJ = quad_invJ; in DMPlexRefineRegularGetAffineFaceTransforms()
133 …xTransform tr, DMPolytopeType ct, PetscInt *Nc, PetscReal *v0[], PetscReal *J[], PetscReal *invJ[]) in DMPlexRefineRegularGetAffineTransforms() argument
276 if (invJ) *invJ = seg_invJ; in DMPlexRefineRegularGetAffineTransforms()
282 if (invJ) *invJ = tri_invJ; in DMPlexRefineRegularGetAffineTransforms()
288 if (invJ) *invJ = quad_invJ; in DMPlexRefineRegularGetAffineTransforms()
294 if (invJ) *invJ = tet_invJ; in DMPlexRefineRegularGetAffineTransforms()
300 if (invJ) *invJ = hex_invJ; in DMPlexRefineRegularGetAffineTransforms()
/petsc/src/dm/impls/plex/tests/ !
H A Ddmplexcomputecellgeometryfem.F9035 …PetscReal :: J(tria_qpts*numDim**2), invJ(tria_qpts*numDim**2), v(tria_qpts*numDim), detJ(tria_qpt…
42 PetscCallA(DMPlexComputeCellGeometryFEM(dm, cStart, q, v, J, invJ, detJ, ierr))
62 …PetscReal :: J(quad_qpts*numDim**2), invJ(quad_qpts*numDim**2), v(quad_qpts*numDim), detJ(quad_qpt…
76 PetscCallA(DMPlexComputeCellGeometryFEM(dm, cStart, q, v, J, invJ, detJ, ierr))
H A Dex8.c17 PetscReal *v0, *J, *invJ, *detJ; /* FEM data */ member
56 …s->v0, numCells * dim * dim, &options->J, numCells * dim * dim, &options->invJ, numCells, &options… in ProcessOptions()
67 …tionsRealArray("-invJ", "Input inverse Jacobian for each cell", "ex8.c", options->invJ, &n, &flg)); in ProcessOptions()
74 PetscCall(PetscFree4(options->v0, options->J, options->invJ, options->detJ)); in ProcessOptions()
75 options->v0 = options->J = options->invJ = options->detJ = NULL; in ProcessOptions()
148 PetscReal v0[3], J[9], invJ[9], detJ; in CheckFEMGeometry() local
152 PetscCall(DMPlexComputeCellGeometryFEM(dm, cell, NULL, v0, J, invJ, &detJ)); in CheckFEMGeometry()
168 …invJ[i * spaceDim + j], invJEx[i * spaceDim + j]) < 10 * PETSC_SMALL, PETSC_COMM_SELF, PETSC_ERR_P… in CheckFEMGeometry()
636 PetscReal *invJ = PetscSafePointerPlusOffset(user.invJ, c * dim * dim); in main() local
642 …PetscCall(CheckCell(user.dm, c + cStart, PETSC_FALSE, v0, J, invJ, detJ, centroid, normal, vol, NU… in main()
[all …]
/petsc/src/dm/impls/plex/tutorials/ !
H A Dex3f90.F9027 PetscReal, target, dimension(dim*dim) :: invJ
48 pinvJ => invJ
/petsc/src/dm/field/impls/ds/ !
H A Ddmfieldds.c246 PetscReal *v, *J, *invJ, *detJ; in DMFieldEvaluate_DS() local
269 …herMax * dimR, &v, gatherMax * dimR * dimR, &J, gatherMax * dimR * dimR, &invJ, gatherMax, &detJ)); in DMFieldEvaluate_DS()
300 PetscCall(DMPlexComputeCellGeometryFEM(field->dm, c, quad, v, J, invJ, detJ)); in DMFieldEvaluate_DS()
325 …for (e = 0; e < dimR; e++) vs[d] += invJ[dimR * dimR * p + e * dimR + d] * cD[(nc * p + g) * dimR … in DMFieldEvaluate_DS()
338 …for (e = 0; e < dimR; e++) v[d] += invJ[dimR * dimR * p + e * dimR + d] * cD[(nc * p + g) * dimR +… in DMFieldEvaluate_DS()
356 …for (e = 0; e < dimR; e++) vs[d] += invJ[dimR * dimR * p + e * dimR + d] * cH[(nc * dimR * p + g) … in DMFieldEvaluate_DS()
366 …for (e = 0; e < dimR; e++) vs[d] += invJ[dimR * dimR * p + e * dimR + d] * cH[((nc * p + g) * dimR… in DMFieldEvaluate_DS()
380 …for (e = 0; e < dimR; e++) v[d] += invJ[dimR * dimR * p + e * dimR + d] * cH[(nc * dimR * p + g) *… in DMFieldEvaluate_DS()
388 …for (e = 0; e < dimR; e++) v[d] += invJ[dimR * dimR * p + e * dimR + d] * cH[((nc * p + g) * dimR … in DMFieldEvaluate_DS()
440 PetscCall(PetscFree4(v, J, invJ, detJ)); in DMFieldEvaluate_DS()
[all …]
/petsc/include/petsc/private/ !
H A Ddmpleximpl.h707 static inline void DMPlex_Invert2D_Internal(PetscReal invJ[], PetscReal J[], PetscReal detJ) in DMPlex_Invert2D_Internal() argument
712 invJ[0] = invDet * J[3]; in DMPlex_Invert2D_Internal()
713 invJ[1] = -invDet * J[1]; in DMPlex_Invert2D_Internal()
714 invJ[2] = -invDet * J[2]; in DMPlex_Invert2D_Internal()
715 invJ[3] = invDet * J[0]; in DMPlex_Invert2D_Internal()
719 static inline void DMPlex_Invert3D_Internal(PetscReal invJ[], PetscReal J[], PetscReal detJ) in DMPlex_Invert3D_Internal() argument
724 invJ[0 * 3 + 0] = invDet * (J[1 * 3 + 1] * J[2 * 3 + 2] - J[1 * 3 + 2] * J[2 * 3 + 1]); in DMPlex_Invert3D_Internal()
725 invJ[0 * 3 + 1] = invDet * (J[0 * 3 + 2] * J[2 * 3 + 1] - J[0 * 3 + 1] * J[2 * 3 + 2]); in DMPlex_Invert3D_Internal()
726 invJ[0 * 3 + 2] = invDet * (J[0 * 3 + 1] * J[1 * 3 + 2] - J[0 * 3 + 2] * J[1 * 3 + 1]); in DMPlex_Invert3D_Internal()
727 invJ[1 * 3 + 0] = invDet * (J[1 * 3 + 2] * J[2 * 3 + 0] - J[1 * 3 + 0] * J[2 * 3 + 2]); in DMPlex_Invert3D_Internal()
[all …]
H A Dpetscfeimpl.h286 …cInt dimRef, const PetscReal xi0[], const PetscReal v0[], const PetscReal invJ[], const PetscReal … in CoordinatesRealToRef() argument
292 for (e = 0; e < dimReal; ++e) xi[d] += invJ[d * dimReal + e] * (x[e] - v0[e]); in CoordinatesRealToRef()
369 …etscReal basisDer[], const PetscScalar x[], PetscInt dim, const PetscReal invJ[], const PetscReal … in PetscFEFreeInterpolateGradient_Static() argument
385 …for (g = 0; g < dim; ++g) realSpaceDer[d] += invJ[g * dim + d] * basisDer[((q * Nb + f) * Nc + fc)… in PetscFEFreeInterpolateGradient_Static()
/petsc/src/dm/impls/swarm/tests/ !
H A Dex2.c178 PetscReal *centroid, *coords, *xi0, *v0, *J, *invJ, detJ; in CreateParticles() local
216 PetscCall(PetscMalloc5(dim, &centroid, dim, &xi0, dim, &v0, dim * dim, &J, dim * dim, &invJ)); in CreateParticles()
223 PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, v0, J, invJ, &detJ)); /* affine */ in CreateParticles()
241 PetscCall(PetscFree5(centroid, xi0, v0, J, invJ)); in CreateParticles()
271 PetscReal *v0, *J, *invJ, detJ, *coords, *xi0; in CreateParticles_Shape() local
297 PetscCall(PetscMalloc4(dim, &xi0, dim, &v0, dim * dim, &J, dim * dim, &invJ)); in CreateParticles_Shape()
304 PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, v0, J, invJ, &detJ)); in CreateParticles_Shape()
314 PetscCall(PetscFree4(xi0, v0, J, invJ)); in CreateParticles_Shape()
482 …s, Nq, &fegeom.detJ, coordDim * coordDim * Nq, &fegeom.J, coordDim * coordDim * Nq, &fegeom.invJ)); in InterpolateGradient()
497 …PetscCall(DMPlexComputeCellGeometryFEM(dm, cell, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ)… in InterpolateGradient()
[all …]
/petsc/src/snes/utils/dm/ !
H A Ddminterpolatesnes.c373 PetscReal v0, J, invJ, detJ, xir[1]; in DMInterpolate_Segment_Private() local
379 PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, &v0, &J, &invJ, &detJ)); in DMInterpolate_Segment_Private()
381 xir[0] = invJ * PetscRealPart(coords[p] - v0); in DMInterpolate_Segment_Private()
400 PetscReal *v0, *J, *invJ, detJ; in DMInterpolate_Triangle_Private() local
404 PetscCall(PetscMalloc3(ctx->dim, &v0, ctx->dim * ctx->dim, &J, ctx->dim * ctx->dim, &invJ)); in DMInterpolate_Triangle_Private()
407 PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, v0, J, invJ, &detJ)); in DMInterpolate_Triangle_Private()
413 …for (PetscInt f = 0; f < ctx->dim; ++f) xi[d] += invJ[d * ctx->dim + f] * 0.5 * PetscRealPart(coor… in DMInterpolate_Triangle_Private()
419 PetscCall(PetscFree3(v0, J, invJ)); in DMInterpolate_Triangle_Private()
428 PetscReal *v0, *J, *invJ, detJ; in DMInterpolate_Tetrahedron_Private() local
434 PetscCall(PetscMalloc3(ctx->dim, &v0, ctx->dim * ctx->dim, &J, ctx->dim * ctx->dim, &invJ)); in DMInterpolate_Tetrahedron_Private()
[all …]
/petsc/src/dm/dt/fe/impls/basic/ !
H A Dfebasic.c218 fegeom.invJ = &cgeom->invJ[e * Np * dE * dE]; in PetscFEIntegrate_Basic()
230 fegeom.invJ = &cgeom->invJ[(e * Np + q) * dE * dE]; in PetscFEIntegrate_Basic()
238 PetscCall(DMPrintCellMatrix(e, "invJ", dim, dim, fegeom.invJ)); in PetscFEIntegrate_Basic()
310 fegeom.invJ = NULL; in PetscFEIntegrateBd_Basic()
320 fegeom.invJ = &fgeom->invJ[e * Np * dE * dE]; in PetscFEIntegrateBd_Basic()
325 cgeom.invJ = &fgeom->suppInvJ[0][e * Np * dE * dE]; in PetscFEIntegrateBd_Basic()
337 fegeom.invJ = &fgeom->invJ[(e * Np + q) * dE * dE]; in PetscFEIntegrateBd_Basic()
342 cgeom.invJ = &fgeom->suppInvJ[0][(e * Np + q) * dE * dE]; in PetscFEIntegrateBd_Basic()
350 PetscCall(DMPrintCellMatrix(e, "invJ", dim, dim, fegeom.invJ)); in PetscFEIntegrateBd_Basic()
455 PetscCall(DMPrintCellMatrix(e, "invJ", dE, dE, fegeom.invJ)); in PetscFEIntegrateResidual_Basic()
[all …]
/petsc/src/ts/utils/dmplexlandau/kokkos/ !
H A Dlandau.kokkos.cxx185 … invJ = new Kokkos::View<… in LandauKokkosStaticDataSet() local
186 …SData_d->invJ … in LandauKokkosStaticDataSet()
203 Kokkos::deep_copy(*invJ, h_invJ); in LandauKokkosStaticDataSet()
292 auto invJ = static_cast<Kokkos::View<PetscReal *, Kokkos::LayoutLeft> *>(SData_d->invJ); in LandauKokkosStaticDataClear() local
293 delete invJ; in LandauKokkosStaticDataClear()
409 … = *static_cast<Kokkos::View<PetscReal *, Kokkos::LayoutLeft> *>(SData_d->invJ); // use Kokkos vec… in LandauKokkosJacobian()
552 … const PetscReal *const invJ = &invJe[myQi * dim * dim]; // b_elem_idx: batch element index in LandauKokkosJacobian() local
565 …fdf_k(b_id, d + 1, idx) = 0.0; e < dim; ++e) d_fdf_k(b_id, d + 1, idx) += invJ[e * dim + d] * refS… in LandauKokkosJacobian()
606 const PetscReal *invJ = &d_invJ_k(jpidx_glb * dim * dim); in LandauKokkosJacobian() local
681 g2(d, fieldA, myQi) += invJ[d * dim + d2] * gg2(d2, fieldA, myQi); in LandauKokkosJacobian()
[all …]
/petsc/src/dm/impls/swarm/ !
H A Dswarmpic.c1011 PetscReal v0[3], J[9], invJ[9], detJ, detJp = 2. / (gmax[0] - gmin[0]), xr[3], den; in DMSwarmComputeLocalSize() local
1014 PetscCall(DMPlexComputeCellGeometryFEM(dm, cell, NULL, v0, J, invJ, &detJ)); in DMSwarmComputeLocalSize()
1167 PetscReal v0[3], J[9], invJ[9], detJ; in DMSwarmInitializeCoordinates() local
1171 PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, v0, J, invJ, &detJ)); in DMSwarmInitializeCoordinates()
1310 PetscReal *xi, *v0, *J, *invJ, detJ = 1.0, v0ref[3] = {-1.0, -1.0, -1.0}; in DMProjectFieldLocal_Swarm() local
1327 PetscCall(PetscMalloc3(dim, &v0, dim * dim, &J, dim * dim, &invJ)); in DMProjectFieldLocal_Swarm()
1335 fegeom.invJ = invJ; in DMProjectFieldLocal_Swarm()
1368 PetscCall(DMPlexComputeCellGeometryFEM(dmIn, cell, NULL, v0, J, invJ, &detJ)); in DMProjectFieldLocal_Swarm()
1384 CoordinatesRealToRef(dE, dim, fegeom.xi, fegeom.v, fegeom.invJ, xr, &xi[p * dim]); in DMProjectFieldLocal_Swarm()
1413 PetscCall(PetscFree3(v0, J, invJ)); in DMProjectFieldLocal_Swarm()
H A Dswarm.c477 PetscReal *xi, *v0, *J, *invJ, detJ = 1.0, v0ref[3] = {-1.0, -1.0, -1.0}; in DMSwarmComputeMassMatrix_Private() local
490 PetscCall(PetscMalloc3(dim, &v0, dim * dim, &J, dim * dim, &invJ)); in DMSwarmComputeMassMatrix_Private()
595 PetscCall(DMPlexComputeCellGeometryFEM(dmf, cell, NULL, v0, J, invJ, &detJ)); in DMSwarmComputeMassMatrix_Private()
606 CoordinatesRealToRef(dim, dim, v0ref, v0, invJ, xr, &xi[j * dim]); in DMSwarmComputeMassMatrix_Private()
634 PetscCall(PetscFree3(v0, J, invJ)); in DMSwarmComputeMassMatrix_Private()
700 PetscReal *xi, *v0, *J, *invJ, detJ = 1.0, v0ref[3] = {-1.0, -1.0, -1.0}; in DMSwarmComputeMassMatrixSquare_Private() local
713 PetscCall(PetscMalloc3(cdim, &v0, cdim * cdim, &J, cdim * cdim, &invJ)); in DMSwarmComputeMassMatrixSquare_Private()
821 PetscCall(DMPlexComputeCellGeometryFEM(dmf, cell, NULL, v0, J, invJ, &detJ)); in DMSwarmComputeMassMatrixSquare_Private()
832 CoordinatesRealToRef(cdim, cdim, v0ref, v0, invJ, xr, &xi[p * cdim]); in DMSwarmComputeMassMatrixSquare_Private()
867 PetscCall(PetscFree3(v0, J, invJ)); in DMSwarmComputeMassMatrixSquare_Private()
[all …]
/petsc/src/dm/dt/dualspace/interface/ !
H A Ddualspace.c1899 …PetscCall(PetscDTAltVPullbackMatrix(dim, dim, isInverse ? fegeom->J : fegeom->invJ, dsp->k, Jstar)… in PetscDualSpaceTransform()
1956 vals[(v * Nc + c) * dim] *= fegeom->invJ[0]; in PetscDualSpaceTransformGradient()
1959 …DMPlex_MultTranspose2DReal_Internal(fegeom->invJ, 1, &vals[(v * Nc + c) * dim], &vals[(v * Nc + c)… in PetscDualSpaceTransformGradient()
1962 …DMPlex_MultTranspose3DReal_Internal(fegeom->invJ, 1, &vals[(v * Nc + c) * dim], &vals[(v * Nc + c)… in PetscDualSpaceTransformGradient()
1971 …for (c = 0; c < Nc; ++c) DMPlex_MultTransposeReal_Internal(fegeom->invJ, dim, dE, 1, &vals[(v * Nc… in PetscDualSpaceTransformGradient()
2000 …DMPlex_MultTranspose2DReal_Internal(fegeom->invJ, dim, &vals[v * Nc * dim + d], &vals[v * Nc * dim… in PetscDualSpaceTransformGradient()
2003 …DMPlex_MultTranspose3DReal_Internal(fegeom->invJ, dim, &vals[v * Nc * dim + d], &vals[v * Nc * dim… in PetscDualSpaceTransformGradient()
2018 … DMPlex_Mult2DReal_Internal(fegeom->invJ, dim, &vals[v * Nc * dim + d], &vals[v * Nc * dim + d]); in PetscDualSpaceTransformGradient()
2021 … DMPlex_Mult3DReal_Internal(fegeom->invJ, dim, &vals[v * Nc * dim + d], &vals[v * Nc * dim + d]); in PetscDualSpaceTransformGradient()
2089 vals[(v * Nc + c) * dim * dim] *= PetscSqr(fegeom->invJ[0]); in PetscDualSpaceTransformHessian()
[all …]
/petsc/src/snes/tutorials/ !
H A Dex16.c282 void DeformationGradient(Field *ex, PetscInt qi, PetscInt qj, PetscInt qk, PetscScalar *invJ, Petsc… in DeformationGradient() argument
296 TensorVector(invJ, &grad[3 * bidx], lgrad); in DeformationGradient()
311 …tscInt qk, PetscInt ii, PetscInt jj, PetscInt kk, PetscInt fld, PetscScalar *invJ, PetscScalar *dF) in DeformationGradientJacobian() argument
319 TensorVector(invJ, &grad[3 * bidx], lgrad); in DeformationGradientJacobian()
465 PetscScalar invJ[9]; in FormElementJacobian() local
489 InvertTensor(J, invJ, &vol); in FormElementJacobian()
491 DeformationGradient(ex, qi, qj, qk, invJ, F); in FormElementJacobian()
502 TensorVector(invJ, &grad[3 * bidx], lgrad); in FormElementJacobian()
530 DeformationGradientJacobian(qi, qj, qk, i, j, k, l, invJ, dF); in FormElementJacobian()
542 TensorVector(invJ, &grad[3 * bidx], lgrad); in FormElementJacobian()
/petsc/src/ts/tests/ !
H A Dex27.c67 PetscReal *centroid, *coords, *velocity, *xi0, *v0, *J, *invJ, detJ, *vals; in SetInitialCoordinates() local
86 PetscCall(PetscMalloc5(dim, &centroid, dim, &xi0, dim, &v0, dim * dim, &J, dim * dim, &invJ)); in SetInitialCoordinates()
97 PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, v0, J, invJ, &detJ)); /* affine */ in SetInitialCoordinates()
127 PetscCall(PetscFree5(centroid, xi0, v0, J, invJ)); in SetInitialCoordinates()
H A Dex28.c64 PetscReal *centroid, *coords, *xi0, *v0, *J, *invJ, detJ, *vals; in SetInitialCoordinates() local
80 PetscCall(PetscMalloc5(dim, &centroid, dim, &xi0, dim, &v0, dim * dim, &J, dim * dim, &invJ)); in SetInitialCoordinates()
96 PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, v0, J, invJ, &detJ)); /* affine */ in SetInitialCoordinates()
114 PetscCall(PetscFree5(centroid, xi0, v0, J, invJ)); in SetInitialCoordinates()
/petsc/src/tao/tutorials/ !
H A Dex3.c123 PetscReal v0[3], J[9], invJ[9], detJ; in CreateMesh() local
125 PetscCall(DMPlexComputeCellGeometryFEM(*dm, 0, NULL, v0, J, invJ, &detJ)); in CreateMesh()
128 PetscCall(DMPlexComputeCellGeometryFEM(*dm, 0, NULL, v0, J, invJ, &detJ)); in CreateMesh()
/petsc/src/dm/dt/fe/impls/opencl/ !
H A Dfeopencl.c604 … for (d = 0; d < dim * dim; ++d) f_invJ[c * dim * dim + d] = (float)cgeom->invJ[c * dim * dim + d]; in PetscFEIntegrateResidual_OpenCL()
625 …for (d = 0; d < dim * dim; ++d) d_invJ[c * dim * dim + d] = (double)cgeom->invJ[c * dim * dim + d]; in PetscFEIntegrateResidual_OpenCL()
649 for (d = 0; d < dim * dim; ++d) r_invJ[c * dim * dim + d] = cgeom->invJ[c * dim * dim + d]; in PetscFEIntegrateResidual_OpenCL()

12