include("Hamiltonian.jl")

T=Float32

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

for (n,N) in [(2,16), (3,8)]
    println("\n$n-body system with N=$N")
    n_image=0
    for L::T in 5.0:9.0
        print("L=$L: ")
        s=system{T}(3,n,N,L)
        print("Constructing H...")
        H=Hamiltonian{T}(s,V_test,0.0,n_image,cpu_tensor)
        print("Diagonalizing...")
        evals,_,_ = eig(H,5)
        println(real.(evals))
    end
end