442   PetscInt        embedDim;  in DMPlexLocatePoint_Simplex_2D_Internal()  local
445   PetscCall(DMGetCoordinateDim(dm, &embedDim));  in DMPlexLocatePoint_Simplex_2D_Internal()
447   for (PetscInt j = 0; j < embedDim; ++j) x[j] = PetscRealPart(point[j]);  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()
458   const PetscInt embedDim = 2;  in DMPlexClosestPoint_Simplex_2D_Internal()  local
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()
489   PetscInt           crossings = 0, numCoords, embedDim;  in DMPlexLocatePoint_Quad_2D_Linear_Internal()  local
494   embedDim = numCoords / 4;  in DMPlexLocatePoint_Quad_2D_Linear_Internal()
498     PetscReal x_i = PetscRealPart(coords[faces[2 * f + 0] * embedDim + 0]);  in DMPlexLocatePoint_Quad_2D_Linear_Internal()
499     PetscReal y_i = PetscRealPart(coords[faces[2 * f + 0] * embedDim + 1]);  in DMPlexLocatePoint_Quad_2D_Linear_Internal()
500     PetscReal x_j = PetscRealPart(coords[faces[2 * f + 1] * embedDim + 0]);  in DMPlexLocatePoint_Quad_2D_Linear_Internal()
501     PetscReal y_j = PetscRealPart(coords[faces[2 * f + 1] * embedDim + 1]);  in DMPlexLocatePoint_Quad_2D_Linear_Internal()
574   const PetscInt  embedDim = 3;  in DMPlexLocatePoint_Simplex_3D_Internal()  local
584 …xi   = invJ[0 * embedDim + 0] * (x - v0[0]) + invJ[0 * embedDim + 1] * (y - v0[1]) + invJ[0 * embe…  in DMPlexLocatePoint_Simplex_3D_Internal()
585 …eta  = invJ[1 * embedDim + 0] * (x - v0[0]) + invJ[1 * embedDim + 1] * (y - v0[1]) + invJ[1 * embe…  in DMPlexLocatePoint_Simplex_3D_Internal()
586 …zeta = invJ[2 * embedDim + 0] * (x - v0[0]) + invJ[2 * embedDim + 1] * (y - v0[1]) + invJ[2 * embe…  in DMPlexLocatePoint_Simplex_3D_Internal()