Fix for MultivariateGaussianDist for 1d cases

If using the MultivariateGaussianDist to set a one dimensional prior, e.g., using it to set a multi-modal distribution on a single parameter, the sample method currently fails. This MR fixes the failure by making sure that the samples that are created are 1d arrays.