#20693: Sage crashes when computing newforms
---------------------------------+------------------------
       Reporter:  ehlen          |        Owner:
           Type:  defect         |       Status:  new
       Priority:  critical       |    Milestone:  sage-7.3
      Component:  modular forms  |   Resolution:
       Keywords:                 |    Merged in:
        Authors:                 |    Reviewers:
Report Upstream:  N/A            |  Work issues:
         Branch:                 |       Commit:
   Dependencies:                 |     Stopgaps:
---------------------------------+------------------------

Comment (by ehlen):

 Replying to [comment:16 pbruin]:
 > From the traceback it appears that NTL runs out of primes for its FFT-
 based `ZZX_XGCD` function.  We should catch the NTL error in
 `NumberFieldElement._invert_c_()` and use an alternative extended GCD
 implementation if that happens.

 Indeed, after finally managing to get sage run in gdb on OSX, I get
 {{{
 #4  0x00000001076cc353 in NTL::UseFFTPrime(long) () from
 /Users/stephan/sage/local/lib/libntl.25.dylib
 #5  0x00000001077bfaf8 in
 NTL::zz_pContext::zz_pContext(NTL::INIT_FFT_STRUCT const&, long) ()
    from /Users/stephan/sage/local/lib/libntl.25.dylib
 #6  0x00000001077bfd07 in NTL::zz_p::FFTInit(long) () from
 /Users/stephan/sage/local/lib/libntl.25.dylib
 #7  0x0000000107743dba in NTL::resultant(NTL::ZZ&, NTL::ZZX const&,
 NTL::ZZX const&, long) ()
    from /Users/stephan/sage/local/lib/libntl.25.dylib
 #8  0x0000000107749594 in NTL::XGCD(NTL::ZZ&, NTL::ZZX&, NTL::ZZX&,
 NTL::ZZX const&, NTL::ZZX const&, long) ()
    from /Users/stephan/sage/local/lib/libntl.25.dylib
 #9  0x000000018b253179 in
 
__pyx_f_4sage_5rings_12number_field_20number_field_element_18NumberFieldElement__invert_c_(__pyx_obj_4sage_5rings_12number_field_20number_field_element_NumberFieldElement*,
 NTL::ZZX*, NTL::ZZ*) ()
 }}}

 Is this a bug in NTL? Is it documented anywhere that this might fail? Can
 we know in advance that this will happen so that we don't have to run into
 _invert_c_c()?

--
Ticket URL: <http://trac.sagemath.org/ticket/20693#comment:17>
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 https://groups.google.com/group/sage-trac.
For more options, visit https://groups.google.com/d/optout.

Reply via email to