Learn how easy it is to sync an existing GitHub or Google Code repo to a SourceForge project!

## #195 partfrac, keepfloat: true ERROR

closed
nobody
5
2008-10-19
2002-12-08
Stavros Macrakis
No

partfrac(1/(x-1.0)/(x+1.0),x),keepfloat:true;

gives

Error: Out of bignum stack space

## Discussion

• Robert Dodier
2006-07-01

• labels: --> 460522

• Robert Dodier
2006-09-09

Logged In: YES
user_id=501686

Observed in Maxima 5.9.3cvs. Error message differs depending
on Lisp implementation (all on Linux):

--------------------------
Clisp 2.38:
(%i1) partfrac(1/(x-1.0)/(x+1.0),x),keepfloat:true;
Inverse of zero divisor?

--------------------------
SBCL 0.9.9:
(%i1) partfrac(1/(x-1.0)/(x+1.0),x),keepfloat:true;

Maxima encountered a Lisp error:

Error during processing of --eval option "(cl-user::run)":

The value 1.0 is not of type FIXNUM.

--------------------------
GCL 2.6.7:
(%i1) partfrac(1/(x-1.0)/(x+1.0),x),keepfloat:true;

Factor ran out of primes.

--------------------------
SBCL, GCL, and Clisp all give same result when keepfloat =
false:

(%i6) partfrac (1 / (x - 1.0) / (x + 1.0), x), keepfloat :
false;

`rat' replaced -1.0 by -1/1 = -1.0

`rat' replaced 1.0 by 1/1 = 1.0
(%o6) 1/(2*(x-1))-1/(2*(x+1))

• Robert Dodier
2006-09-09

• labels: 460522 --> Lisp Core - Polynomials

• Dan Gildea
2008-10-18

• status: open --> pending

• Dan Gildea
2008-10-18

OK now with sbcl:

(%i3) partfrac(1/(x-1.0)/(x+1.0),x),keepfloat:true;
(%o3) 0.5/(x-1)-0.5/(x+1)

Can someone check clisp and gcl?

• Harald Geyer
2008-10-19

• status: pending --> closed

• Harald Geyer
2008-10-19

I can confirm that it work's as expected on clisp and gcl too. (Using some 5.16.0cvs version)

Thus closing this report.