On 3/5/12 3:11 PM, Edwin Woollett wrote:
> One message from slatec is not surpressed when
> using quad_control:
> ----------------------------------------
> (%i1) display2d:false$
>
> (%i2) float(bessel_i(1,10^4));
>
> ***MESSAGE FROM ROUTINE DBESI1 IN LIBRARY SLATEC.
> ***FATAL ERROR, PROG ABORTED, TRACEBACK REQUESTED
> * ABS(X) SO BIG I1 OVERFLOWS
> * ERROR NUMBER = 2
> * ***END OF MESSAGE
I think that's because it's a fatal error and dbesi is expecting the
program to abort. Of course, we don't.
I think it would be better to handle this in the Lisp wrapper and check
for overflows there instead of mucking around with the slatec fortran
routines.
Ray