Karim BELABAS on Wed, 4 Aug 1999 09:37:07 +0200 (MET DST)


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

Re: [2.0.16] polylog(m,1.0)


[Daniel Duparc:]
> polylog(2,1)
> polylog(2,1.0)
> polylog(3,1.0)
> etc...
> gives the error
> "non positive argument in mplog".
> in version 2.0.16 for Linux 2.0.37
> on an amd k6-2. 
> 
> As far as I remember, the version in use
> in june 1998 was perfectly correct
> (I verified a small library I wrote for
> polylog in MuPAD against this version)

A typo introduced 1 year ago. Types t_REAL were never recognized
by gcmp1 / gcmp_1

  Karim.

*** src/basemath/gen2.c.orig	Fri Jul 30 14:04:36 1999
--- src/basemath/gen2.c	Tue Aug  3 15:52:06 1999
***************
*** 435,441 ****
        if (signe(x) > 0 && expo(x)==0 && x[2]==HIGHBIT)
        {
          long i,lx = lg(x);
!         for (i=2; i<lx; i++)
            if (x[i]) return 0;
          return 1;
        }
--- 435,441 ----
        if (signe(x) > 0 && expo(x)==0 && x[2]==HIGHBIT)
        {
          long i,lx = lg(x);
!         for (i=3; i<lx; i++)
            if (x[i]) return 0;
          return 1;
        }
***************
*** 480,486 ****
        if (signe(x) < 0 && expo(x)==0 && x[2]==HIGHBIT)
        {
          long i,lx = lg(x);
!         for (i=2; i<lx; i++)
            if (x[i]) return 0;
          return 1;
        }
--- 480,486 ----
        if (signe(x) < 0 && expo(x)==0 && x[2]==HIGHBIT)
        {
          long i,lx = lg(x);
!         for (i=3; i<lx; i++)
            if (x[i]) return 0;
          return 1;
        }
__
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://hasse.mathematik.tu-muenchen.de/ntsw/pari/