Hi, I'm trying to evaluate integrals (symbolically or numerically) over polyhedral regions, using the new Polyhedron tools in Mathematica 12. From the documentation:
Polyhedron can be used a geometric region and a graphics primitive.
However, if I actually try to do this, it doesn't seem to work:
region = RandomPolyhedron[{"ConvexHull", 14}]
NIntegrate[x + y^2, {x, y, z} \[Element] region]
>NIntegrate::ilim: Invalid integration variable or limit(s) in {x,y,z}\[Element]Polyhedron[...]
It works fine with meshes:
R = ConvexHullMesh[RandomReal[1, {50, 3}]];
NIntegrate[x + y^2, {x, y, z} \[Element] R]
> 0.494771
Is this expected behavior or a bug? Alternatively, is there a way to query the vertices/faces of a Polyhedron so I might pass that information to BoundaryMeshRegion to convert it to a Mesh?
