use the vectorized implementation of marginal distribution computation