From ec79ed21ee9fc9b48e368387638d471dc2ee32be Mon Sep 17 00:00:00 2001
From: Jutho <jutho.haegeman@ugent.be>
Date: Tue, 23 Oct 2018 14:35:44 +0200
Subject: [PATCH] make space for eigenvectors complex valued

---
 src/eigsolve/arnoldi.jl | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/eigsolve/arnoldi.jl b/src/eigsolve/arnoldi.jl
index f763c48..1f55b86 100644
--- a/src/eigsolve/arnoldi.jl
+++ b/src/eigsolve/arnoldi.jl
@@ -110,7 +110,7 @@ function eigsolve(A, x₀, howmany::Int, which::Selector, alg::Arnoldi)
         [B*v for v in cols(V)]
     end
     residuals = let r = residual(fact)
-        [mul!(similar(r), r, last(v)) for v in cols(V)]
+        [mul!(similar(r, complex(eltype(r))), r, last(v)) for v in cols(V)]
     end
     normres = [abs(last(v)) for v in cols(V)]