Testing script
This commit is contained in:
parent
ded9e3a162
commit
443aaae391
|
|
@ -1,51 +0,0 @@
|
||||||
{
|
|
||||||
"cells": [
|
|
||||||
{
|
|
||||||
"cell_type": "code",
|
|
||||||
"execution_count": 1,
|
|
||||||
"metadata": {},
|
|
||||||
"outputs": [],
|
|
||||||
"source": [
|
|
||||||
"include(\"Hamiltonian.jl\")\n",
|
|
||||||
"\n",
|
|
||||||
"println(\"Running with \",Threads.nthreads(),\" thread(s)\")\n",
|
|
||||||
"\n",
|
|
||||||
"T=Float32\n",
|
|
||||||
"\n",
|
|
||||||
"function V_test(r2)\n",
|
|
||||||
" return -4*exp(-r2/4)\n",
|
|
||||||
"end\n",
|
|
||||||
"\n",
|
|
||||||
"n = 2\n",
|
|
||||||
"N = 16\n",
|
|
||||||
"println(\"\\n$n-body system with N=$N\")\n",
|
|
||||||
"\n",
|
|
||||||
"for L::T in 5.0:9.0\n",
|
|
||||||
" println(\"L=$L\")\n",
|
|
||||||
" println(\"Constructing Hamiltonian\")\n",
|
|
||||||
" s=system{T}(3,n,N,L,0.5,all)\n",
|
|
||||||
" @time H=Hamiltonian{T}(s,V_test,0,0,cpu_tensor)\n",
|
|
||||||
" println(\"Solving eigenvalues\")\n",
|
|
||||||
" @time evals,_,_ = eig(H,5)\n",
|
|
||||||
" println(real.(evals))\n",
|
|
||||||
"end"
|
|
||||||
]
|
|
||||||
}
|
|
||||||
],
|
|
||||||
"metadata": {
|
|
||||||
"kernelspec": {
|
|
||||||
"display_name": "Julia 1.9.0",
|
|
||||||
"language": "julia",
|
|
||||||
"name": "julia-1.9"
|
|
||||||
},
|
|
||||||
"language_info": {
|
|
||||||
"file_extension": ".jl",
|
|
||||||
"mimetype": "application/julia",
|
|
||||||
"name": "julia",
|
|
||||||
"version": "1.9.0"
|
|
||||||
},
|
|
||||||
"orig_nbformat": 4
|
|
||||||
},
|
|
||||||
"nbformat": 4,
|
|
||||||
"nbformat_minor": 2
|
|
||||||
}
|
|
||||||
|
|
@ -0,0 +1,23 @@
|
||||||
|
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,all)
|
||||||
|
@time H=Hamiltonian{T}(s,V_test,0,0,cpu_tensor)
|
||||||
|
println("Solving eigenvalues")
|
||||||
|
@time evals,_,_ = eig(H,5)
|
||||||
|
println(real.(evals))
|
||||||
|
end
|
||||||
Loading…
Reference in New Issue