Karim Belabas on Sun, 01 Feb 2004 16:13:05 +0100 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: subst bug ? |
* Bill Allombert [2004-02-01 11:35]: > ? subst(a+O(b),a,c) > %1 = c [...] > It looks like a bug. Fixed in CVS. Was due to inconsistent handling of 0 t_SERs and 0 components of polynomials. I fixed a large number of these yesterday and today, but probably more remain [ symptom: look for gcmp0, and check whether it shouldn't be isexactzero instead ]. From the CHANGES files 16- problems with 0 t_SERs: O(y^5)/(1+x) + O(x^3) --> O(x^3) O(y^5)/(1+x+O(x^3)) --> O(y^5) + O(y^5)*x + O(y^5)*x^2 + O(x^3) O(y) / (1+x) --> O(y) O(y) * 1/(1+x) --> O(y) subst(x+O(b),x,a) --> a 17- (0. * x) / x --> 0.E-28 * x 18- content(0.*x) --> 0 Karim. -- Karim Belabas Tel: (+33) (0)1 69 15 57 48 Dep. de Mathematiques, Bat. 425 Fax: (+33) (0)1 69 15 60 19 Universite Paris-Sud http://www.math.u-psud.fr/~belabas/ F-91405 Orsay (France) http://pari.math.u-bordeaux.fr/ [PARI/GP]