test_clebsch_gordan.py 文件源码

def test_clebsch_gordan3():
    j_1 = S(3)/2
    j_2 = S(3)/2
    m = S(3)
    j = S(3)
    m_1 = S(3)/2
    m_2 = S(3)/2
    assert clebsch_gordan(j_1, j_2, j, m_1, m_2, m) == 1


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

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