23077: fix jacobi_P numerics

sagemath · May 26, 2017 · b5c0864 · b5c0864
1 parent cb88513
commit b5c0864
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/sage/functions/orthogonal_polys.py b/src/sage/functions/orthogonal_polys.py
@@ -1942,7 +1942,10 @@ def _evalf_(self, n, a, b, x, **kwds):
         prec = the_parent.precision()
         BF = CBF(prec+5)
         ret = BF(x).jacobi_P(BF(n), BF(a), BF(b))
-        return the_parent(ret)
+        try:
+            return the_parent(ret)
+        except (TypeError, ValueError):
+            return the_parent.complex_field()(ret)
 
 jacobi_P = Func_jacobi_P()