Multivariate Gamma returns a single value

[samuel-verburg]

Description

Only one value is returned when sampling from a multivariate Gamma distribution.

Example to reproduce

import cuqi
from cuqi.distribution import Gamma
d = Gamma(1, 1e-4, geometry=50)
ds = d.sample()
ds

Out[1]:
CUQIarray: NumPy array wrapped with geometry.

Geometry:
_DefaultGeometry1D(50,)

Parameters:
True

Array:
CUQIarray(680.24372138)