use math instead of np.math

irnl_rdt_correction/equation_system.py

@@ -7,6 +7,7 @@
 
 """
 import logging
+import math
 from typing import Dict, List, Sequence, Set, Tuple
 
 import numpy as np
@@ -354,7 +355,7 @@ def get_elements_integral(rdt: RDT, ip: int, optics: Optics, feeddown: int) -> f
             iksl_err = 1j*errors_df.loc[elements, f"K{n_mad:d}SL"]
 
             k_sum += ((kl_opt + kl_err + iksl_opt + iksl_err) *
-                      (dx_idy**q) / np.math.factorial(q))
+                      (dx_idy**q) / math.factorial(q))
 
         # note the minus sign before the sum!
         integral += -sum(np.real(i_pow(lm) * k_sum.to_numpy()) * (side_sign * betax * betay).to_numpy())