| modified bessell functions [message #14786] |
Tue, 30 March 1999 00:00  |
Michael Asten
Messages: 53
Registered: March 1999
|
Member |
|
|
Any one know of routines to compute modified bessell functions In(z) and
Kn(z),
where n=0 or 1, and z is a complex argument. The idl library routine
beselli(x,n)
provides In for real arguments only, (but if fed a complex argument it
seems to cheerily ignore the imaginary part without comment!)
Thanks and regards,
Michael Asten
|
|
|
|
| Re: modified bessell functions [message #14960 is a reply to message #14786] |
Fri, 09 April 1999 00:00  |
Michael Asten
Messages: 53
Registered: March 1999
|
Member |
|
|
Just for the archive record, I'm posting an answer to my own question.
I found siutable fortran77 routines which are free to use conditional on
acknowledgement (given below) and translated line-for-line to idl5.2.
As a matter of interest, the runtime comparision for 10000 calls with
complex argument is 11 sec in idl vs 3 sec in Lahey fortran 90 (standard
opt). I'm quite impressed; I expected greater overheads for idl considering
the fortran loops and goto statments were left intact.
Of course, another answer to my question would be to use pv-wave; - the IMSL
library within pv-wave has the desired complex bessel functions built in!
Regards,
Michael Asten
Michael Asten wrote:
> Any one know of routines to compute modified bessell functions In(z) and
> Kn(z),
> where n=0 or 1, and z is a complex argument. The idl library routine
> beselli(x,n)
> provides In for real arguments only, (but if fed a complex argument it
> seems to cheerily ignore the imaginary part without comment!)
>
> Thanks and regards,
> Michael Asten
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
