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

Seperate .jl files instead of notebook

This commit is contained in:
ysyapa 2023-05-23 19:11:36 +00:00
parent d3726c725b
commit 700d2b4308
2 changed files with 37 additions and 65 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

37
3b_res_from_paper.jl Normal file
View File

@ -0,0 +1,37 @@
# 10.1007/s00601-020-01550-8
# Fig. 7
# E_R = 4.18(8)
#./En.run -d 3 -n 3 -N 16 -c pot=v_shifted_gauss,v0=2.0,r=1.5,a=3.0 -c n_eig=20 -c which=li -c tol=1e-6 -L 16 -c phi=0.3 -v
include("Hamiltonian.jl")
mode = gpu_cutensor
T = Float32 # single-precision mode
using Plots
V_gauss(r2) =
2 * exp(-((sqrt(r2) - 3) / 1.5) ^ 2)
d = 3
n = 3
N = 14
L = 16
μ = 0.5
n_imag = 0
for ϕ::T in 0.2:0.05:0.4
s = system{T}(d, n, N, L)
H = Hamiltonian{T}(s, V_gauss, ϕ, n_imag, mode)
@time evals, _, info = eig(H, 20)
print(info.numops, " operations")
display(evals)
scatter(real.(evals), imag.(evals); legend=false)
xlabel!("Re E")
ylabel!("Im E")
xlims!(0, 6)
ylims!(-0.6, 0)
savefig("temp/phi$(Int(round(ϕ * 100))).png")
end

65
calculations.ipynb
View File

@ -1,65 +0,0 @@
{
"cells": [
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"# 10.1007/s00601-020-01550-8\n",
"# Fig. 7\n",
"# E_R = 4.18(8)\n",
"\n",
"#./En.run -d 3 -n 3 -N 16 -c pot=v_shifted_gauss,v0=2.0,r=1.5,a=3.0 -c n_eig=20 -c which=li -c tol=1e-6 -L 16 -c phi=0.3 -v\n",
"\n",
"include(\"Hamiltonian.jl\")\n",
"mode = gpu_cutensor\n",
"T = Float32 # single-precision mode\n",
"\n",
"using Plots\n",
"\n",
"V_gauss(r2) =\n",
" 2 * exp(-((sqrt(r2) - 3) / 1.5) ^ 2)\n",
"\n",
"d = 3\n",
"n = 3\n",
"N = 14\n",
"L = 16\n",
"μ = 0.5\n",
"n_imag = 0\n",
"\n",
"for ϕ::T in 0.2:0.05:0.4\n",
" s = system{T}(d, n, N, L)\n",
" H = Hamiltonian{T}(s, V_gauss, ϕ, n_imag, mode)\n",
" @time evals, _, info = eig(H, 20)\n",
" \n",
" print(info.numops, \" operations\")\n",
" display(evals)\n",
"\n",
" scatter(real.(evals), imag.(evals); legend=false)\n",
" xlabel!(\"Re E\")\n",
" ylabel!(\"Im E\")\n",
" xlims!(0, 6)\n",
" ylims!(-0.6, 0)\n",
" savefig(\"temp/phi$(Int(round(ϕ * 100))).png\")\n",
"end"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Julia 1.8.5",
"language": "julia",
"name": "julia-1.8"
},
"language_info": {
"file_extension": ".jl",
"mimetype": "application/julia",
"name": "julia",
"version": "1.8.5"
},
"orig_nbformat": 4
},
"nbformat": 4,
"nbformat_minor": 2
}