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

add test: wigner3j orthogonality

Browse source
This commit is contained in:
Jutho Haegeman 2019-07-29 12:25:35 +02:00
parent 6e0c8bfc0a
commit e5efd21bdf
1 changed files with 24 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

25
test/runtests.jl
View file

@ -27,12 +27,10 @@ end
ind1 = 1
for m1 in -j1:j1, m2 in -j2:j2
ind2 = 1
for j3 in abs(j1-j2):(j1+j2)
for m3 in -j3:j3
@inbounds for j3 in abs(j1-j2):(j1+j2), m3 in -j3:j3
M[ind1,ind2] = clebschgordan(j1,m1,j2,m2,j3,m3)
ind2 += 1
end
end
ind1 += 1
end
@test M'*M ≈ one(M)
@ -57,6 +55,27 @@ end
end
end
@testset "wigner3j: test orthogonality relations" begin
# equivalent to Clebsch-Gordan orthogonality, now test using Float32
for j1 in smalljlist, j2 in smalljlist
d1::Int = 2*j1+1
d2::Int = 2*j2+1
M = zeros(Float32, (d1*d2, d1*d2))
ind2 = 1
for m1 in -j1:j1, m2 in -j2:j2
ind1 = 1
@inbounds for j3 in abs(j1-j2):(j1+j2), m3 in -j3:j3
d3::Int = 2*j3+1
M[ind1,ind2] += sqrt(d3) * wigner3j(Float32, j1, j2, j3, m1, m2, m3)
ind1 += 1
end
ind2 += 1
end
@test M'*M ≈ one(M) # orthogonality relation type 1
@test M*M' ≈ one(M) # orthogonality relation type 2
end
end
# test 6j
#----------
@testset "wigner6j: test orthogonality" begin