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.,
[
423
]
and references therein).