New ComplexScaling-FV calculation

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)

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)

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