def _rho_kq(rho, j, k, q): v = 0j for m1 in arange(-j, j+1): for m2 in arange(-j, j+1): v += (-1)**(j - m1 - q) * clebsch(j, j, k, m1, -m2, q) * rho.data[m1 + j, m2 + j] return v