I’ve just finished PolyhedralCubature.jl, a package allowing to evaluate the multiple integral of a function over a convex polytope, given by linear inequalities such as:

It remains to improve the doc. I also implemented this method in R and Python.

It would be nice for users who cannot afford a heavy dependency. For example, we have similar integration features in our GeoStats.jl stack and would benefit from a well-tested package for integration over polytopes.

In Python and R, I get the V-representation with the help of the C library cddlib. I think it is also used by Polyhedra, via CDDLib.jl. I’m not aware of another available implementation.

Got it. I wish I had the time to read the original papers describing the methods. I am assuming that it is an optimization problem that can be formulated and solved with native Julia packages.

Perhaps. And I don’t need the sense V -> H. Moreover I need H -> V in the polytope case only, without rays and lines. So maybe that’s not so difiicult.