Similar issue for

distrib = MultivariateNormal(loc=torch.zeros(input_dimensionality), covariance_matrix=torch.diag(torch.ones(input_dimensionality[0])))
distrib.sample((N,)).cuda()

for N=5 and input_dimensionality=16K