23 lines
454 B
Julia
23 lines
454 B
Julia
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 [16]
|
|
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 |