From c38e4c8c0e1185954d98429baa8151e608cae3b3 Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <ysyapa@ncsu.edu>
Date: Mon, 14 Aug 2023 10:46:22 -0400
Subject: [PATCH 1/2] Ignore calculation outputs

---
 .gitignore | 6 ++++++
 1 file changed, 6 insertions(+)

diff --git a/.gitignore b/.gitignore
index e4f443c..1fbef14 100644
--- a/.gitignore
+++ b/.gitignore
@@ -1,3 +1,9 @@
+# Calculation outputs
+*.dat
+*.csv
+*.hdf5
+*.out
+
 # Temporary and scratch files
 temp/
 scratch/

From 7e419aee266fa9bebf0cd112ebf8e75352680938 Mon Sep 17 00:00:00 2001
From: ysyapa <ysyapa@ncsu.edu>
Date: Fri, 18 Aug 2023 21:36:09 +0000
Subject: [PATCH 2/2] Increase krylovdim and check convergence

---
 Hamiltonian.jl | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/Hamiltonian.jl b/Hamiltonian.jl
index 324346c..3be10f3 100644
--- a/Hamiltonian.jl
+++ b/Hamiltonian.jl
@@ -134,8 +134,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 * 2)
-    resonances || info.converged < levels && throw(error("Not enough convergence")) # don't check convergence for resonances
+    evals, evecs, info = eigsolve(H, x₀, levels, resonances ? :LI : :SR; ishermitian = H.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
         CUDA.reclaim()