Integral on `SimpleMesh` not yet implemented

[mikeingold]

Ideas sought. A simple version of this might be possible by simply decomposing the mesh into a set of Triangles and then returning the sum of their integrals, i.e.:

function _integral_2d(f, mesh::SimpleMesh, rule::IntegrationRule)
    triangle = _some_function_mesh_to_triangles(mesh)
    integrand(triangle) = _integral(f, triangle, rule)
    return sum(integrand, triangles)
end

[mikeingold]

I think this might just be as simple as something like:

function integral(f, mesh::SimpleMesh, ...)
    facets = mesh |> simplexify |> elements
    integrand(facet) = integral(f, facet, ...)
    return sum(integrand, facets)
end

which would seem to also be compatible with other meshes like subtypes of Grid.