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

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
--