diff --git a/ComplexScaling-FV-bound.jl b/ComplexScaling-FV-bound.jl index f9b6891..2208c61 100644 --- a/ComplexScaling-FV-bound.jl +++ b/ComplexScaling-FV-bound.jl @@ -7,15 +7,15 @@ V_gauss(r2) = d = 3 n = 2 -N = 32 -L = 16 +N = 30 +L = 6 open("ComplexScaling-FV-bound.dat", "w") do f for ϕ = range(0.0, 0.5, length=11) println("Calculating ϕ=", ϕ) s = system{T}(d, n, N, L) H = Hamiltonian{T}(s, V_gauss, ϕ, 0, mode) - @time evals, _, info = eig(H, 5, resonances = false) + @time evals, _, info = eig(H, 10, resonances = false) dataline = vcat([ϕ], hcat(real.(evals), imag.(evals))'[:]) println(f, join(dataline, '\t'))