  We discuss an algorithm for Tchakaloff-like compression of Quasi-MonteCarlo
(QMC) volume/surface integration on union of balls (multibubbles). The key
tools are Davis-Wilhelmsen theorem on the so-called Tchakaloff sets for
positive linear functionals on polynomial spaces, and Lawson-Hanson algorithm
for NNLS. We provide the corresponding Matlab package together with several
examples.
