How to implement differentiable inverse cdf of bivariate Gaussian mixture model?

I’m looking for a way to implement the inverse cumulative distribution function of bivariate Gaussian mixture model. Is there a straightforward way to do this in PyTorch?