Lines Matching refs:v0
425 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()
476 cpoint[0] = J[0 * embedDim + 0] * xi + J[0 * embedDim + 1] * eta + v0[0]; in DMPlexClosestPoint_Simplex_2D_Internal()
477 cpoint[1] = J[1 * embedDim + 0] * xi + J[1 * embedDim + 1] * eta + v0[1]; in DMPlexClosestPoint_Simplex_2D_Internal()
576 PetscReal v0[3], J[9], invJ[9], detJ; in DMPlexLocatePoint_Simplex_3D_Internal() local
583 PetscCall(DMPlexComputeCellGeometryFEM(dm, c, NULL, v0, J, invJ, &detJ)); in DMPlexLocatePoint_Simplex_3D_Internal()
584 …i = invJ[0 * embedDim + 0] * (x - v0[0]) + invJ[0 * embedDim + 1] * (y - v0[1]) + invJ[0 * embed… in DMPlexLocatePoint_Simplex_3D_Internal()
585 …ta = invJ[1 * embedDim + 0] * (x - v0[0]) + invJ[1 * embedDim + 1] * (y - v0[1]) + invJ[1 * embed… in DMPlexLocatePoint_Simplex_3D_Internal()
586 …eta = invJ[2 * embedDim + 0] * (x - v0[0]) + invJ[2 * embedDim + 1] * (y - v0[1]) + invJ[2 * embed… in DMPlexLocatePoint_Simplex_3D_Internal()
1668 static PetscErrorCode DMPlexComputePointGeometry_Internal(DM dm, PetscInt e, PetscReal v0[], PetscR… in DMPlexComputePointGeometry_Internal() argument
1682 if (v0) { in DMPlexComputePointGeometry_Internal()
1683 for (d = 0; d < dim; d++) v0[d] = PetscRealPart(coords[off + d]); in DMPlexComputePointGeometry_Internal()
1798 static PetscErrorCode DMPlexComputeLineGeometry_Internal(DM dm, PetscInt e, PetscReal v0[], PetscRe… in DMPlexComputeLineGeometry_Internal() argument
1813 if (v0) { in DMPlexComputeLineGeometry_Internal()
1814 for (d = 0; d < dim; d++) v0[d] = PetscRealPart(coords[d]); in DMPlexComputeLineGeometry_Internal()
1835 if (v0) { in DMPlexComputeLineGeometry_Internal()
1836 for (d = 0; d < dim; d++) v0[d] = PetscRealPart(coords[d]); in DMPlexComputeLineGeometry_Internal()
1851 if (v0) { in DMPlexComputeLineGeometry_Internal()
1852 for (d = 0; d < dim; d++) v0[d] = PetscRealPart(coords[d]); in DMPlexComputeLineGeometry_Internal()
1868 static PetscErrorCode DMPlexComputeTriangleGeometry_Internal(DM dm, PetscInt e, PetscReal v0[], Pet… in DMPlexComputeTriangleGeometry_Internal() argument
1883 if (v0) { in DMPlexComputeTriangleGeometry_Internal()
1884 for (d = 0; d < dim; d++) v0[d] = PetscRealPart(coords[d]); in DMPlexComputeTriangleGeometry_Internal()
1907 if (v0) { in DMPlexComputeTriangleGeometry_Internal()
1908 for (d = 0; d < dim; d++) v0[d] = PetscRealPart(coords[d]); in DMPlexComputeTriangleGeometry_Internal()
2078 static PetscErrorCode DMPlexComputeTetrahedronGeometry_Internal(DM dm, PetscInt e, PetscReal v0[], … in DMPlexComputeTetrahedronGeometry_Internal() argument
2090 if (v0) { in DMPlexComputeTetrahedronGeometry_Internal()
2091 for (d = 0; d < dim; d++) v0[d] = PetscRealPart(coords[d]); in DMPlexComputeTetrahedronGeometry_Internal()
2372 PetscReal xi0[3] = {-1., -1., -1.}, v0[3], J0[9], detJ0; in DMPlexComputeCellGeometryFEM_Implicit() local
2392 if (Nq) PetscCall(DMPlexComputeLineGeometry_Internal(dm, cell, v0, J0, NULL, &detJ0)); in DMPlexComputeCellGeometryFEM_Implicit()
2396 if (Nq) PetscCall(DMPlexComputeTriangleGeometry_Internal(dm, cell, v0, J0, NULL, &detJ0)); in DMPlexComputeCellGeometryFEM_Implicit()
2408 if (Nq) PetscCall(DMPlexComputeTetrahedronGeometry_Internal(dm, cell, v0, J0, NULL, &detJ0)); in DMPlexComputeCellGeometryFEM_Implicit()
2424 …for (i = 0; i < Nq; i++) CoordinatesRefToReal(coordDim, dim, xi0, v0, J0, &points[dim * i], &v[coo… in DMPlexComputeCellGeometryFEM_Implicit()
2482 PetscErrorCode DMPlexComputeCellGeometryAffineFEM(DM dm, PetscInt cell, PetscReal v0[], PetscReal J… in DMPlexComputeCellGeometryAffineFEM() argument
2485 PetscCall(DMPlexComputeCellGeometryFEM_Implicit(dm, cell, NULL, v0, J, invJ, detJ)); in DMPlexComputeCellGeometryAffineFEM()
3895 PetscReal detJ, *v0, *J, *invJ; in DMPlexCoordinatesToReference() local
3897 PetscCall(DMGetWorkArray(dm, dimC + 2 * dimC * dimC, MPIU_REAL, &v0)); in DMPlexCoordinatesToReference()
3898 J = &v0[dimC]; in DMPlexCoordinatesToReference()
3900 PetscCall(DMPlexComputeCellGeometryAffineFEM(dm, cell, v0, J, invJ, &detJ)); in DMPlexCoordinatesToReference()
3904 CoordinatesRealToRef(dimC, dimR, x0, v0, invJ, &realCoords[dimC * i], &refCoords[dimR * i]); in DMPlexCoordinatesToReference()
3906 PetscCall(DMRestoreWorkArray(dm, dimC + 2 * dimC * dimC, MPIU_REAL, &v0)); in DMPlexCoordinatesToReference()
3978 PetscReal detJ, *v0, *J; in DMPlexReferenceToCoordinates() local
3980 PetscCall(DMGetWorkArray(dm, dimC + 2 * dimC * dimC, MPIU_REAL, &v0)); in DMPlexReferenceToCoordinates()
3981 J = &v0[dimC]; in DMPlexReferenceToCoordinates()
3982 PetscCall(DMPlexComputeCellGeometryAffineFEM(dm, cell, v0, J, NULL, &detJ)); in DMPlexReferenceToCoordinates()
3986 CoordinatesRefToReal(dimC, dimR, xi0, v0, J, &refCoords[dimR * i], &realCoords[dimC * i]); in DMPlexReferenceToCoordinates()
3988 PetscCall(DMRestoreWorkArray(dm, dimC + 2 * dimC * dimC, MPIU_REAL, &v0)); in DMPlexReferenceToCoordinates()