#16812: use FLINT to speed up Chebyshev T polynomial creation
-------------------------------------+-------------------------------------
Reporter: rws | Owner:
Type: enhancement | Status: needs_work
Priority: major | Milestone: sage-6.4
Component: symbolics | Resolution:
Keywords: flint, speedup | Merged in:
Authors: Ralf Stephan | Reviewers:
Report Upstream: N/A | Work issues:
Branch: | Commit:
u/rws/use_flint_to_speed_up_chebyshev_t_polynomial_creation|
61972b5574632e914e178fe211062d0dd0e35b87
Dependencies: | Stopgaps:
-------------------------------------+-------------------------------------
Comment (by rws):
I would also say there is a bug in mpmath because it cannot evaluate
`chebyshev_T(10^6, 0.1)` without throwing a `mpmath.NoConvergence` while
Sage's recursive algorithm is happy:
{{{
sage: chebyshev_T(10^6, 0.1, algorithm="recursive")
0.636384327171504
sage: timeit('chebyshev_T(10^6, 0.1, algorithm="recursive")')
625 loops, best of 3: 229 µs per loop
}}}
The result agrees with Wolfram to the 9th place which has
`0.636384327247008190135495327481582817644932547700768253402752...`
--
Ticket URL: <http://trac.sagemath.org/ticket/16812#comment:23>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-trac.
For more options, visit https://groups.google.com/d/optout.
