As so many people have asked me the C++ code for calculating besel functions of complex order and complex arguemnt, I now submit them to benifit more people ...
I am trying to use this function. It seems there is a file "cgamma.m" that is missing. Can you help?
Hi again, I needed to use this function intensively. By shifting the series summation into the mex file, one can gain significant speed-up (for my work, about 50x faster).
Many thanks for this.
Concerning the previous comment, it seems that the tolerance can be adjusted in cbesselj.m ('tol'). Still needs some tweaking to allow one to call the function with independent arrays for nu and x (i.e. both not scalar).
It's a very good initiative but a precision of two or so digits is not good enough for my purposes, and I would think this goes for many other users too.
If you compare the result with mathematica which can set the precision digits. This code has very big error. Thanks for sharing though.
Thank you very much for sharing.
Hongxue, thank you for sharing this.