From d3a25b2ccf442602bf237e8019cd5035c284d9d7 Mon Sep 17 00:00:00 2001
From: ysyapa <ysyapa@ncsu.edu>
Date: Sat, 21 Oct 2023 14:57:45 -0400
Subject: [PATCH] New ComplexScaling-FV calculation

---
 ...ng-FV-res.jl => ComplexScaling-FV-P-res.jl |  2 +-
 ...und.jl => ComplexScaling-FV-S-bound-phi.jl |  2 +-
 ComplexScaling-FV-S-res-phi.jl                | 24 +++++++++++++++++++
 3 files changed, 26 insertions(+), 2 deletions(-)
 rename ComplexScaling-FV-res.jl => ComplexScaling-FV-P-res.jl (91%)
 rename ComplexScaling-FV-bound.jl => ComplexScaling-FV-S-bound-phi.jl (90%)
 create mode 100644 ComplexScaling-FV-S-res-phi.jl

diff --git a/ComplexScaling-FV-res.jl b/ComplexScaling-FV-P-res.jl
similarity index 91%
rename from ComplexScaling-FV-res.jl
rename to ComplexScaling-FV-P-res.jl
index 2d7f869..f8ddfb8 100644
--- a/ComplexScaling-FV-res.jl
+++ b/ComplexScaling-FV-P-res.jl
@@ -11,7 +11,7 @@ N = 96
 ϕ = pi/6
 n_imag = 1
 
-open("ComplexScaling-FV-res.dat", "w") do f
+open("ComplexScaling-FV-P-res.dat", "w") do f
     for L = range(20, 35, length=16)
         println("Calculating L=", L)
         s = system{T}(d, n, N, L)
diff --git a/ComplexScaling-FV-bound.jl b/ComplexScaling-FV-S-bound-phi.jl
similarity index 90%
rename from ComplexScaling-FV-bound.jl
rename to ComplexScaling-FV-S-bound-phi.jl
index f852c44..fdba365 100644
--- a/ComplexScaling-FV-bound.jl
+++ b/ComplexScaling-FV-S-bound-phi.jl
@@ -11,7 +11,7 @@ N = 30
 L = 6
 n_imag = 1
 
-open("ComplexScaling-FV-bound.dat", "w") do f
+open("ComplexScaling-FV-S-bound-phi.dat", "w") do f
     for ϕ = range(0.0, 0.5, length=11)
         println("Calculating ϕ=", ϕ)
         s = system{T}(d, n, N, L)
diff --git a/ComplexScaling-FV-S-res-phi.jl b/ComplexScaling-FV-S-res-phi.jl
new file mode 100644
index 0000000..afe1f32
--- /dev/null
+++ b/ComplexScaling-FV-S-res-phi.jl
@@ -0,0 +1,24 @@
+include("Hamiltonian.jl")
+mode = cpu_tensor
+T = Float32 # single-precision mode
+
+V_gauss(r2) =
+    2 * exp(- ((sqrt(r2)-3)/1.5) ^ 2)
+
+d = 3
+n = 2
+N = 96
+L = 30
+n_imag = 1
+
+open("ComplexScaling-FV-S-res-phi.dat", "w") do f
+    for ϕ = range(0.1, 0.6, length=26)
+        println("Calculating ϕ=", ϕ)
+        s = system{T}(d, n, N, L)
+        H = Hamiltonian{T}(s, V_gauss, ϕ, n_imag, mode)
+        @time evals, _, info = eig(H, 40, resonances = true)
+
+        dataline = vcat([ϕ], hcat(real.(evals), imag.(evals))'[:])
+        println(f, join(dataline, '\t'))
+    end
+end