From 3d29c323d0700b0c3e3a1fc4575c2dcd9e1d22ab Mon Sep 17 00:00:00 2001
From: ysyapa <ysyapa@ncsu.edu>
Date: Fri, 25 Aug 2023 21:03:25 +0000
Subject: [PATCH] New calculation for paper

---
 ComplexScaling-FV-bound.jl                    | 23 +++++++++++++++++++
 ...xScaling-FV.jl => ComplexScaling-FV-res.jl |  2 +-
 2 files changed, 24 insertions(+), 1 deletion(-)
 create mode 100644 ComplexScaling-FV-bound.jl
 rename ComplexScaling-FV.jl => ComplexScaling-FV-res.jl (91%)

diff --git a/ComplexScaling-FV-bound.jl b/ComplexScaling-FV-bound.jl
new file mode 100644
index 0000000..f9b6891
--- /dev/null
+++ b/ComplexScaling-FV-bound.jl
@@ -0,0 +1,23 @@
+include("Hamiltonian.jl")
+mode = cpu_tensor
+T = Float32 # single-precision mode
+
+V_gauss(r2) =
+    -10 * exp(-(sqrt(r2)) ^ 2)
+
+d = 3
+n = 2
+N = 32
+L = 16
+
+open("ComplexScaling-FV-bound.dat", "w") do f
+    for ϕ = range(0.0, 0.5, length=11)
+        println("Calculating ϕ=", ϕ)
+        s = system{T}(d, n, N, L)
+        H = Hamiltonian{T}(s, V_gauss, ϕ, 0, mode)
+        @time evals, _, info = eig(H, 5, resonances = false)
+
+        dataline = vcat([ϕ], hcat(real.(evals), imag.(evals))'[:])
+        println(f, join(dataline, '\t'))
+    end
+end
diff --git a/ComplexScaling-FV.jl b/ComplexScaling-FV-res.jl
similarity index 91%
rename from ComplexScaling-FV.jl
rename to ComplexScaling-FV-res.jl
index 61f4acb..e3b1cfe 100644
--- a/ComplexScaling-FV.jl
+++ b/ComplexScaling-FV-res.jl
@@ -10,7 +10,7 @@ n = 2
 N = 96
 ϕ = pi/6
 
-open("ComplexScaling-FV.dat", "w") do f
+open("ComplexScaling-FV-res.dat", "w") do f
     for L = range(20, 35, length=16)
         println("Calculating L=", L)
         s = system{T}(d, n, N, L)