From 073972d0eb61e812ef6c48403348d75a7aef9d46 Mon Sep 17 00:00:00 2001 From: Nuwan Yapa Date: Tue, 12 Mar 2024 19:41:16 -0400 Subject: [PATCH] Critical error in Moshinsky transform fixed --- ho_basis_3body.jl | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/ho_basis_3body.jl b/ho_basis_3body.jl index 70f9fc7..0eb3d4e 100644 --- a/ho_basis_3body.jl +++ b/ho_basis_3body.jl @@ -29,9 +29,9 @@ println("Constructing KE matrices") println("Constructing PE matrices") 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 "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 "V2 via transform" V2 = Moshinsky_mat * V_relative +@time "V relative" V_relative = V1 + sp_V_matrix(V_l, n2s, l2s) +@time "Moshinsky brackets" U = Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ) +@time "V2" V2 = U' * V_relative * U println("Calculating spectrum") @time "H" H = T1 + T2 + V1 + V2