  We present a brief survey on the compression of discrete measures by
Caratheodory-Tchakaloff Subsampling, its implementation by Linear or Quadratic
Programming and the application to multivariate polynomial Least Squares. We
also give an algorithm that computes the corresponding Caratheodory-Tchakaloff
(CATCH) points and weights for polynomial spaces on compact sets and manifolds
in 2D and 3D.
