Lines Matching refs:detJ
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()
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()
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()
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()
1668 …metry_Internal(DM dm, PetscInt e, PetscReal v0[], PetscReal J[], PetscReal invJ[], PetscReal *detJ) in DMPlexComputePointGeometry_Internal() argument
1686 *detJ = 1.; in DMPlexComputePointGeometry_Internal()
1798 …metry_Internal(DM dm, PetscInt e, PetscReal v0[], PetscReal J[], PetscReal invJ[], PetscReal *detJ) in DMPlexComputeLineGeometry_Internal() argument
1808 *detJ = 0.0; in DMPlexComputeLineGeometry_Internal()
1828 DMPlex_Det3D_Internal(detJ, J); in DMPlexComputeLineGeometry_Internal()
1829 if (invJ) DMPlex_Invert3D_Internal(invJ, J, *detJ); in DMPlexComputeLineGeometry_Internal()
1845 DMPlex_Det2D_Internal(detJ, J); in DMPlexComputeLineGeometry_Internal()
1846 if (invJ) DMPlex_Invert2D_Internal(invJ, J, *detJ); in DMPlexComputeLineGeometry_Internal()
1856 *detJ = J[0]; in DMPlexComputeLineGeometry_Internal()
1868 …metry_Internal(DM dm, PetscInt e, PetscReal v0[], PetscReal J[], PetscReal invJ[], PetscReal *detJ) in DMPlexComputeTriangleGeometry_Internal() argument
1878 *detJ = 0.0; in DMPlexComputeTriangleGeometry_Internal()
1894 DMPlex_Det3D_Internal(detJ, J0); in DMPlexComputeTriangleGeometry_Internal()
1903 if (invJ) DMPlex_Invert3D_Internal(invJ, J, *detJ); in DMPlexComputeTriangleGeometry_Internal()
1915 DMPlex_Det2D_Internal(detJ, J); in DMPlexComputeTriangleGeometry_Internal()
1917 if (invJ) DMPlex_Invert2D_Internal(invJ, J, *detJ); in DMPlexComputeTriangleGeometry_Internal()
1923 …cInt Nq, const PetscReal points[], PetscReal v[], PetscReal J[], PetscReal invJ[], PetscReal *detJ) in DMPlexComputeRectangleGeometry_Internal() argument
1940 *detJ = 0.0; in DMPlexComputeRectangleGeometry_Internal()
1957 DMPlex_Det3D_Internal(detJ, J0); in DMPlexComputeRectangleGeometry_Internal()
1966 if (invJ) DMPlex_Invert3D_Internal(invJ, J, *detJ); in DMPlexComputeRectangleGeometry_Internal()
1979 DMPlex_Det2D_Internal(detJ, J); in DMPlexComputeRectangleGeometry_Internal()
1981 if (invJ) DMPlex_Invert2D_Internal(invJ, J, *detJ); in DMPlexComputeRectangleGeometry_Internal()
2065 DMPlex_Det3D_Internal(&detJ[i], &J[i * dim * dim]); in DMPlexComputeRectangleGeometry_Internal()
2066 if (invJ) DMPlex_Invert3D_Internal(&invJ[i * dim * dim], &J[i * dim * dim], detJ[i]); in DMPlexComputeRectangleGeometry_Internal()
2068 DMPlex_Det2D_Internal(&detJ[i], &J[i * dim * dim]); in DMPlexComputeRectangleGeometry_Internal()
2069 if (invJ) DMPlex_Invert2D_Internal(&invJ[i * dim * dim], &J[i * dim * dim], detJ[i]); in DMPlexComputeRectangleGeometry_Internal()
2078 …metry_Internal(DM dm, PetscInt e, PetscReal v0[], PetscReal J[], PetscReal invJ[], PetscReal *detJ) in DMPlexComputeTetrahedronGeometry_Internal() argument
2089 *detJ = 0.0; in DMPlexComputeTetrahedronGeometry_Internal()
2101 DMPlex_Det3D_Internal(detJ, J); in DMPlexComputeTetrahedronGeometry_Internal()
2103 if (invJ) DMPlex_Invert3D_Internal(invJ, J, *detJ); in DMPlexComputeTetrahedronGeometry_Internal()
2108 …cInt Nq, const PetscReal points[], PetscReal v[], PetscReal J[], PetscReal invJ[], PetscReal *detJ) in DMPlexComputeHexahedronGeometry_Internal() argument
2120 *detJ = 0.0; in DMPlexComputeHexahedronGeometry_Internal()
2131 DMPlex_Det3D_Internal(detJ, J); in DMPlexComputeHexahedronGeometry_Internal()
2133 if (invJ) DMPlex_Invert3D_Internal(invJ, J, *detJ); in DMPlexComputeHexahedronGeometry_Internal()
2215 DMPlex_Det3D_Internal(&detJ[i], &J[i * dim * dim]); in DMPlexComputeHexahedronGeometry_Internal()
2216 if (invJ) DMPlex_Invert3D_Internal(&invJ[i * dim * dim], &J[i * dim * dim], detJ[i]); in DMPlexComputeHexahedronGeometry_Internal()
2224 …cInt Nq, const PetscReal points[], PetscReal v[], PetscReal J[], PetscReal invJ[], PetscReal *detJ) in DMPlexComputeTriangularPrismGeometry_Internal() argument
2237 *detJ = 0.0; in DMPlexComputeTriangularPrismGeometry_Internal()
2248 DMPlex_Det3D_Internal(detJ, J); in DMPlexComputeTriangularPrismGeometry_Internal()
2250 if (invJ) DMPlex_Invert3D_Internal(invJ, J, *detJ); in DMPlexComputeTriangularPrismGeometry_Internal()
2356 DMPlex_Det3D_Internal(&detJ[i], &J[i * dim * dim]); in DMPlexComputeTriangularPrismGeometry_Internal()
2357 if (invJ) DMPlex_Invert3D_Internal(&invJ[i * dim * dim], &J[i * dim * dim], detJ[i]); in DMPlexComputeTriangularPrismGeometry_Internal()
2365 … PetscInt cell, PetscQuadrature quad, PetscReal *v, PetscReal *J, PetscReal *invJ, PetscReal *detJ) in DMPlexComputeCellGeometryFEM_Implicit() argument
2387 PetscCall(DMPlexComputePointGeometry_Internal(dm, cell, v, J, invJ, detJ)); in DMPlexComputeCellGeometryFEM_Implicit()
2393 else PetscCall(DMPlexComputeLineGeometry_Internal(dm, cell, v, J, invJ, detJ)); in DMPlexComputeCellGeometryFEM_Implicit()
2397 else PetscCall(DMPlexComputeTriangleGeometry_Internal(dm, cell, v, J, invJ, detJ)); in DMPlexComputeCellGeometryFEM_Implicit()
2400 …Call(DMPlexComputeRectangleGeometry_Internal(dm, cell, PETSC_FALSE, Nq, points, v, J, invJ, detJ)); in DMPlexComputeCellGeometryFEM_Implicit()
2404 …cCall(DMPlexComputeRectangleGeometry_Internal(dm, cell, PETSC_TRUE, Nq, points, v, J, invJ, detJ)); in DMPlexComputeCellGeometryFEM_Implicit()
2409 else PetscCall(DMPlexComputeTetrahedronGeometry_Internal(dm, cell, v, J, invJ, detJ)); in DMPlexComputeCellGeometryFEM_Implicit()
2412 PetscCall(DMPlexComputeHexahedronGeometry_Internal(dm, cell, Nq, points, v, J, invJ, detJ)); in DMPlexComputeCellGeometryFEM_Implicit()
2416 … PetscCall(DMPlexComputeTriangularPrismGeometry_Internal(dm, cell, Nq, points, v, J, invJ, detJ)); in DMPlexComputeCellGeometryFEM_Implicit()
2426 if (detJ) { in DMPlexComputeCellGeometryFEM_Implicit()
2427 for (i = 0; i < Nq; i++) detJ[i] = detJ0; in DMPlexComputeCellGeometryFEM_Implicit()
2482 …ryAffineFEM(DM dm, PetscInt cell, PetscReal v0[], PetscReal J[], PetscReal invJ[], PetscReal *detJ) in DMPlexComputeCellGeometryAffineFEM() argument
2485 PetscCall(DMPlexComputeCellGeometryFEM_Implicit(dm, cell, NULL, v0, J, invJ, detJ)); in DMPlexComputeCellGeometryAffineFEM()
2489 …scInt point, PetscQuadrature quad, PetscReal v[], PetscReal J[], PetscReal invJ[], PetscReal *detJ) in DMPlexComputeCellGeometryFEM_FE() argument
2561 if (!detJ && !invJ) continue; in DMPlexComputeCellGeometryFEM_FE()
2576 if (detJ) detJ[q] = detJt; in DMPlexComputeCellGeometryFEM_FE()
2578 …} else PetscCheck(!detJ && !invJ, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Need J to compute invJ or d… in DMPlexComputeCellGeometryFEM_FE()
2612 …scInt cell, PetscQuadrature quad, PetscReal v[], PetscReal J[], PetscReal invJ[], PetscReal detJ[]) in DMPlexComputeCellGeometryFEM() argument
2619 PetscAssertPointer(detJ, 7); in DMPlexComputeCellGeometryFEM()
2637 … (dim != cdim)) PetscCall(DMPlexComputeCellGeometryFEM_Implicit(dm, cell, quad, v, J, invJ, detJ)); in DMPlexComputeCellGeometryFEM()
2638 else PetscCall(DMPlexComputeCellGeometryFEM_FE(dm, fe, cell, quad, v, J, invJ, detJ)); in DMPlexComputeCellGeometryFEM()
3895 PetscReal detJ, *v0, *J, *invJ; in DMPlexCoordinatesToReference() local
3900 PetscCall(DMPlexComputeCellGeometryAffineFEM(dm, cell, v0, J, invJ, &detJ)); in DMPlexCoordinatesToReference()
3978 PetscReal detJ, *v0, *J; in DMPlexReferenceToCoordinates() local
3982 PetscCall(DMPlexComputeCellGeometryAffineFEM(dm, cell, v0, J, NULL, &detJ)); in DMPlexReferenceToCoordinates()