symbolic: actually fix sinc

Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>

diff --git a/symbolic.py b/symbolic.py
index 611bbd575fedc228c05bd14f57fe996ff025bf95..25c4a686ea3836af5d9a63a8fbf3dd21884135cc 100644 (file)
--- a/symbolic.py
+++ b/symbolic.py
@@ -245,7 +245,9 @@ def gen_C():
 def get_python():
   # https://github.com/sympy/sympy/issues/13642
   # "lambdify sinc gives wrong answer!"
+  out = q_sqparm
   sinc_fixed = Function('sinc_fixed')
   implemented_function(sinc_fixed, lambda x: np.sinc(x/np.pi))
+  out = out..subs(sinc,sinc_fixed)
   p = list(map(eval,params))
-  return lambdify(p, q_sqparm)
+  return lambdify(p, out)