fix: ensure k=0 is complex if it needs to be

steven-murray · May 31, 2020 · c40925b · c40925b
1 parent c4f712b
commit c40925b
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/hankel/hankel.py b/hankel/hankel.py
@@ -197,8 +197,8 @@ def integrand(r):
                         return f(r).imag * safe_power(r, lim_r_pow)
 
                     int_res_cmplx = quad(integrand, 0, np.inf)
-                    ret_0 = (int_res[0] + int_res_cmplx[0]) * int_fac
-                    err_0 = (int_res[1] + int_res_cmplx[1]) * int_fac
+                    ret_0 = (int_res[0] + 1j * int_res_cmplx[0]) * int_fac
+                    err_0 = (int_res[1] + 1j * int_res_cmplx[1]) * int_fac
 
                 else:
                     ret_0 = int_res[0] * int_fac