Karim Belabas on Mon, 06 Aug 2007 12:24:09 +0200 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: bug in current CVS |
* Walter Neumann [2007-08-06 05:33]: > Using the current CVS with GMP: [...] > I get the following: > > ? f(a)=vecsort(subst(nfgaloisconj(x^4+1,,a),x,100)) > ? f(1) > *** nfgaloisconj: bug in PARI/GP (Segmentation Fault), please report Yes, or simply 'nfgaloisconj(x^4+1,,1)'. Igor Schein reported a similar bug as #652 <http://pari.math.u-bordeaux.fr/cgi-bin/bugreport.cgi?bug=652>. This is due to a recent change I made to the gerepileupto() libpari function (call lower-level gerepileuptoint / gerepiluptoleaf if possible), which didn't take into account the widespread use of an undocumented hack : the possibility to use sentinels. Namely, gerepileupto(av, q) was functionally equivalent to avma = av when q was not on the PARI stack. It needn't even be a well formed GEN, and this was exploited here with q = NULL. I have fixed my previous patch, and documented the hack, promoting it to a useful feature. Thanks for spotting this! K.B. -- Karim Belabas Tel: (+33) (0)5 40 00 26 17 Universite Bordeaux 1 Fax: (+33) (0)5 40 00 69 50 351, cours de la Liberation http://www.math.u-bordeaux.fr/~belabas/ F-33405 Talence (France) http://pari.math.u-bordeaux.fr/ [PARI/GP] `