From 305e3759469f46e3e5faa08bc26f0b594908c63e Mon Sep 17 00:00:00 2001
From: Daniel Karrasch <daniel.karrasch@posteo.de>
Date: Fri, 29 Jul 2022 17:31:57 +0200
Subject: [PATCH] Avoid broadcasting over a `Q` matrix

---
 src/RandomQuantum.jl | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/src/RandomQuantum.jl b/src/RandomQuantum.jl
index 579e33c..d855021 100644
--- 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