#4000: [with patch, needs work] Implement QQ['x'] via Flint ZZ['x'] +
denominator
------------------------------+---------------------------------------------
Reporter: malb | Owner: somebody
Type: enhancement | Status: new
Priority: major | Milestone: sage-wishlist
Component: basic arithmetic | Keywords:
Reviewer: | Author: Martin Albrecht
Merged: |
------------------------------+---------------------------------------------
Comment(by mhansen):
Replying to [comment:29 spancratz]:
> As said, the above three lines of code extracted from the ``qqbar.py``
doctests still cause a problem for me. I've chased it down for the last
three hours now, and the following code breaks on my setup:
>
> {{{
> sage: R.<x> = QQ[]
> sage: f = 422826864750/4773824138704099*x^18 -
8134231405059/9547648277408198*x^16 + 11311262264874/4773824138704099*x^14
- 12814039341867/4773824138704099*x^12 -
8509019074752/4773824138704099*x^10 + 707815020483605/9547648277408198*x^8
- 1781974116019893/4773824138704099*x^6+
1316925435907659/4773824138704099*x^4 -
1088322011947813/9547648277408198*x^2 - 1/2*x +
1289415905296105/4773824138704099
> sage: g = -76937/62774*x^19 - 30011/62774*x^18 + 144945/31387*x^17 +
174999/62774*x^16 - 377075/31387*x^15 - 354028/31387*x^14 +
929437/62774*x^13 + 983229/62774*x^12 - 725164/31387*x^11 -
984029/31387*x^10 + 945031/62774*x^9 + 1132829/31387*x^8 +
277343/31387*x^7 - 1107925/62774*x^6 - 432756/31387*x^5 - 23909/62774*x^4
+ 202423/31387*x^3 + 167709/31387*x^2 - 10729/31387*x - 47216/31387
> sage: f(g)
> }}}
>
> I've upload a complete log of the session to
[url]http://sage.pastebin.com/m7757deba[/url]. I am happy also re-
implement polynomial composition using FLINT, it should be a lot faster
than the generic code for this anyway. (Idea: To compose F = f/d with G
= g/e, where f, g are in ZZ[] and d, e are integers, first "rescale" F by
1/e --- this method is implemented already --- and then compose the new
polynomial with g. There is a FLINT function for the last part.)
However, I don't know how or where the generic code is implemented in
SAGE.
>
> Sebastian
This does not crash for me on 64-bit linux with both FLINT 1.3.0 and FLINT
1.5.0. You should try the FLINT 1.5.0 spkg at
http://sage.math.washington.edu/home/mhansen/flint-1.5.0.spkg.
--Mike
> sage: f = 422826864750/4773824138704099*x^18 -
8134231405059/9547648277408198*x^16 + 11311262264874/4773824138704099*x^14
- 12814039341867/4773824138704099*x^12 -
8509019074752/4773824138704099*x^10 + 707815020483605/9547648277408198*x^8
- 1781974116019893/4773824138704099*x^6+
1316925435907659/4773824138704099*x^4 -
1088322011947813/9547648277408198*x^2 - 1/2*x +
1289415905296105/4773824138704099
> sage: g = -76937/62774*x^19 - 30011/62774*x^18 + 144945/31387*x^17 +
174999/62774*x^16 - 377075/31387*x^15 - 354028/31387*x^14 +
929437/62774*x^13 + 983229/62774*x^12 - 725164/31387*x^11 -
984029/31387*x^10 + 945031/62774*x^9 + 1132829/31387*x^8 +
277343/31387*x^7 - 1107925/62774*x^6 - 432756/31387*x^5 - 23909/62774*x^4
+ 202423/31387*x^3 + 167709/31387*x^2 - 10729/31387*x - 47216/31387
> sage: f(g)
> }}}
>
> I've upload a complete log of the session to
[url]http://sage.pastebin.com/m7757deba[/url]. I am happy also re-
implement polynomial composition using FLINT, it should be a lot faster
than the generic code for this anyway. (Idea: To compose F = f/d with G
= g/e, where f, g are in ZZ[] and d, e are integers, first "rescale" F by
1/e --- this method is implemented already --- and then compose the new
polynomial with g. There is a FLINT function for the last part.)
However, I don't know how or where the generic code
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4000#comment:30>
Sage <http://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
-~----------~----~----~----~------~----~------~--~---
