#1624 jacobi_p & float overflow

open
Barton Willis
5
2012-11-18
2009-03-20
Barton Willis
No

Using SBCL:

(%i1) jacobi_p(300,42, 12,10.1);

Maxima encountered a Lisp error:

arithmetic error FLOATING-POINT-OVERFLOW signalled

Discussion

  • Barton Willis
    Barton Willis
    2009-04-13

    We have jacobi_p(300,42, 12,10.1) = 4.85995022683574b405. I suppose Maxima could automatically switch to big floats when floating point overflow happens, but I don't know
    how that would work.

     
  • Raymond Toy
    Raymond Toy
    2010-03-29

    I don't think Maxima should do anything about this. There's some expectation that floating-point arithmetic should be fast and we're willing to trade that for things like overflow and underflow. Automatically changing to bigfloats could cause a huge slowdown.