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

Modified Jacobi

This commit is contained in:
Nuwan Yapa 2024-04-10 14:51:26 -04:00
parent fa00c75949
commit 3b3f096913
1 changed files with 2 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
ho_basis_3body.jl
View File

@ -11,9 +11,8 @@ Va = -2
Ra = 2 Ra = 2
μ1 = m * 1/2 μ1 = m * 1/2
μ2 = m * 2/3 μ2 = m * 2/3 * 4
c = sqrt(2) c = sqrt(2)
c2 = 2
println("No of threads = ", Threads.nthreads()) println("No of threads = ", Threads.nthreads())
@ -27,7 +26,7 @@ println("Basis size = ", length(Es))
println("Constructing KE matrices") println("Constructing KE matrices")
@time "T1" T1 = sp_T_matrix(n1s, l1s; mask=mask1, ω=ω, μ=μ1) @time "T1" T1 = sp_T_matrix(n1s, l1s; mask=mask1, ω=ω, μ=μ1)
@time "T2" T2 = sp_T_matrix(n2s, l2s; mask=mask2, ω=ω, μ=c2^2 * μ2) @time "T2" T2 = sp_T_matrix(n2s, l2s; mask=mask2, ω=ω, μ=μ2)
println("Constructing PE matrices") println("Constructing PE matrices")
V1_elem(l, n1, n2) = Va * V_Gaussian(Ra, l, n1, n2; ω=ω) V1_elem(l, n1, n2) = Va * V_Gaussian(Ra, l, n1, n2; ω=ω)