Karim BELABAS on Wed, 18 Oct 2000 13:49:41 +0200 (MET DST) |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: bug: (1+O(x^2))/x -> ** not enough memory |
[Xavier:] > with pari-2.0.20 (beta) > > ? (1+O(x^2))/x > %1 = x^-1 + O(x) > > with the CVS version > > ? (1+O(x^2))/x > *** not enough memory > > Looks like it's a 64bits only bug. Bug starts with line 1876 of gen1.c. > At this point, PARI thinks lgef(x) is 2147483648, so the next command > > y = (GEN) gpmalloc(l*sizeof(long)); > > line 330 of gen2.c (with l being lgef(x)) exhaust all the available > memory. It's a very recently introduced bug, meant to correct another one (related to multivariate power series). I thought there was a typo here and that lgef really meant lg, but it didn't (mistook x for y, the power series for the polynomial...). It's corrected. Sorry, Karim. -- Karim Belabas email: Karim.Belabas@math.u-psud.fr Dep. de Mathematiques, Bat. 425 Universite Paris-Sud Tel: (00 33) 1 69 15 57 48 F-91405 Orsay (France) Fax: (00 33) 1 69 15 60 19 -- PARI/GP Home Page: http://www.parigp-home.de/