gpuArray with besselj for non-integral order values

I'm trying to use gpuArray with the besselj function. However, it doesn't work with non-integral values. I get the following error: 'Error using gpuArray/besselj The NU input (first argument) must be a positive integral value.'

The requirement for the non-gpu besselj is that the order must be real but doesn't need to be an integer. Anybody knows why this is not implemented with the gpuArray version?

Thanks.

[sbstn]

Forgot to mention that I'm trying to run something like this: besselj(1.5,ones(1,3,'gpuArray')).

Answer 1

[Joss Knight]

This is a limitation of the NVIDIA libraries MATLAB uses to compute the bessel functions.

[sbstn]

I was specifically trying to evaluate spherical bessel functions using gpuArrays. I haven't been able to find a fast MATLAB implementation for these functions. Since this is the slowest part of the computation in my application then I will have to implement a CUDA version of these. Thanks for your answer.

[Joss Knight]

I took your post as an enhancement request so we'll see whether that eventually filters through to an implementation. You should put in a request to NVIDIA to improve their math libraries too!