#195 partfrac, keepfloat: true ERROR

closed
nobody
5
2008-10-19
2002-12-08
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.

     

Log in to post a comment.

Get latest updates about Open Source Projects, Conferences and News.

Sign up for the SourceForge newsletter:

JavaScript is required for this form.





No, thanks