From: Ian Jackson <ijackson@chiark.greenend.org.uk>
Date: Sun, 26 Nov 2017 14:38:29 +0000 (+0000)
Subject: symbolic: actually fix sinc
X-Git-Url: https://www.chiark.greenend.org.uk/ucgi/~ian/git?a=commitdiff_plain;h=64748dfc33b41daf25a6abb495506e3fc5e93e32;p=moebius3.git

symbolic: actually fix sinc

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

diff --git a/symbolic.py b/symbolic.py
index 611bbd5..25c4a68 100644
--- 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)