Michael Somos on Mon, 26 Aug 2002 20:20:01 -0400


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

gp: forbidden multiplication t_INT * t_RFRAC


This is a bug in division which may have been triggered by not reducing
properly a rational fraction. Note that the first example without the
minus sign is okay. The second example with the '1/Ser(...)' is okay also.
However, the third example which first does the 'Ser(...)' and then the
reciprocal triggers the bug.

==============================================================================
Reading GPRC: /u/home/somos/.gprc ...Done.

          GP/PARI CALCULATOR Version 2.2.4 (development CHANGES-1.491)
                 UltraSparc (MicroSparc kernel) 32-bit version
              (readline v2.2 enabled, extended help not available)

                       Copyright (C) 2002 The PARI Group

PARI/GP is free software, covered by the GNU General Public License, and 
comes WITHOUT ANY WARRANTY WHATSOEVER.

Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.

   realprecision = 28 significant digits
   seriesprecision = 16 significant terms
   format = g0.28

parisize = 4000000, primelimit = 500000
? [1/y,1/(y+z)]
%1 = [1/y, 1/(y + z)]
? Ser(%/%[1])
%2 = 1 + y/(y + z)*x + O(x^2)
? 1/%
%3 = 1 - y/(y + z)*x + O(x^2)
? [-1/y,1/(y+z)]
%4 = [-1/y, 1/(y + z)]
? 1/Ser(%/%[1])
%5 = 1 + y/(y + z)*x + O(x^2)
? [-1/y,1/(y+z)]
%6 = [-1/y, 1/(y + z)]
? Ser(%/%[1])
%7 = 1 + y/(-y - z)*x + O(x^2)
? 1/%
  ***   forbidden multiplication t_INT * t_RFRAC.
? quit
Goodbye!
==============================================================================

Shalom, Michael