Tropical sampling from Feynman measures

We introduce an algorithm that samples a set of loop momenta distributed as a given Feynman integrand. The algorithm uses the tropical sampling method and can be applied to evaluate phase-space-type integrals efficiently. We provide an implementation, momtrop, and apply it to a series of relevant integrals from the loop-tree duality framework. Compared to naive sampling methods, we observe convergence speedups by factors of more than ${10}^6$.