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