Walter Neumann on Mon, 7 Apr 2003 11:01:05 -0400 (EDT)


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: algdep broken in CVS?


There still seems to be a bug in algdep:

Some time ago Karim wrote

>2.2.5 uses PSLQ which either
>  * returns a polynomial which evaluates to 0 at the input precision [ I
>have just changed that. It used to be "at realprecision" which led to
>weird results if you changed the precision after approximating your
>algebraic number. Now PSLQ is insensitive to realprecision
>
>or
>
>  * returns a real number B and guarantees that no polynomial of height
>less than B can evaluate to 0.

The following example contradicts this:

? \p 50
   realprecision = 57 significant digits (50 digits displayed)
?  a=sqrt(2)+sqrt(3)-5^(1/3)
%1 = 1.4362884232652753529760261931717103356444220186443
? algdep(a,12)
%2 = 235*x^12 - 479*x^11 + 361*x^10 - 278*x^9 + 1252*x^8 - 1922*x^7 +
337*x^6 + 204*x^5 - 549*x^4 - 1096*x^3 + 1491*x^2 - 113*x + 1361
? x=a
%3 = 1.4362884232652753529760261931717103356444220186443
? eval(%2)
%4 = 2.171420994059481315 E-37

--walter