1465   PetscScalar     *funcVal, *interpolant;  in DMPlexComputeL2DiffLocal()  local
1504 …PetscCall(PetscMalloc6(numComponents, &funcVal, numComponents, &interpolant, coordDim * (Nq + 1), …  in DMPlexComputeL2DiffLocal()
1550         PetscCall(PetscArrayzero(funcVal, Nc));  in DMPlexComputeL2DiffLocal()
1551         ierr = (*funcs[field])(coordDim, time, gcoords, Nc, funcVal, ctx);  in DMPlexComputeL2DiffLocal()
1555           PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMPlexComputeL2DiffLocal()
1557 …TransformApply_Internal(dm, &coords[coordDim * q], PETSC_FALSE, Nc, funcVal, funcVal, dm->transfor…  in DMPlexComputeL2DiffLocal()
1565 …Part(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]), (double)PetscRealPart(interpolant[fc]…  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 …PetscCall(PetscMalloc6(numComponents, &funcVal, coordDim * (Nq + 1), &coords, coordDim * coordDim …  in DMComputeL2GradientDiff_Plex()
1659         PetscCall(PetscArrayzero(funcVal, Nc));  in DMComputeL2GradientDiff_Plex()
1660         ierr = (*funcs[field])(coordDim, time, gcoords, n, Nc, funcVal, ctx);  in DMComputeL2GradientDiff_Plex()
1664           PetscCall(PetscFree6(funcVal, coords, fegeom.J, fegeom.invJ, interpolant, fegeom.detJ));  in DMComputeL2GradientDiff_Plex()
1666 …TransformApply_Internal(dm, &coords[coordDim * q], PETSC_FALSE, Nc, funcVal, funcVal, dm->transfor…  in DMComputeL2GradientDiff_Plex()
1679 …f %g\n", c, field, fc, (double)(PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fege…  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()
1813           for (fc = 0; fc < Nc; ++fc) funcVal[fc] = 0.;  in DMComputeL2FieldDiff_Plex()
1814           ierr = (*funcs[fields[f]])(dE, time, gcoords, Nc, funcVal, ctx);  in DMComputeL2FieldDiff_Plex()
1818             PetscCall(PetscFree6(funcVal, interpolant, coords, fegeom.detJ, fegeom.J, fegeom.invJ));  in DMComputeL2FieldDiff_Plex()
1820 …xBasisTransformApply_Internal(dm, &coords[dE * q], PETSC_FALSE, Nc, funcVal, funcVal, dm->transfor…  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()
1950           PetscCall((*funcs[field])(coordDim, time, &coords[q * coordDim], Nc, funcVal, ctx));  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()