#10803: critical bug in real_roots
--------------------------------------+-------------------------------------
       Reporter:  zimmerma            |         Owner:  jason, jkantor
           Type:  defect              |        Status:  needs_review  
       Priority:  critical            |     Milestone:  sage-5.1      
      Component:  numerical           |    Resolution:                
       Keywords:  real roots, sd40.5  |   Work issues:                
Report Upstream:  N/A                 |     Reviewers:                
        Authors:  Paul Zimmermann     |     Merged in:                
   Dependencies:                      |      Stopgaps:                
--------------------------------------+-------------------------------------

Comment (by zimmerma):

 I believe I've found the reason for the bug: the function
 {{{down_degree}}} decreases the
 degree of the (interval) polynomial, but keeps the value {{{self.usign}}}.
 With the example in the description and {{{seed=1}}}, {{{down_degree}}} is
 called with:
 {{{
 enter down_degree, self= <IBP: ((34519150034, 25826190301, 17130808432,
 9120071648, 2724478006, -830629355, 3694263) + [0 .. 21299)) * 2^-19 over
 [3921/16384 .. 8017/32768]; level 2; slope_err [-136.91992834614365 ..
 136.91992834614365]> usign= 1
 }}}
 and the polynomial is reduced to:
 {{{
 coeffs= (67412088, 41929949, 16504578, -2435470, -4640) n= 5
 }}}
 where it is clear that {{{usign=1}}} is wrong, since the sign of
 {{{-4640}}} is -1.

 Paul

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

Reply via email to