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

Small changes

This commit is contained in:
Nuwan Yapa 2024-03-08 18:35:19 -05:00
parent 8bcb5c8c12
commit 4175152a8d
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
ho_basis_3body.jl
View File

@ -2,7 +2,7 @@ using Arpack, SparseArrays
include("ho_basis.jl") include("ho_basis.jl")
include("p_space.jl") include("p_space.jl")
E_max = 10 E_max = 20
ω = 1.0 ω = 1.0
Λ = 0 Λ = 0
m = 1.0 m = 1.0
@ -35,6 +35,6 @@ V_l(l, n1, n2) = Va * V_Gaussian(Ra, l, n1, n2; ω=ω) + Vb * V_Gaussian(Rb, l,
println("Calculating spectrum") println("Calculating spectrum")
@time "H" H = T1 + T2 + V1 + V2 @time "H" H = T1 + T2 + V1 + V2
@time "Eigenvalues" evals, _ = eigs(H, nev=5, which=:SM, maxiter=10000, tol=1e-6, ritzvec=false) @time "Eigenvalues" evals, _ = eigs(H, nev=3, ncv=30, which=:SR, maxiter=5000, tol=1e-5, ritzvec=false, check=1)
display(evals) display(evals)