Instead of Monte Carlo simulations one may use a sufficiently dense grid in the 3D space of binary parameters and integrate over this grid (see, e.g., [423Jump To The Next Citation Point] and references therein).