include("Hamiltonian.jl")
mode = cpu_tensor
T = Float32 # single-precision mode

using Plots

V_gauss(r2) =
    -10 * exp(-(sqrt(r2)) ^ 2)

d = 3
n = 2
N = 32
L = 16
ϕ = pi/5

s = system{T}(d, n, N, L)
H = Hamiltonian{T}(s, V_gauss, ϕ, 0, mode)
@time evals, _, info = eig(H, 20)

print(info.numops, " operations")
display(evals)

scatter(real.(evals), imag.(evals); legend=false)
xlabel!("Re E")
ylabel!("Im E")
xlims!(0, 0.5)
ylims!(-0.5, 0)
savefig("temp/ComplexScaling-FV.png")