#13054: Pernicious bug for algebraic numbers
----------------------------------------------------------------+-----------
       Reporter:  rbeezer                                       |         
Owner:  davidloeffler
           Type:  defect                                        |        
Status:  needs_review 
       Priority:  major                                         |     
Milestone:  sage-5.7     
      Component:  number fields                                 |    
Resolution:               
       Keywords:  sd40.5                                        |   Work 
issues:               
Report Upstream:  Fixed upstream, but not in a stable release.  |     
Reviewers:               
        Authors:  Jeroen Demeyer                                |     Merged 
in:               
   Dependencies:                                                |      
Stopgaps:               
----------------------------------------------------------------+-----------
Changes (by jdemeyer):

  * status:  new => needs_review


Old description:

> PARI's `polred()` returns reducible polynomials.  GP session:
> {{{
> gp> pol = x^4 - 4294967296*x^2 + 54265257667816538374400;
> gp> L = polred(pol);
> gp> factor(L[4])
> %14 =
> [x^2 + 211955648366398871041 2]
> }}}
>
> '''Upstream''': [http://pari.math.u-bordeaux.fr/cgi-
> bin/bugreport.cgi?bug=1395]
>
> This bug causes problem with arithmetic over `QQbar`:
>
> Line 1563 (5.1.beta0) of sage/rings/qqbar.py:
>
> {{{
> rev = parent(best_elt.Mod(pari_poly).modreverse().lift())
> }}}
>
> Here, we are calling `modreverse()` on a element of a subfield. It would
> make a lot of exact linear algebra '''much''' more reliable if this was
> resolved.
>
> '''spkg''':
> [http://boxen.math.washington.edu/home/jdemeyer/spkg/pari-2.5.3.p3.spkg]

New description:

 PARI's `polred()` returns reducible polynomials.  GP session:
 {{{
 gp> pol = x^4 - 4294967296*x^2 + 54265257667816538374400;
 gp> L = polred(pol);
 gp> factor(L[4])
 %14 =
 [x^2 + 211955648366398871041 2]
 }}}

 '''Upstream''': [http://pari.math.u-bordeaux.fr/cgi-
 bin/bugreport.cgi?bug=1395]

 This bug causes problem with arithmetic over `QQbar`:

 Line 1563 (5.1.beta0) of sage/rings/qqbar.py:

 {{{
 rev = parent(best_elt.Mod(pari_poly).modreverse().lift())
 }}}

 Here, we are calling `modreverse()` on a element of a subfield. It would
 make a lot of exact linear algebra '''much''' more reliable if this was
 resolved.

 '''Apply''': [attachment:13054_polredbest.patch]

 '''spkg''':
 [http://boxen.math.washington.edu/home/jdemeyer/spkg/pari-2.5.3.p3.spkg]

--

-- 
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/13054#comment:13>
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?hl=en.
For more options, visit https://groups.google.com/groups/opt_out.


Reply via email to