MAR_mv_MKR function

C++ function to estimate Pitman-Yor multivariate mixtures via marginal sampler - PRODUCT KERNEL