| 8 | |
| 9 | The message is also misleading, since for some inputs it retries `mpmath` anyway: |
| 10 | {{{ |
| 11 | sage: chebyshev_T(100001/2, 2) |
| 12 | --------------------------------------------------------------------------- |
| 13 | NoConvergence Traceback (most recent call last) |
| 14 | <ipython-input-34-9c95a5ff4276> in <module>() |
| 15 | ----> 1 chebyshev_T(Integer(100001)/Integer(2), Integer(2)) |
| 16 | |
| 17 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/sage/functions/ort |
| 18 | 557 return self._eval_(n,x) # Let eval methode decide which is best |
| 19 | 558 else: # Consider OrthogonalPolynomial as symbol |
| 20 | --> 559 return super(OrthogonalPolynomial,self).__call__(n,x) |
| 21 | 560 |
| 22 | 561 |
| 23 | |
| 24 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/sage/symbolic/function.so in sage.symbolic.function.BuiltinFp:8126)() |
| 25 | |
| 26 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/sage/symbolic/function.so in sage.symbolic.function.Function() |
| 27 | |
| 28 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/sage/functions/orth |
| 29 | 496 if not is_Expression(args[-1]): |
| 30 | 497 try: |
| 31 | --> 498 return self._evalf_(*args) |
| 32 | 499 except AttributeError: |
| 33 | 500 pass |
| 34 | |
| 35 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/sage/functions/orthogonal_polys.pyc in _evalf_(self, *args, |
| 36 | 651 from sage.libs.mpmath.all import chebyt as mpchebyt |
| 37 | 652 |
| 38 | --> 653 return mpcall(mpchebyt,args[0],args[-1],prec = precision) |
| 39 | 654 |
| 40 | 655 def _maxima_init_evaled_(self, *args): |
| 41 | |
| 42 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/sage/libs/mpmath/utils.so in sage.libs.mpmath.utils.call (sage/libs/ |
| 43 | |
| 44 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/sage/libs/mpmath/ext_main.so in sage.libs.mpmath.ext_main.wrapped_specfun.__call__ (sage/lit_main.c:17447)() |
| 45 | |
| 46 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/mpmath/functions/orthogonal.pyc in chebyt(ctx, n, x, **kwarg |
| 47 | 444 if (not x) and ctx.isint(n) and int(ctx._re(n)) % 2 == 1: |
| 48 | 445 return x * 0 |
| 49 | --> 446 return ctx.hyp2f1(-n,n,(1,2),(1-x)/2, **kwargs) |
| 50 | 447 |
| 51 | 448 @defun_wrapped |
| 52 | |
| 53 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/mpmath/functions/hypergeometric.pyc in hyp2f1(ctx, a, b, c, z, **kwargs) |
| 54 | 248 @defun |
| 55 | 249 def hyp2f1(ctx,a,b,c,z,**kwargs): |
| 56 | --> 250 return ctx.hyper([a,b],[c],z,**kwargs) |
| 57 | 251 |
| 58 | 252 @defun |
| 59 | |
| 60 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/mpmath/functions/hypergeometric.pyc in hyper(ctx, a_s, b_s, z, **kwargs) |
| 61 | 223 elif q == 0: return ctx._hyp1f0(a_s[0][0], z) |
| 62 | 224 elif p == 2: |
| 63 | --> 225 if q == 1: return ctx._hyp2f1(a_s, b_s, z, **kwargs) |
| 64 | 226 elif q == 2: return ctx._hyp2f2(a_s, b_s, z, **kwargs) |
| 65 | 227 elif q == 3: return ctx._hyp2f3(a_s, b_s, z, **kwargs) |
| 66 | |
| 67 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/mpmath/functions/hypergeometric.pyc in _hyp2f1(ctx, a_s, b_s, z, **kwargs) |
| 68 | 442 if absz <= 0.8 or (ctx.isint(a) and a <= 0 and a >= -1000) or \ |
| 69 | 443 (ctx.isint(b) and b <= 0 and b >= -1000): |
| 70 | --> 444 return ctx.hypsum(2, 1, (atype, btype, ctype), [a, b, c], z, **kwargs) |
| 71 | 445 |
| 72 | 446 orig = ctx.prec |
| 73 | |
| 74 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/mpmath/ctx_mp.pyc in hypsum(ctx, p, q, flags, coeffs, z, accurate_small, **kwargs) |
| 75 | 640 mag_dict = {} |
| 76 | 641 zv, have_complex, magnitude = summator(coeffs, v, prec, wp, \ |
| 77 | --> 642 epsshift, mag_dict, **kwargs) |
| 78 | 643 cancel = -magnitude |
| 79 | 644 jumps_resolved = True |
| 80 | |
| 81 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/mpmath/libmp/libhyper.pyc in _hypsum(coeffs, z, prec, wp, epsshift, magnitude_check, **kwargs) |
| 82 | 319 def _hypsum(coeffs, z, prec, wp, epsshift, magnitude_check, **kwargs): |
| 83 | 320 return hypsum_internal(p, q, param_types, ztype, coeffs, z, |
| 84 | --> 321 prec, wp, epsshift, magnitude_check, kwargs) |
| 85 | 322 |
| 86 | 323 return "(none)", _hypsum |
| 87 | |
| 88 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/sage/libs/mpmath/ext_main.so in sage.libs.mpmath.ext_main.hypsum_internal (sage/libs/mpmath/ext_main.c:24850)() |
| 89 | |
| 90 | /usr/local/src/sage-5.13.beta1/local/lib/python2.7/site-packages/sage/libs/mpmath/ext_impl.so in sage.libs.mpmath.ext_impl.MPF_hypsum (sage/libs/mpmath/ext_impl.c:20614)() |
| 91 | |
| 92 | NoConvergence: Hypergeometric series converges too slowly. Try increasing maxterms. |
| 93 | }}} |
| 94 | }}} |