From bc8b026d4d35128bc6822614971fc6dd03e471ad Mon Sep 17 00:00:00 2001 From: ysyapa Date: Tue, 5 Sep 2023 23:14:56 +0000 Subject: [PATCH] 2-body works --- Hamiltonian.jl | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/Hamiltonian.jl b/Hamiltonian.jl index 8d23a73..db5c174 100644 --- a/Hamiltonian.jl +++ b/Hamiltonian.jl @@ -141,7 +141,8 @@ function eig(H::Hamiltonian{T}, levels::Int; resonances = !H.hermitian)::Tuple{V x₀ = CUDA.rand(Complex{T}, vectorDims(H)...) synchronize() end - evals, evecs, info = eigsolve(H, x₀, levels, resonances ? :LI : :SR; ishermitian = H.hermitian, tol = tolerance, krylovdim = levels * 4) + KrylovKit_hermitian = H.hermitian && H.s.sym == all + evals, evecs, info = eigsolve(H, x₀, levels, resonances ? :LI : :SR; ishermitian = KrylovKit_hermitian, tol = tolerance, krylovdim = levels * 4) info.converged < levels && throw(error("Not enough convergence")) if H.hermitian evals = real.(evals) end if H.mode == gpu_cutensor # to avoid possible GPU memory leak