Changes between Version 1 and Version 2 of Ticket #9706, comment 90


Ignore:
Timestamp:
12/06/13 10:23:07 (8 years ago)
Author:
jdemeyer
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • Ticket #9706, comment 90

    v1 v2  
    1 Also, you should never simply `print` stuff. Either delete those prints or use a [http://docs.python.org/2/library/warnings.html#available-functions Python warning] or get rid of `mpmath` altogether (since it doesn't seem to work well anyway):
     1Also, you should never simply `print` stuff. Either delete those prints or use a [http://docs.python.org/2/library/warnings.html#available-functions Python warning].
    22{{{
    33sage: chebyshev_T(10^6,RealField(256)(2))
     
    661.764019427245793968639371137094247875315949668035854027376792158135922267593e571947
    77}}}
     8
     9The message is also misleading, since for some inputs it retries `mpmath` anyway:
     10{{{
     11sage: chebyshev_T(100001/2, 2)
     12---------------------------------------------------------------------------
     13NoConvergence                             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
     92NoConvergence: Hypergeometric series converges too slowly. Try increasing maxterms.
     93}}}
     94}}}