test_clebsch_gordan.py 文件源码

def test_clebsch_gordan5():
    j_1 = S(5)/2
    j_2 = S(1)
    m = S(7)/2
    j = S(7)/2
    m_1 = S(5)/2
    m_2 = 1
    assert clebsch_gordan(j_1, j_2, j, m_1, m_2, m) == 1


    j_1 = S(5)/2
    j_2 = S(1)
    m = S(5)/2
    j = S(5)/2
    m_1 = S(5)/2
    m_2 = 0
    assert clebsch_gordan(j_1, j_2, j, m_1, m_2, m) == sqrt(5)/sqrt(7)

    j_1 = S(5)/2
    j_2 = S(1)
    m = S(3)/2
    j = S(3)/2
    m_1 = S(1)/2
    m_2 = 1
    assert clebsch_gordan(j_1, j_2, j, m_1, m_2, m) == 1/sqrt(15)