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

n_imag=1

This commit is contained in:
ysyapa 2023-09-05 16:18:19 +00:00
parent 547be9fa98
commit 922807eca1
2 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
ComplexScaling-FV-bound.jl
View File

@ -9,12 +9,13 @@ d = 3
n = 2 n = 2
N = 30 N = 30
L = 6 L = 6
n_imag = 1
open("ComplexScaling-FV-bound.dat", "w") do f open("ComplexScaling-FV-bound.dat", "w") do f
for ϕ = range(0.0, 0.5, length=11) for ϕ = range(0.0, 0.5, length=11)
println("Calculating ϕ=", ϕ) println("Calculating ϕ=", ϕ)
s = system{T}(d, n, N, L) s = system{T}(d, n, N, L)
H = Hamiltonian{T}(s, V_gauss, ϕ, 0, mode) H = Hamiltonian{T}(s, V_gauss, ϕ, n_imag, mode)
@time evals, _, info = eig(H, 10, resonances = false) @time evals, _, info = eig(H, 10, resonances = false)
dataline = vcat([ϕ], hcat(real.(evals), imag.(evals))'[:]) dataline = vcat([ϕ], hcat(real.(evals), imag.(evals))'[:])

3
ComplexScaling-FV-res.jl
View File

@ -9,12 +9,13 @@ d = 3
n = 2 n = 2
N = 96 N = 96
ϕ = pi/6 ϕ = pi/6
n_imag = 1
open("ComplexScaling-FV-res.dat", "w") do f open("ComplexScaling-FV-res.dat", "w") do f
for L = range(20, 35, length=16) for L = range(20, 35, length=16)
println("Calculating L=", L) println("Calculating L=", L)
s = system{T}(d, n, N, L) s = system{T}(d, n, N, L)
H = Hamiltonian{T}(s, V_gauss, ϕ, 0, mode) H = Hamiltonian{T}(s, V_gauss, ϕ, n_imag, mode)
@time evals, _, info = eig(H, 40) @time evals, _, info = eig(H, 40)
dataline = vcat([L], hcat(real.(evals), imag.(evals))'[:]) dataline = vcat([L], hcat(real.(evals), imag.(evals))'[:])