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

Move EC calculations to a folder

This commit is contained in:
Nuwan Yapa 2024-06-27 18:20:19 -04:00
parent 449c93817a
commit d1a76c1909
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
3body_resonance_EC.jl → calculations/3body_R2R_EC.jl
View File

@ -1,5 +1,5 @@
using DelimitedFiles, Plots using DelimitedFiles, Plots
include("3body_resonance.jl") include("../3body_resonance.jl")
current_E = 5.9673 - 0.0006im current_E = 5.9673 - 0.0006im
@ -19,8 +19,8 @@ training_vecs = Vector{ComplexF64}[]
for c in training_c for c in training_c
println("Training for c = $c") println("Training for c = $c")
H = H0 + c .* V local H = H0 + c .* V
evals, evecs = eigs(H, nev=3, ncv=24, which=:LI, maxiter=5000, tol=1e-5, ritzvec=true, check=1) local evals, evecs = eigs(H, nev=3, ncv=24, which=:LI, maxiter=5000, tol=1e-5, ritzvec=true, check=1)
global current_E = nearest(evals, current_E) global current_E = nearest(evals, current_E)
push!(training, current_E) push!(training, current_E)
@ -34,8 +34,8 @@ V_EC = transpose(EC_basis) * V * EC_basis
for c in extrapolating_c for c in extrapolating_c
println("Extrapolating for c = $c") println("Extrapolating for c = $c")
H = H0 + c .* V local H = H0 + c .* V
evals, evecs = eigs(H, nev=3, ncv=24, which=:LI, maxiter=5000, tol=1e-5, ritzvec=true, check=1) local evals, evecs = eigs(H, nev=3, ncv=24, which=:LI, maxiter=5000, tol=1e-5, ritzvec=true, check=1)
global current_E = nearest(evals, current_E) global current_E = nearest(evals, current_E)
push!(exact, current_E) push!(exact, current_E)