make space for eigenvectors complex valued

Jutho · Oct 23, 2018 · ec79ed2 · ec79ed2
1 parent 0c55176
commit ec79ed2
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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)]