#13314: segmentation fault with roots over ComplexField
----------------------------------------------------------------+-----------
       Reporter:  zimmerma                                      |         
Owner:  burcin      
           Type:  defect                                        |        
Status:  needs_review
       Priority:  major                                         |     
Milestone:  sage-5.3    
      Component:  calculus                                      |    
Resolution:              
       Keywords:                                                |   Work 
issues:              
Report Upstream:  Fixed upstream, but not in a stable release.  |     
Reviewers:              
        Authors:  Paul Zimmermann, Jeroen Demeyer               |     Merged 
in:              
   Dependencies:                                                |      
Stopgaps:              
----------------------------------------------------------------+-----------
Changes (by {'newvalue': u'Paul Zimmermann, Jeroen Demeyer', 'oldvalue': u'Paul 
Zimmermann'}):

  * author:  Paul Zimmermann => Paul Zimmermann, Jeroen Demeyer


Old description:

> this was reported to me by Benjamin Dadoun:
> {{{
> sage: CC=ComplexField(128)
> sage: z=polygen(CC, 'z')
> sage: p=340282366920938463463374607431768211457*z^84 -
> 446662175392119146634315708974074167296*z^80 +
> 279189840125316599983095153558720348160*z^76 -
> 110522928999596165272398266482588385280*z^72 +
> 31087466142071580858231009754176552960*z^68 -
> 6606701181888932299640001500276588544*z^64 +
> 1101219301431275611150985571298443264*z^60 -
> 147498447216213728441335387395194880*z^56 +
> 16134143217737051761322600788131840*z^52 -
> 1456690062677843407408776850964480*z^48 +
> 109261914698248556477361418141696*z^44 -
> 6829504668478678587441535778816*z^40 +
> 355736441058167833075415777280*z^36 - 15393873200469425036030115840*z^32
> + 549832294277433090197422080*z^28 - 16038265915694612515651584*z^24 +
> 375931794881970935169024*z^20 - 6911153168471496130560*z^16 +
> 95997158355731169280*z^12 - 947428403886210560*z^8 + 5921977682886976*z^4
> - 17626570956801
> sage: p.roots(ring=CC)
> ---------------------------------------------------------------------------
> RuntimeError                              Traceback (most recent call
> last)
>
> /users/caramel/bdadoun/Desktop/ceval/<ipython console> in <module>()
>
> /usr/local/sage-5.1-linux-64bit-
> ubuntu_12.04_lts-x86_64-Linux/local/lib/python2\
> .7/site-packages/sage/rings/polynomial/polynomial_element.so in
> sage.rings.poly\
> nomial.polynomial_element.Polynomial.roots
> (sage/rings/polynomial/polynomial_el\
> ement.c:31093)()
>
> /usr/local/sage-5.1-linux-64bit-
> ubuntu_12.04_lts-x86_64-Linux/local/lib/python2\
> .7/site-packages/sage/libs/pari/gen.so in sage.libs.pari.gen.gen.polroots
> (sage\
> /libs/pari/gen.c:35318)()
>
> RuntimeError: Segmentation fault
> }}}
>
> '''spkg''':
> [http://boxen.math.washington.edu/home/jdemeyer/spkg/pari-2.5.2.p1.spkg]
> (diff: [attachment:pari-2.5.2.p1.diff])

New description:

 this was reported to me by Benjamin Dadoun:
 {{{
 sage: CC=ComplexField(128)
 sage: z=polygen(CC, 'z')
 sage: p=340282366920938463463374607431768211457*z^84 -
 446662175392119146634315708974074167296*z^80 +
 279189840125316599983095153558720348160*z^76 -
 110522928999596165272398266482588385280*z^72 +
 31087466142071580858231009754176552960*z^68 -
 6606701181888932299640001500276588544*z^64 +
 1101219301431275611150985571298443264*z^60 -
 147498447216213728441335387395194880*z^56 +
 16134143217737051761322600788131840*z^52 -
 1456690062677843407408776850964480*z^48 +
 109261914698248556477361418141696*z^44 -
 6829504668478678587441535778816*z^40 + 355736441058167833075415777280*z^36
 - 15393873200469425036030115840*z^32 + 549832294277433090197422080*z^28 -
 16038265915694612515651584*z^24 + 375931794881970935169024*z^20 -
 6911153168471496130560*z^16 + 95997158355731169280*z^12 -
 947428403886210560*z^8 + 5921977682886976*z^4 - 17626570956801
 sage: p.roots(ring=CC)
 ---------------------------------------------------------------------------
 RuntimeError                              Traceback (most recent call
 last)

 /users/caramel/bdadoun/Desktop/ceval/<ipython console> in <module>()

 /usr/local/sage-5.1-linux-64bit-
 ubuntu_12.04_lts-x86_64-Linux/local/lib/python2\
 .7/site-packages/sage/rings/polynomial/polynomial_element.so in
 sage.rings.poly\
 nomial.polynomial_element.Polynomial.roots
 (sage/rings/polynomial/polynomial_el\
 ement.c:31093)()

 /usr/local/sage-5.1-linux-64bit-
 ubuntu_12.04_lts-x86_64-Linux/local/lib/python2\
 .7/site-packages/sage/libs/pari/gen.so in sage.libs.pari.gen.gen.polroots
 (sage\
 /libs/pari/gen.c:35318)()

 RuntimeError: Segmentation fault
 }}}

 '''spkg''':
 [http://boxen.math.washington.edu/home/jdemeyer/spkg/pari-2.5.2.p1.spkg]
 (diff: [attachment:pari-2.5.2.p1.diff])

 '''apply''' [attachment:13314_polroots_test.patch]

--

Comment:

 Added a regression test ;-)  Needs review.

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