Hi.
I would to solve an optical waveguide problem, in which I need to numerically calculate multiple times Bessel functions with complex argument. The corresponding code works, but it is quite slow since it takes Mathematica ages to numerically get the value of the Bessel function.
The following code ilustrates the problem :
Table[N[BesselJ[m, x + I]], {m, 100}, {x, 100}]; // AbsoluteTiming
This takes about 21 second, using Matlab below 1 second with identical precision. It there a way to speed this up to reasonable times (i.e., speed up the Bessel function calculation (incl. complex argument))? I found the following link, but I' m not an expert in C, so I' m not sure what to do : https://mathematica.stackexchange.com/questions/84982/how-to-improve-performance-of-besselj-to-the-level-of-gsl
Thanks in advance,
Markus