Ilya Zakharevich on Wed, 03 Feb 2010 14:46:23 +0100


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

Re: PARI stable release 2.3.5 (for real now)


On Mon, Feb 01, 2010 at 05:56:27PM -0800, Ilya Zakharevich wrote:
> On Fri, Jan 29, 2010 at 11:09:50PM +0100, Bill Allombert wrote:
> > I put an hopefully final snapshot at
> > <http://pari.math.u-bordeaux.fr/~bill/pari-2.3.5-pre2.tar.gz>

> > Please test!

> ..././pari-2.3.5/src/language/es.c: In function 'os_read':
> ..././pari-2.3.5/src/language/es.c:2835: warning: ignoring return value
>   of 'read', declared with attribute warn_unused_result

> Otherwise fine.

In fact, I compiled GP/PARI now (above is from Math::Pari), and see a
couple of problems:

  *) startup copyright notice mentions 2006;

  *) intnum() is still borken (same as in 3.0):

    ? intnum(Z=[0,-1/2],[[1],-3/2],1/sqrt(Z)*(1+Z)^(-3/2))
    %1 = 7.150393303262919781690330041 + 1018.215229301285883563073358 I

     (I presume this is integral from 0 to infinity...)
 
  *) moreover, putting some printf()s inside, one can see that the
     expression is calculated for negative Z:

    ? intnum(Z=[0,-1/2],[[1],-3/2],if(print(Z),1,1)/sqrt(Z)*(1+Z)^(-3/2))

     [output skipped].  [*]

Thanks anyway,
Ilya


 [*] P.S. I see that some functions do not take expression-sequences,
     as intnum() in the example above.  (This is why I need this ugly
     `if'-hack.)  I thought this was fixed long ago...

     P.P.S. Is there any simpler way to put print() inside intnum()?