25 lines
575 B
Julia
25 lines
575 B
Julia
include("Hamiltonian.jl")
|
|
mode = cpu_tensor
|
|
T = Float32 # single-precision mode
|
|
|
|
V_gauss(r2) =
|
|
2 * exp(- ((sqrt(r2)-3)/1.5) ^ 2)
|
|
|
|
d = 3
|
|
n = 2
|
|
N = 96
|
|
L = 30
|
|
n_imag = 1
|
|
|
|
open("ComplexScaling-FV-S-res-phi.dat", "w") do f
|
|
for ϕ = range(0.1, 0.6, length=26)
|
|
println("Calculating ϕ=", ϕ)
|
|
s = system{T}(d, n, N, L)
|
|
H = Hamiltonian{T}(s, V_gauss, ϕ, n_imag, mode)
|
|
@time evals, _, info = eig(H, 40, resonances = true)
|
|
|
|
dataline = vcat([ϕ], hcat(real.(evals), imag.(evals))'[:])
|
|
println(f, join(dataline, '\t'))
|
|
end
|
|
end
|