How can i use GPU Computing for Besselk & Besseli function and integral ？

The form of my function and its integral ：

Function = @( xx_in , a , b ,u, k ) besselk( 0 , sqrt( ( a - xx_in ) .^ 2 + b  ) .* sqrt( u ) ) + k .* besseli( 0 , sqrt( ( a - xx_in ) .^ 2 + b  ) .* sqrt( u ) )  ;   
intmyfun( ixxx , : ) = integral( @( xx_in ) Function( xx_in , A_Fis_roll , A_Comb_roll_r ,  u_f_roll , ksi_2_roll ) , Limt_U_Lock ,    Limt_D_Lock , 'ArrayValued' , true ) ;   %*The parameter ：a , b ,u , k is (1 * 5.8e6 )double 

The Besselk and Besseli is common functions in study of fulid or gas seepage mechanism ，if i have to use GPU computing in matlab to make my code run faster ，Is there any alternative methods？

Thanks a lot！！！

Answer 1

[Edric Ellis]

Unfortunately, the Bessel functions are not yet supported on the GPU, so there is no straightforward workaround for this part. Instead of integral, you might be able to use trapz instead.

[Wang]

Thank you very much for your answer!