Fixed Haar sampling

Author: dlfivefifty

src/HaarMeasure.jl

@@ -22,13 +22,13 @@ import Base.rand
 #TODO implement O(n^2) method
 function rand(W::Haar, n::Integer, doCorrection::Integer)
     beta = W.beta
-    M=rand(Ginibre(n, beta))
+    M=rand(Ginibre(beta,n))
     q,r=qr(M)
     if doCorrection==0
         q
     elseif doCorrection==1 
         if beta==1
-            L = sign(rand(n)-0.5)
+            L = sign(rand(n).-0.5)
         elseif beta==2
             L = exp(im*rand(n)*2pi)
         elseif beta==4
@@ -53,7 +53,7 @@ end
 #By default, always do piecewise correction
 #For most applications where you use the HaarMatrix as a similarity transform
 #it doesn't matter, but better safe than sorry... let the user choose else
-rand(W::Haar) = rand(W, 1)
+rand(W::Haar,n::Integer) = rand(W,n, 1)
 
 #A utility method to check if the piecewise correction is needed
 #This checks the R part of the QR factorization; if correctly done,