#244 Taylor internal error (rat problem?)

closed
nobody
6
2007-09-02
2003-03-03
No

taylor(asin(sin(x)+e),x,0,2) =>
Quotient by a polynomial of higher degree

taylor(log(sqrt(e*%e^x+1)+e),x,0,2) =>
Quotient by a polynomial of higher degree

Maxima version: 5.9.0
Maxima build date: 19:10 2/9/2003
host type: i686-pc-mingw32
lisp-implementation-type: Kyoto Common Lisp
lisp-implementation-version: GCL-2-5.0

Discussion

  • Stavros Macrakis

    Logged In: YES
    user_id=588346

    Fateman suggested:

    > Sometimes changing the GCD algorithm to one without
    > a bug in it helps. gcd:subres for example.

    Thanks for the tip. Here are results with the different GCD
    algorithms. Subres/oldgcd is the default in 5.9.0, though the
    documentation says the default is Spmod/zgcd. The doc
    also claims that subres is the "new algorithm", though
    internally, confusingly, the Subres option calls oldgcd; has
    the mapping between names and functions changed?

    --------------

    taylor(asin(sin(x)+e),x,0,2),gcd:XXX;

    EZ OK =ezgcd2
    SPMOD OK =zgcd
    FALSE OK
    MOD OK =newgcd (MOD option in code, not documented)

    SUBRES error =oldgcd
    RED error =oldgcd

    EEZ "EEZGCD is invalid as a function"
    SPHEN "SPHGCD is invalid as a function" (SPHEN option
    in code, not documented)

     
  • Stavros Macrakis

    • priority: 5 --> 6
     
  • Stavros Macrakis

    Logged In: YES
    user_id=588346

    gcd:'algebraic (which calls subres internally) also fails.

     
  • Robert Dodier

    Robert Dodier - 2006-04-09
    • labels: --> Lisp Core - Taylor
     
  • Dan Gildea

    Dan Gildea - 2007-07-31

    Logged In: YES
    user_id=1797506
    Originator: NO

    see 1606731

     
  • Dan Gildea

    Dan Gildea - 2007-09-02
    • status: open --> closed
     
  • Dan Gildea

    Dan Gildea - 2007-09-02

    Logged In: YES
    user_id=1797506
    Originator: NO

    fixed by 1606731

     

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