testeigvaljpdf(d::GaussianLaguerre{β}, lambda})

Author: apkille

src/GaussianEnsembles.jl

@@ -289,7 +289,7 @@ function eigvaljpdf(d::GaussianLaguerre{β}, lambda::Vector{Eigenvalue}) where {
      c *= gamma(1 + β/2)/(gamma(1 + β*j/2)*gamma(z))
   end
 
-  Prod = prod(lambda.^(a-p)) #Calculate Laguerre product term
+  Prod = prod(lambda.^(d.a-p)) #Calculate Laguerre product term
   Energy = sum(lambda)/2 #Calculate argument of exponential
   return c * VandermondeDeterminant(lambda, β) * Prod * exp(-Energy)
 end

test/GaussianEnsembles.jl

@@ -47,8 +47,8 @@ for (β, T, N) in [(1, Real, n), (2, Complex, n), (4, Complex, 2n)]
         #vd = RandomMatrices.VandermondeDeterminant(vals, β)
         #@test isa(vd, Real)
 
-        #ed = eigvaljpdf(d, vals)
-        #@test isa(ed, Real)
+        ed = eigvaljpdf(d, vals)
+        @test isa(ed, Real)
     end
 end