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.

1 个评论

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

请先登录,再进行评论。

 采纳的回答

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

2 个评论

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.
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!

请先登录,再进行评论。

更多回答(0 个)

类别

在 帮助中心 和 File Exchange 中查找有关 Bessel functions 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by