include("Hamiltonian.jl")

println("Running with ",Threads.nthreads()," thread(s)")

T=Float32

function V_test(r2)
    return -4*exp(-r2/4)
end

n = 2
N = 16
println("\n$n-body system with N=$N")

for L::T in 5.0:9.0
    println("L=$L")
    println("Constructing Hamiltonian")
    s=system{T}(3,n,N,L,0.5,A1)
    @time H=Hamiltonian{T}(s,V_test,0,0,cpu_tensor)
    println("Solving eigenvalues")
    @time evals,_,_ = eig(H,5)
    println(real.(evals))
end