From e52d5aae837345016d77f7e7a8cf3e58b6a2a0c9 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?=C5=81ukasz=20Pawela?= Date: Thu, 22 Aug 2019 18:52:27 +0200 Subject: [PATCH] fix wigner ensemble --- randommatrices/src/wigner.jl | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/randommatrices/src/wigner.jl b/randommatrices/src/wigner.jl index 0c7186c..5261619 100644 --- a/randommatrices/src/wigner.jl +++ b/randommatrices/src/wigner.jl @@ -2,16 +2,18 @@ export WignerEnsemble struct WignerEnsemble{β} <: QIContinuousMatrixDistribution d::Int + g::GinibreEnsemble{β} function WignerEnsemble{β}(d::Int) where β β == 4 && mod(d, 2) == 1 ? throw(ArgumentError("Dim must even")) : () - new(d) + g = GinibreEnsemble{β}(d) + new(d, g) end end WignerEnsemble(d::Int) = WignerEnsemble{2}(d) function rand(rng::AbstractRNG, w::WignerEnsemble{β}) where β - z = rand(rng, GinibreEnsemble{β}(w.d)) + z = rand(rng, w.g) (z + z') / 2sqrt(2β * w.d) end \ No newline at end of file