Unify the Cholesky mvnorm with the mvnorm based on correlations.

The unified family should have an interface such as:

bamlss_mvnorm(k = 2, type = c("Cholesky", "cor"), ...)