From b21242ad4976975582a87acb3647c691bd012440 Mon Sep 17 00:00:00 2001 From: ysyapa Date: Fri, 25 Aug 2023 21:20:14 +0000 Subject: [PATCH] Double krylovdim again --- Hamiltonian.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Hamiltonian.jl b/Hamiltonian.jl index 3be10f3..a206677 100644 --- a/Hamiltonian.jl +++ b/Hamiltonian.jl @@ -134,7 +134,7 @@ 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) + evals, evecs, info = eigsolve(H, x₀, levels, resonances ? :LI : :SR; ishermitian = H.hermitian, tol = tolerance, krylovdim = levels * 8) 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