Multivariate Poisson Distribution Generated via Reduction from Independent Poisson Variates

Let's say that we are given a k number of random variables following Poisson distribution that are individually dependent and which forms multivariate Poisson distribution. We particularly dealt with a method of creating random numbers that satisfies the covariance matrix, where the elements of covariance matrix are parameters forming a multivariate Poisson distribution. To create such random numbers, we propose a new algorithm based on the method reducing the number of parameter set and deal with its relationship to the Park et al.(1996) algorithm used in creating multivariate Bernoulli random numbers.