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