Bessel function scaling + limited range?

[Kris Kuhlman]

On Feb 4, 2008 3:07 PM, John W. Eaton <jwe at bevo.che.wisc.edu> wrote:

> so if these functions are subject to the ACM copyright, then I don't
> think we can use them in Octave since the "non-commercial use only"
> clause makes the license incompatible with the GPL.
>

So maybe this is an uninformed question or comment, but I will state what
seems to be to be the obvious.  Since Matlab uses the same Amos Fortran
library that Octave does (according to the "help besseli" command), can
there really be a non-commercial-use license attached with it?

Secondly, if Octave uses the same Amos Fortran libraries as Matlab, why does
it's precision only go to 1.0E+30, while Matlab's goes much higher? (this
question is sort of the point of my initial question in this thread, but I
suppose moving Bessel function support over to GSL is a good idea anyway.)

Since it sounds like the move to GSL will not happen soon, could it be that
maybe the Amos library is configured or compiled wrong?

I have attached a small f90 program and its output (named screen) that
drives the f90 amos code I posted the earlier link to.  It gives me correct
results well beyond where Octave stops.  I have tried to do the same with
the amos f77 code in the libcruft directory, but I get junk.  I am not sure
why.

Curious,

Kris


-- 

Kristopher L. Kuhlman
PhD candidate Department of Hydrology & Water Resources
University of Arizona
(520) 621-1380
http://www.u.arizona.edu/~kuhlman/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://www.cae.wisc.edu/pipermail/bug-octave/attachments/20080206/6cda0c2a/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test_amos_f90.f90
Type: text/x-fortran
Size: 977 bytes
Desc: not available
Url : https://www.cae.wisc.edu/pipermail/bug-octave/attachments/20080206/6cda0c2a/attachment.bin 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: screen
Type: application/octet-stream
Size: 4142 bytes
Desc: not available
Url : https://www.cae.wisc.edu/pipermail/bug-octave/attachments/20080206/6cda0c2a/attachment.obj