1504 …l, numComponents, &interpolant, coordDim * (Nq + 1), &coords, Nq, &fegeom.detJ, coordDim * coordDi…  in DMPlexComputeL2DiffLocal()
1512 …  PetscCall(DMPlexComputeCellGeometryFEM(dm, c, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ));  in DMPlexComputeL2DiffLocal()
1542         qgeom.detJ     = &fegeom.detJ[q];  in DMPlexComputeL2DiffLocal()
1543 …detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for element %" P…  in DMPlexComputeL2DiffLocal()
1555           PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMPlexComputeL2DiffLocal()
1565 …ble)(PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]), (double)PetscR…  in DMPlexComputeL2DiffLocal()
1566           elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q];  in DMPlexComputeL2DiffLocal()
1576   PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMPlexComputeL2DiffLocal()
1621 …coordDim * coordDim * Nq, &fegeom.invJ, numComponents * coordDim, &interpolant, Nq, &fegeom.detJ));  in DMComputeL2GradientDiff_Plex()
1628 …  PetscCall(DMPlexComputeCellGeometryFEM(dm, c, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ));  in DMComputeL2GradientDiff_Plex()
1651         qgeom.detJ     = &fegeom.detJ[q];  in DMComputeL2GradientDiff_Plex()
1652 …detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for element %" P…  in DMComputeL2GradientDiff_Plex()
1664           PetscCall(PetscFree6(funcVal, coords, fegeom.J, fegeom.invJ, interpolant, fegeom.detJ));  in DMComputeL2GradientDiff_Plex()
1679 …ield, fc, (double)(PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q])));  in DMComputeL2GradientDiff_Plex()
1680           elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q];  in DMComputeL2GradientDiff_Plex()
1690   PetscCall(PetscFree6(funcVal, coords, fegeom.J, fegeom.invJ, interpolant, fegeom.detJ));  in DMComputeL2GradientDiff_Plex()
1747 …totNc, &funcVal, totNc, &interpolant, dE * (Nq + 1), &coords, Nq, &fegeom.detJ, dE * dE * Nq, &feg…  in DMComputeL2FieldDiff_Plex()
1765 …cCall(DMPlexComputeCellGeometryFEM(dm, cone[0], quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ));  in DMComputeL2FieldDiff_Plex()
1767 …etscCall(DMPlexComputeCellGeometryFEM(dm, cell, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ));  in DMComputeL2FieldDiff_Plex()
1805           qgeom.detJ     = &fegeom.detJ[q];  in DMComputeL2FieldDiff_Plex()
1806 …detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for cell %" Pets…  in DMComputeL2FieldDiff_Plex()
1818             PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMComputeL2FieldDiff_Plex()
1829 …          (double)(PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q])));  in DMComputeL2FieldDiff_Plex()
1830 …          elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q];  in DMComputeL2FieldDiff_Plex()
1845     PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMComputeL2FieldDiff_Plex()
1916 …funcVal, numComponents, &interpolant, coordDim * Nq, &coords, Nq, &fegeom.detJ, coordDim * coordDi…  in DMPlexComputeL2DiffVec()
1923 …  PetscCall(DMPlexComputeCellGeometryFEM(dm, c, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ));  in DMPlexComputeL2DiffVec()
1948           qgeom.detJ     = &fegeom.detJ[q];  in DMPlexComputeL2DiffVec()
1949 …detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for element %" P…  in DMPlexComputeL2DiffVec()
1954             PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMPlexComputeL2DiffVec()
1963 …          elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q];  in DMPlexComputeL2DiffVec()
1973   PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMPlexComputeL2DiffVec()
2037 …(Nq + 1), &coords, cdim * cdim * Nq, &fegeom.J, cdim * cdim * Nq, &fegeom.invJ, Nq, &fegeom.detJ));  in DMPlexComputeL2FluxDiffVecLocal()
2044 …  PetscCall(DMPlexComputeCellGeometryFEM(dm, c, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ));  in DMPlexComputeL2FluxDiffVecLocal()
2054       qgeom.detJ     = &fegeom.detJ[q];  in DMPlexComputeL2FluxDiffVecLocal()
2056 …detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for element %" P…  in DMPlexComputeL2FluxDiffVecLocal()
2064 … fluxElemDiff += PetscSqr(PetscRealPart(interpolant[fc] - minterpolant[fc])) * wt * fegeom.detJ[q];  in DMPlexComputeL2FluxDiffVecLocal()
2072   PetscCall(PetscFree6(interpolant, minterpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMPlexComputeL2FluxDiffVecLocal()
2188 …PetscCall(PetscMalloc6(Nc * 2, &valsum, Nc, &interpolant, cdim * Nq, &coords, Nq, &fegeom.detJ, cd…  in DMPlexComputeClementInterpolant()
2206 …etscCall(DMPlexComputeCellGeometryFEM(dm, cell, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ));  in DMPlexComputeClementInterpolant()
2224           const PetscReal wt = quadWeights[q] * fegeom.detJ[q];  in DMPlexComputeClementInterpolant()
2230           qgeom.detJ     = &fegeom.detJ[q];  in DMPlexComputeClementInterpolant()
2231 …detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for element %" P…  in DMPlexComputeClementInterpolant()
2255   PetscCall(PetscFree6(valsum, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMPlexComputeClementInterpolant()
2321 …ordDim * numComponents, &interpolant, coordDim * Nq, &coords, Nq, &fegeom.detJ, coordDim * coordDi…  in DMPlexComputeGradientClementInterpolant()
2338 …etscCall(DMPlexComputeCellGeometryFEM(dm, cell, quad, coords, fegeom.J, fegeom.invJ, fegeom.detJ));  in DMPlexComputeGradientClementInterpolant()
2361           qgeom.detJ     = &fegeom.detJ[q];  in DMPlexComputeGradientClementInterpolant()
2362 …detJ[q] > 0.0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Invalid determinant %g for element %" P…  in DMPlexComputeGradientClementInterpolant()
2368 …= 0; d < coordDim; ++d) grad[fc * coordDim + d] += interpolant[fc * dim + d] * wt * fegeom.detJ[q];  in DMPlexComputeGradientClementInterpolant()
2370           vol += quadWeights[q * qNc] * fegeom.detJ[q];  in DMPlexComputeGradientClementInterpolant()
2397   PetscCall(PetscFree6(gradsum, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMPlexComputeGradientClementInterpolant()
3180   PetscReal   *x, *v0, *J, *invJ, detJ;  in DMPlexComputeInterpolatorGeneral()  local
3238         PetscCall(DMPlexComputeCellGeometryFEM(dmf, cell, NULL, v0, J, invJ, &detJ));  in DMPlexComputeInterpolatorGeneral()
3353   PetscReal   *x, *v0, *J, *invJ, detJ;  in DMPlexComputeMassMatrixGeneral()  local
3405       PetscCall(DMPlexComputeCellGeometryFEM(dmf, cell, NULL, v0, J, invJ, &detJ));  in DMPlexComputeMassMatrixGeneral()
3486       PetscCall(DMPlexComputeCellGeometryFEM(dmf, cell, NULL, v0, J, invJ, &detJ));  in DMPlexComputeMassMatrixGeneral()
3525 …0][j * Nc + c] * Tfine->T[0][(ccell * numFIndices + i) * Nc + c] * qweights[ccell * Nc + c] * detJ;  in DMPlexComputeMassMatrixGeneral()
3537               for (c = 0; c < Nc; ++c) elemMat[j] += 1.0 * 1.0 * qweights[ccell * Nc + c] * detJ;  in DMPlexComputeMassMatrixGeneral()