Avoid broadcasting over a Q matrix

BBN-Q · Jul 29, 2022 · 305e375 · 305e375
1 parent 39f7496
commit 305e375
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/RandomQuantum.jl b/src/RandomQuantum.jl
@@ -130,10 +130,10 @@ end
 
 function rand(dist::ClosedHaarEnsemble)
     X = rand(GinibreEnsemble(dist.dim))
-    Q,_ = qr(X)
+    Q = qr(X).Q * I # yields dense square matrix
     d = diag(Q)
-    d = d ./ abs.(d)
-    Q = Q ./ d
+    d ./= abs.(d)
+    Q ./= d
     return Q
 end