yapanuwan
/
DVR-jl
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Identified two discrepancies

This commit is contained in:
ysyapa 2023-03-19 14:46:00 -04:00
parent f7fc232d8b
commit 0f53e4e719
1 changed files with 45 additions and 31 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

76
tester.ipynb
View File

@ -6,7 +6,45 @@
"metadata": {},
"outputs": [],
"source": [
"# ./En.run -d 3 -n 3 -e 5 -c eps=0 -c pot=v_gauss,v0=-4,r=2 -N 6 -L 5:14 -c n_imag=1\n",
"# ./En.run -d 3 -n 3 -e 5 -c pot=v_gauss,v0=-4,r=2 -N 8 -L 4 -c n_imag=0\n",
"# Calculating...\n",
"# -> N = 8\n",
"# -> L = 4\n",
"# -> Spectrum = {-6.07632,-3.81486,-3.71969,-3.71968,-3.38263}...\n",
"# Done.\n",
"# -> Time used = 12s\n",
"\n",
"include(\"CPU.jl\")\n",
"\n",
"T=Float32\n",
"\n",
"function V_test(r2::T)::T\n",
" return -4*exp(-r2/4)\n",
"end\n",
"\n",
"N=8\n",
"L=4.0f0\n",
"n_imag=0\n",
"\n",
"H=HOperator{T}(V_test,3,3,N,L,0.5f0,n_imag)\n",
"@time evals,evecs,info=eig(H,5)\n",
"print(info.numops,\" operations : \")\n",
"println(evals)\n"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"# ./En.run -d 3 -n 3 -e 5 -c pot=v_gauss,v0=-4,r=2 -N 6 -L 14 -c n_imag=0\n",
"# Calculating...\n",
"# -> N = 6\n",
"# -> L = 14\n",
"# -> Spectrum = {-8.49538,-3.35492,-3.34356,-3.32830,-3.07909}...\n",
"# Done.\n",
"# -> Time used = 0.229049s\n",
"\n",
"include(\"CPU.jl\")\n",
"\n",
@ -17,37 +55,13 @@
"end\n",
"\n",
"N=6\n",
"for L::T in 5.0:14.0\n",
" H=HOperator{T}(V_test,3,3,N,L,0.5f0,1)\n",
" @time evals,evecs,info=eig(H,5)\n",
" print(info.numops,\" operations : \")\n",
" println(evals)\n",
"end"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"# ./En.run -d 3 -n 2 -e 5 -c eps=0 -c pot=v_gauss,v0=-4,r=2 -N 32 -L 5:14 -c n_imag=0\n",
"L=14.0f0\n",
"n_imag=0\n",
"\n",
"include(\"CPU.jl\")\n",
"\n",
"T=Float64\n",
"\n",
"function V_test(r2::T)::T\n",
" return -4*exp(-r2/4)\n",
"end\n",
"\n",
"N=32\n",
"for L::T in 5.0:14.0\n",
" H=HOperator{T}(V_test,3,2,N,L)\n",
" @time evals,evecs,info=eig(H,5)\n",
" print(info.numops,\" operations : \")\n",
" println(evals)\n",
"end"
"H=HOperator{T}(V_test,3,3,N,L,0.5f0,n_imag)\n",
"@time evals,evecs,info=eig(H,5)\n",
"print(info.numops,\" operations : \")\n",
"println(evals)"
]
}
],