#4000: Implement QQ['x'] via Flint ZZ['x'] + denominator
--------------------------------------------------------------------------+-
Reporter: malb |
Owner: somebody
Type: enhancement |
Status: needs_review
Priority: major |
Milestone: sage-4.5.2
Component: basic arithmetic |
Keywords:
Author: Sebastian Pancratz, Martin Albrecht |
Upstream: N/A
Reviewer: John Cremona, Martin Albrecht, Alex Ghitza, Harald Schilly |
Merged:
Work_issues: |
--------------------------------------------------------------------------+-
Comment(by schilly):
Hi, william's nagbot reminded me of this here and I installed these last 3
patches on my 4.5.2 RC build. I tried to replicate your timings, but they
are less convincing than I thought:
{{{
sage: S = PolynomialRing(QQ, 'x', implementation="flint")
sage: R = PolynomialRing(QQ, 'y', implementation="NTL")
sage: f = R.random_element(degree=30); timeit('f*f')
625 loops, best of 3: 15.6 µs per loop
sage: g = S.random_element(degree=30); timeit('g*g')
625 loops, best of 3: 14.7 µs per loop
sage: f = R.random_element(degree=300); timeit('f*f')
625 loops, best of 3: 602 µs per loop
sage: g = S.random_element(degree=300); timeit('g*g')
625 loops, best of 3: 350 µs per loop
sage: f = R.random_element(degree=3000); timeit('f*f')
25 loops, best of 3: 15.2 ms per loop
sage: g = S.random_element(degree=3000); timeit('g*g')
25 loops, best of 3: 19.5 ms per loop
sage: f = R.random_element(degree=3000); timeit('f*f')
25 loops, best of 3: 16.3 ms per loop
sage: g = S.random_element(degree=3000); timeit('g*g')
25 loops, best of 3: 14.3 ms per loop
sage: f = R.random_element(degree=30000); timeit('f*f')
5 loops, best of 3: 1.03 s per loop
sage: g = S.random_element(degree=30000); timeit('g*g')
5 loops, best of 3: 1.04 s per loop
sage: f = R.random_element(degree=30000); timeit('f*f')
5 loops, best of 3: 995 ms per loop
sage: g = S.random_element(degree=30000); timeit('g*g')
5 loops, best of 3: 1.09 s per loop
}}}
Maybe I did something wrong?
I also got this doctest error (my RC built w/o any errors) when checking
the rings/polynomial dir.
{{{
File "/scratch/scratch/schilly/sage/sage-4.5.2.rc0/devel/sage-
main/sage/rings/polynomial/polynomial_element.pyx", line 474:
sage: f(x) is f
Expected:
True
Got:
False
}}}
Btw, the doctest failure in the tutorial I reported earlier is fixed.
The system were I did run this is Ubuntu 8.10 32 bit, gcc version 4.3.2
(Ubuntu 4.3.2-1ubuntu12), Intel(R) Core(TM)2 Quad CPU Q9400
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4000#comment:75>
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 post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
