Internal function to calculate the density of multivariate normal distribution.
dmvnrm_arma(x, mean, sigma, logd = FALSE)
Useful links