diff --git a/ho_basis.jl b/ho_basis.jl index bfc0bd9..268e448 100644 --- a/ho_basis.jl +++ b/ho_basis.jl @@ -103,14 +103,14 @@ function Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ) (Elhs, N, L, n, l) = s[i, :] (Erhs, n1, l1, n2, l2) = s[j, :] if Elhs == Erhs && triangle_ineq(L, l, Λ) && triangle_ineq(l1, l2, Λ) - mat[i, j] = (-1)^(n1 + n2 + N + n) * pick_Moshinsky_brakcet(BRAC, n1, l1, n2, l2, N, L, n, l, Λ) + mat[i, j] = (-1)^(n1 + n2 + N + n) * pick_Moshinsky_bracket(BRAC, n1, l1, n2, l2, N, L, n, l, Λ) end end return sparse(mat) end -function pick_Moshinsky_brakcet(BRAC, n1′, l1′, n2′, l2′, n1, l1, n2, l2, Λ) # Efros notation -- don't confuse +function pick_Moshinsky_bracket(BRAC, n1′, l1′, n2′, l2′, n1, l1, n2, l2, Λ) # Efros notation -- don't confuse ϵ = (l1 + l2 - Λ) % 2 NP = (l1′ - l2′ + Λ - ϵ) ÷ 2 MP = (l1′ + l2′ - Λ - ϵ) ÷ 2