tgorordo/WignerSymbols.jl
1
0
Fork 0
mirror of https://github.com/tgorordo/WignerSymbols.jl.git synced 2026-06-13 02:02:14 -07:00
Code Issues Projects Releases Packages Wiki Activity Actions 1

bugfix in Racah W sign when different irreps don't match

Browse source
This commit is contained in:
Jutho Haegeman 2017-10-09 15:28:29 +02:00
parent e56bea0cf9
commit cb09b427c9
1 changed files with 5 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

6
src/WignerSymbols.jl
View file

@ -211,7 +211,11 @@ the `jᵢ`s and `J`s are not integer or halfinteger.
racahW(j₁, j₂, J, j₃, J₁₂, J₂₃) = racahW(Float64, j₁, j₂, J, j₃, J₁₂, J₂₃) racahW(j₁, j₂, J, j₃, J₁₂, J₂₃) = racahW(Float64, j₁, j₂, J, j₃, J₁₂, J₂₃)
function racahW(T::Type{<:AbstractFloat}, j₁, j₂, J, j₃, J₁₂, J₂₃) function racahW(T::Type{<:AbstractFloat}, j₁, j₂, J, j₃, J₁₂, J₂₃)
s = wigner6j(T, j₁, j₂, J₁₂, j₃, J, J₂₃) s = wigner6j(T, j₁, j₂, J₁₂, j₃, J, J₂₃)
return isodd(convert(Int, j₁ + j₂ + j₃ + J)) ? -s : s if !iszero(s) && isodd(convert(Int, j₁ + j₂ + j₃ + J))
return -s
else
return s
end
end end