From d90551577dedba16de6ca0e5989bf6c924587a7a Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <yapanuwan@gmail.com>
Date: Wed, 1 May 2024 16:08:18 -0400
Subject: [PATCH] Error in HO wavefunction fixed

---
 berggren_3body.jl | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/berggren_3body.jl b/berggren_3body.jl
index 5c5a283..1d777ae 100644
--- a/berggren_3body.jl
+++ b/berggren_3body.jl
@@ -60,8 +60,8 @@ function get_W_matrix(basis, basis_HO, μ1ω1, μ2ω2=μ1ω1)
         (i1, i2) = Tuple(idx)
         ((k1, w1), j1, (k2, w2), j2) = basis[i1]
         if j1 == l1s[i2] && j2 == l2s[i2]
-            elem1 = 1/sqrt(sqrt(μ1ω1)) * ho_basis(j1, n1s[i2], 1/sqrt(μ1ω1) * k1) * w1
-            elem2 = 1/sqrt(sqrt(μ2ω2)) * ho_basis(j2, n2s[i2], 1/sqrt(μ2ω2) * k2) * w2
+            elem1 = 1/sqrt(sqrt(μ1ω1)) * (-1)^n1s[i2] * ho_basis(j1, n1s[i2], 1/sqrt(μ1ω1) * k1) * w1
+            elem2 = 1/sqrt(sqrt(μ2ω2)) * (-1)^n2s[i2] * ho_basis(j2, n2s[i2], 1/sqrt(μ2ω2) * k2) * w2
             W[idx] = elem1 * elem2
         end
     end