1465   PetscScalar     *funcVal, *interpolant;  in DMPlexComputeL2DiffLocal()  local
1504 …PetscCall(PetscMalloc6(numComponents, &funcVal, numComponents, &interpolant, coordDim * (Nq + 1), …  in DMPlexComputeL2DiffLocal()
1555           PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMPlexComputeL2DiffLocal()
1558 …ASSID) PetscCall(PetscFEInterpolate_Static((PetscFE)obj, &x[fieldOffset], &qgeom, q, interpolant));  in DMPlexComputeL2DiffLocal()
1559 …TSCFV_CLASSID) PetscCall(PetscFVInterpolate_Static((PetscFV)obj, &x[fieldOffset], q, interpolant));  in DMPlexComputeL2DiffLocal()
1565 …(double)(PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]), (double)Pe…  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()
1589   PetscScalar     *funcVal, *interpolant;  in DMComputeL2GradientDiff_Plex()  local
1621 ….J, coordDim * coordDim * Nq, &fegeom.invJ, numComponents * coordDim, &interpolant, Nq, &fegeom.de…  in DMComputeL2GradientDiff_Plex()
1664           PetscCall(PetscFree6(funcVal, coords, fegeom.J, fegeom.invJ, interpolant, fegeom.detJ));  in DMComputeL2GradientDiff_Plex()
1667 …      PetscCall(PetscFEInterpolateGradient_Static(fe, 1, &x[fieldOffset], &qgeom, q, interpolant));  in DMComputeL2GradientDiff_Plex()
1673             for (d = 0; d < dim; ++d) sum += interpolant[fc * dim + d] * n[d];  in DMComputeL2GradientDiff_Plex()
1674             interpolant[fc] = sum;  in DMComputeL2GradientDiff_Plex()
1679 …tscInt_FMT " diff %g\n", c, field, fc, (double)(PetscSqr(PetscRealPart(interpolant[fc] - funcVal[f…  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()
1735     PetscScalar     *funcVal, *interpolant;  in DMComputeL2FieldDiff_Plex()  local
1747 …PetscCall(PetscCalloc6(totNc, &funcVal, totNc, &interpolant, dE * (Nq + 1), &coords, Nq, &fegeom.d…  in DMComputeL2FieldDiff_Plex()
1818             PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMComputeL2FieldDiff_Plex()
1822 …SCFE_CLASSID) PetscCall(PetscFEInterpolate_Static((PetscFE)obj, &x[fOff], &qgeom, q, interpolant));  in DMComputeL2FieldDiff_Plex()
1823 …d == PETSCFV_CLASSID) PetscCall(PetscFVInterpolate_Static((PetscFV)obj, &x[fOff], q, interpolant));  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()
1879   PetscScalar     *funcVal, *interpolant;  in DMPlexComputeL2DiffVec()  local
1916 …PetscCall(PetscMalloc6(numComponents, &funcVal, numComponents, &interpolant, coordDim * Nq, &coord…  in DMPlexComputeL2DiffVec()
1954             PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMPlexComputeL2DiffVec()
1958 …ASSID) PetscCall(PetscFEInterpolate_Static((PetscFE)obj, &x[fieldOffset], &qgeom, q, interpolant));  in DMPlexComputeL2DiffVec()
1959 …TSCFV_CLASSID) PetscCall(PetscFVInterpolate_Static((PetscFV)obj, &x[fieldOffset], q, interpolant));  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()
2010   PetscScalar     *interpolant, *minterpolant, *earray;  in DMPlexComputeL2FluxDiffVecLocal()  local
2037 …PetscCall(PetscMalloc6(Nc * cdim, &interpolant, mNc * cdim, &minterpolant, cdim * (Nq + 1), &coord…  in DMPlexComputeL2FluxDiffVecLocal()
2059       PetscCall(PetscFEInterpolateGradient_Static(fe, 1, &x[0], &qgeom, q, interpolant));  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()
2151   PetscScalar     *interpolant, *valsum;  in DMPlexComputeClementInterpolant()  local
2188 …PetscCall(PetscMalloc6(Nc * 2, &valsum, Nc, &interpolant, cdim * Nq, &coords, Nq, &fegeom.detJ, cd…  in DMPlexComputeClementInterpolant()
2233           PetscCall(PetscFEInterpolate_Static((PetscFE)obj, &x[foff], &qgeom, q, interpolant));  in DMPlexComputeClementInterpolant()
2234           for (fc = 0; fc < fNc; ++fc) val[foff + fc] += interpolant[fc] * wt;  in DMPlexComputeClementInterpolant()
2255   PetscCall(PetscFree6(valsum, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMPlexComputeClementInterpolant()
2284   PetscScalar     *interpolant, *gradsum;  in DMPlexComputeGradientClementInterpolant()  local
2321 …oc6(coordDim * numComponents * 2, &gradsum, coordDim * numComponents, &interpolant, coordDim * Nq,…  in DMPlexComputeGradientClementInterpolant()
2364 …cCall(PetscFEInterpolateGradient_Static((PetscFE)obj, 1, &x[fieldOffset], &qgeom, q, interpolant));  in DMPlexComputeGradientClementInterpolant()
2368 …for (d = 0; d < coordDim; ++d) grad[fc * coordDim + d] += interpolant[fc * dim + d] * wt * fegeom.…  in DMPlexComputeGradientClementInterpolant()
2397   PetscCall(PetscFree6(gradsum, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMPlexComputeGradientClementInterpolant()