gpuArray with besselj for non-integral order values

Question

sbstn 2015-11-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/257066-gpuarray-with-besselj-for-non-integral-order-values

评论： Joss Knight 2015-12-3

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 个评论
显示 -1更早的评论隐藏 -1更早的评论

sbstn 2015-11-23

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

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Joss Knight 2015-11-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/257066-gpuarray-with-besselj-for-non-integral-order-values#answer_201009

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

2 个评论
显示无隐藏无

sbstn 2015-11-30

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 2015-12-3

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!

请先登录，再进行评论。