#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 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

-- 
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4000#comment:29>
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
-~----------~----~----~----~------~----~------~--~---

Reply via email to