#4835: pari starts up without inititializing enough primes?
---------------------+------------------------------------------------------
 Reporter:  cremona  |        Owner:  tbd
     Type:  defect   |       Status:  new
 Priority:  major    |    Milestone:     
Component:  algebra  |   Resolution:     
 Keywords:           |  
---------------------+------------------------------------------------------
Comment (by cremona):

 I have narrowed it down to here:
 {{{
 sage: K.<z> = CyclotomicField(23)
 sage: pK = K.pari_bnf(certify=False, units=True)
 sage: pK.bnfcertify()

   ***   Warning: large Minkowski bound: certification will be VERY long.
 ---------------------------------------------------------------------------
 PariError                                 Traceback (most recent call
 last)

 /home/john/sage-3.2.2.rc1/<ipython console> in <module>()

 /home/john/sage-3.2.2.rc1/local/lib/python2.5/site-
 packages/sage/libs/pari/gen.so in sage.libs.pari.gen._pari_trap
 (sage/libs/pari/gen.c:38663)()

 PariError: not enough precomputed primes, need primelimit ~  (35)

 sage: pari.default("primelimit")
 500000
 }}}
 I uncommented line 6903 in sage/libs/pari/gen.pyx so I know that
 pari_init() was called with maxprime=500000 (the default) exactly once on
 startup and not again.  So the question in the ticket description seems to
 have a negative answer.

 However in a gp session ("pari -gp", so the same version) I get:
 {{{
 ? bnfcertify(bnfinit(polcyclo(23),1))
   *** bnfcertify: Warning: large Minkowski bound: certification will be
 VERY long.
   *** bnfcertify: not enough precomputed primes, need primelimit ~
 9324407.
 }}}
 so perhaps pari really cannot certify this field with having all primes up
 to 9.3 million, and the problem is that the error message report the wrong
 value for some reason when called from within Sage.  This is confirmed by
 continuing the above Sage session like this:
 {{{
 sage: pari.init_primes(10^7)
 sage: pK.bnfcertify()
   ***   Warning: large Minkowski bound: certification will be VERY long.
 }}}
 followed by a long wait but no error message.

-- 
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/4835#comment:1>
Sage <http://sagemath.org/>
Sage - Open Source Mathematical Software: Building the Car Instead of 
Reinventing the Wheel
--~--~---------~--~----~------------~-------~--~----~
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