New calculation for paper

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

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)