Update for 0.5

Author: dlfivefifty

src/HaarMeasure.jl

@@ -65,18 +65,22 @@ end
 #Zyczkowski and Kus, Random unitary matrices, J. Phys. A: Math. Gen. 27,
 #4235–4245 (1994).
 
+
+if VERSION < v"0.5.0-dev"
+    macro blasfunc(x)
+       return :( $(BLAS.blasfunc(x) ))
+    end
+else
+    import Base.BLAS.@blasfunc
+end
+
+
 using Base.LinAlg: BlasInt
 for (s, elty) in (("dlarfg_", Float64),
                   ("zlarfg_", Complex128))
-
-    if Base.VERSION < v"0.5-"
-	blass = Base.blasfunc(s)
-    else
-	blass = eval(Base, :(@blasfunc($s)))
-    end
     @eval begin
         function larfg!(n::Int, α::Ptr{$elty}, x::Ptr{$elty}, incx::Int, τ::Ptr{$elty})
-	    ccall(($blass, Base.liblapack_name), Void,
+	    ccall((@blasfunc($s), LAPACK.liblapack), Void,
 		  (Ptr{BlasInt}, Ptr{$elty}, Ptr{$elty}, Ptr{BlasInt}, Ptr{$elty}),
                   &n, α, x, &incx, τ)
         end
@@ -124,4 +128,3 @@ function randfast(W::Haar, n::Int)
         error("beta = $beta not implemented")
     end
 end
-

src/InvariantEnsembles.jl

@@ -50,7 +50,7 @@ function legendrepolynomials(n)
   orthonormalpolynomials(1./sqrt(2),k->0,k->sqrt(k^2./(4k^2-1)),[-1.,1.],n)
 end
 function scaledhermitepolynomials(n)
-  orthonormalpolynomials(n^(1/4)/(1.π)^(1/4),k->0,k->sqrt((k)/(2n)),[-3.,3.],n)
+  orthonormalpolynomials(n^(1/4)/(1.0π)^(1/4),k->0,k->sqrt((k)/(2n)),[-3.,3.],n)
 end
 
 
@@ -179,11 +179,11 @@ function InvariantEnsemble(str::AbstractString,n::Integer)
         error("Legendre not implemented")
         #InvariantEnsembles(map(f->pad(f,2n),legendrepolynomials(n)))
     elseif(str == "CoshUnscaled")
-        InvariantEnsembleUnscaled("Cosh",x->cosh(x),[-2acosh(1.n),2acosh(1.n)],n)
+        InvariantEnsembleUnscaled("Cosh",x->cosh(x),[-2acosh(1.0n),2acosh(1.0n)],n)
     elseif(str == "EightUnscaled")
-        InvariantEnsembleUnscaled("Eight",x->x.^8,[-3.n^(1/8),3.n^(1/8)],n)
+        InvariantEnsembleUnscaled("Eight",x->x.^8,[-3.0n^(1/8),3.0n^(1/8)],n)
     elseif(str == "QuarticUnscaled")
-        InvariantEnsembleUnscaled("Quartic",x->x.^4,[-3.n^(1/4),3.n^(1/4)],n)
+        InvariantEnsembleUnscaled("Quartic",x->x.^4,[-3.0n^(1/4),3.0n^(1/4)],n)
     end
 end
 

src/densities/Semicircle.jl

@@ -9,6 +9,7 @@ with mean 0 and variance 1 (radius 2)
 immutable Semicircle{T<:Real} <: ContinuousUnivariateDistribution
     mean::T
     radius::T
+    Semicircle(μ::T,r::T) = new(μ,r)
 end
 Semicircle{T<:Real}(μ::T=0.0, r::T=2.0) = r > 0 ? Semicircle{T}(μ, r) :
     throw(ArgumentError("radius r must be positive, got $r"))

src/densities/TracyWidom.jl

@@ -9,7 +9,7 @@ Hermitian matrix.
 The cdf of Tracy-Widom is given by
 
 ``
-F_2 (s) = lim_{n→∞} Pr(√2 n^{1/6} (λₙ - √(2n) ≤ s) 
+F_2 (s) = lim_{n→∞} Pr(√2 n^{1/6} (λₙ - √(2n) ≤ s)
 ``
 
 References:
@@ -76,13 +76,13 @@ end
 """
 Samples the largest eigenvalue of the n × n GUE matrix
 """
-function rand(d::TracyWidom, n::Integer)
+function rand(d::TracyWidom, n::Int)
     n > 1 || error("Cannot construct $n × $n matrix")
     if n < 100
         k = n
     else #Exploit the fact that the eigenvector is concentrated in the top left corner
 	k = round(Int, n-10*n^(1/3)-1)
-    end	
+    end
     a=randn(n-k+1)
     b=[χ(i) for i=(n-1):-1:k]
     v=eigmax(SymTridiagonal(a, b))