Update tests for multivariateGaussianRand

Now we automatically generate tests for the permitted input forms.

Also fixes a few problems with argument handling that the tests caught.

Added support for passing diagonal covariance matrix as a vector of just
the diagonal enries.

Limitations: tests do not cover the sigma = NxD case (i.e. N distinct
diagonal covariance vectors) right now.