How to create multivariate normal samples in C++ given mean vector and covariance matrix?

There is no torch::distributions::multivariate_normal::MultivariateNormal in C++. So can anyone please tell how to achieve this in C++ frontend.

Thanks in advance.