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

Close

#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
    Robert Dodier
    2006-07-01

    • labels: --> 460522
     
  • Robert Dodier
    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
    Robert Dodier
    2006-09-09

    • labels: 460522 --> Lisp Core - Polynomials
     
  • Dan Gildea
    Dan Gildea
    2008-10-18

    • status: open --> pending
     
  • Dan Gildea
    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
    Harald Geyer
    2008-10-19

    • status: pending --> closed
     
  • Harald Geyer
    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.