From 490e94262bd2797b1db457614973423a8259304b Mon Sep 17 00:00:00 2001
From: ysyapa <ysyapa@ncsu.edu>
Date: Wed, 9 Aug 2023 18:22:05 +0000
Subject: [PATCH] Increase N

---
 ComplexScaling-FV.dat | 0
 ComplexScaling-FV.jl  | 2 +-
 Hamiltonian.jl        | 2 +-
 3 files changed, 2 insertions(+), 2 deletions(-)
 create mode 100644 ComplexScaling-FV.dat

diff --git a/ComplexScaling-FV.dat b/ComplexScaling-FV.dat
new file mode 100644
index 0000000..e69de29
diff --git a/ComplexScaling-FV.jl b/ComplexScaling-FV.jl
index 6d7d0d2..4017bdb 100644
--- a/ComplexScaling-FV.jl
+++ b/ComplexScaling-FV.jl
@@ -7,7 +7,7 @@ V_gauss(r2) =
 
 d = 3
 n = 2
-N = 16
+N = 64
 ϕ = pi/6
 
 open("ComplexScaling-FV.dat", "w") do f
diff --git a/Hamiltonian.jl b/Hamiltonian.jl
index 321c343..f368cf6 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)
+    evals, evecs, info = eigsolve(H, x₀, levels, resonances ? :LI : :SR; ishermitian = H.hermitian, tol = tolerance, krylovdim = levels * 2)
     resonances || info.converged < levels && throw(error("Not enough convergence")) # don't check convergence for resonances
     if H.hermitian evals = real.(evals) end
     if H.mode == gpu_cutensor # to avoid possible GPU memory leak