diff --git a/ho_basis_3body.jl b/ho_basis_3body.jl index 9afa77f..8e20791 100644 --- a/ho_basis_3body.jl +++ b/ho_basis_3body.jl @@ -11,9 +11,8 @@ Va = -2 Ra = 2 μ1 = m * 1/2 -μ2 = m * 2/3 +μ2 = m * 2/3 * 4 c = sqrt(2) -c2 = 2 println("No of threads = ", Threads.nthreads()) @@ -27,7 +26,7 @@ println("Basis size = ", length(Es)) println("Constructing KE matrices") @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") V1_elem(l, n1, n2) = Va * V_Gaussian(Ra, l, n1, n2; ω=ω)