Is there an efficient way to generate multinomial random variables in parallel?

Question

numpy.random has the following function to generate multinomial random samples.

multinomial(n, p, size)

But I wonder if there is an e