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

Critical error in Moshinsky transform fixed

This commit is contained in:
Nuwan Yapa 2024-03-12 19:41:16 -04:00
parent 3c72bbc7b6
commit 073972d0eb
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
ho_basis_3body.jl
View File

@ -29,9 +29,9 @@ println("Constructing KE matrices")
println("Constructing PE matrices") println("Constructing PE matrices")
V_l(l, n1, n2) = Va * V_Gaussian(Ra, l, n1, n2; ω=ω) + Vb * V_Gaussian(Rb, l, n1, n2; ω=ω) V_l(l, n1, n2) = Va * V_Gaussian(Ra, l, n1, n2; ω=ω) + Vb * V_Gaussian(Rb, l, n1, n2; ω=ω)
@time "V1" V1 = sp_V_matrix(V_l, n1s, l1s) @time "V1" V1 = sp_V_matrix(V_l, n1s, l1s)
@time "V_relative" V_relative = V1 + sp_V_matrix(V_l, n2s, l2s) @time "V relative" V_relative = V1 + sp_V_matrix(V_l, n2s, l2s)
@time "Moshinsky brackets" Moshinsky_mat = Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ) @time "Moshinsky brackets" U = Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ)
@time "V2 via transform" V2 = Moshinsky_mat * V_relative @time "V2" V2 = U' * V_relative * U
println("Calculating spectrum") println("Calculating spectrum")
@time "H" H = T1 + T2 + V1 + V2 @time "H" H = T1 + T2 + V1 + V2