fix for new sympy version

Author: harrisonlabollita

dosrixs/dosrixs.py

@@ -81,10 +81,11 @@ def build_core_states(edge:EDGE) -> list[YlmExpansion]:
     loop = int(j*2)
     core_states:list[YlmExpansion] = []
     for tmj in [x for x in np.arange(-loop, loop+1, 2) if x != 0]:
-        mj= 0.5*tmj
+        mj= float(0.5*tmj)
         data = {}
         for m in range(-l, l+1):
             for spin in [-s, s]:
+                print(l, m, s, spin, j, mj)
                 data[(m, spin2idx(spin))] = gaunt_sympy(l, m, s, spin, j, mj)
         core_states.append(YlmExpansion(l=l, data=data))
     return core_states

dosrixs/utils.py

@@ -55,4 +55,4 @@ def gaunt(l1:int=1, l2:int=1, l3:int=2, m1:int=0, m2:int=0, m3:int=0) -> float:
 def gaunt_sympy(l1:float, m1:float, l2:float, m2:float, l3:float, m3:float) -> complex:
     from sympy.physics.quantum.cg import CG
     """wrapper around sympy.physics.quantum.cg"""
-    return complex(CG(l1, m1, l2, m2, l3, m3).doit())
+    return complex(CG(l1, m1, l2, m2, l3, m3).doit())