The function ` rand(d::Wigner{2}, n::Int)` returns incorrect matrix???

I think that the following code is missing a factor of 2 dividing when the function returns.
```
function rand(d::Wigner{2}, n::Int)
A = randn(n, n) + im*randn(n, n)
normalization = √(4*n)
return Hermitian((A + A') / normalization)
end
```
That is, replace this line `return Hermitian((A + A') / normalization)` with this other one `return Hermitian( (A + A') / 2 / normalization )`.

Thanks and regards

Pablo

Locate the `rand(d::Wigner{2}, n::Int)` entry point and inspect how the generated Hermitian matrix is normalized. Verify the expected scaling for the Wigner matrix and compare it with the reported factor-of-two change. Done means the function returns correctly normalized matrices and relevant tests pass.

Written by the indexing model from the issue text.