Export .dat file

ComplexScaling-FV.jl

@@ -2,27 +2,22 @@ include("Hamiltonian.jl")
 mode = cpu_tensor
 T = Float32 # single-precision mode
 
-using Plots
-
 V_gauss(r2) =
     -10 * exp(-(sqrt(r2)) ^ 2)
 
 d = 3
 n = 2
-N = 32
-L = 16
-ϕ = pi/5
+N = 16
+ϕ = pi/6
 
-s = system{T}(d, n, N, L)
-H = Hamiltonian{T}(s, V_gauss, ϕ, 0, mode)
-@time evals, _, info = eig(H, 20)
+open("ComplexScaling-FV.dat", "w") do f
+    for L = range(16, 32, length=33)
+        println("Calculating L=", L)
+        s = system{T}(d, n, N, L)
+        H = Hamiltonian{T}(s, V_gauss, ϕ, 0, mode)
+        @time evals, _, info = eig(H, 40)
 
-print(info.numops, " operations")
-display(evals)
-
-scatter(real.(evals), imag.(evals); legend=false)
-xlabel!("Re E")
-ylabel!("Im E")
-xlims!(0, 0.5)
-ylims!(-0.5, 0)
-savefig("temp/ComplexScaling-FV.png")
+        dataline = vcat([L], hcat(real.(evals), imag.(evals))'[:])
+        println(f, join(dataline, '\t'))
+    end
+end