22 lines
448 B
Julia
22 lines
448 B
Julia
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
|