• Woden
  • New member
  • Offline

Topic: Complex arguments to Bessel functions?

Hi all,

I'm using simply fortran for a university project and have run into a slight issue; I'm using the intrinsic Bessel function (BESSEL_JN) in my code which have thus far been inputting real arguments. However I now need to feed complex arguments to the function but the intrinsic function will only accept real arguments.

Is there a way to include or add a Bessel function into the library which can handle complexes? Or any other solutions?

Thanks

  • jeff
  • Administrator
  • Offline

Re: Complex arguments to Bessel functions?

The Fortran 2008 standard only specifies that these intrinsics be provided for REAL arguments currently.  However, I believe the SLATEC library provides this functionality.  I should have a version compiled for Simply Fortran that I can provide you via email to get you up and running.  I'll test it out tomorrow morning and get the library to you.

Jeff Armstrong
Approximatrix, LLC
  • Woden
  • New member
  • Offline

Re: Complex arguments to Bessel functions?

That is very greatly appreciated! Thanks very much Jeff.

  • jeff
  • Administrator
  • Offline

Re: Complex arguments to Bessel functions?

I just sent an email to your email address listed in your forum profile.  Let me know if the attached project works out in your case.

Jeff Armstrong
Approximatrix, LLC