Index exceeds size of array in gpukernelcore_p4.cu?

[helange23]

I am trying to build a pytorch interface based on this implementation and might have found a possible bug?

Note that in order to get the code to work 10 years after it has been written, I have changed some things but mostly only removed the SDK dependent checks but most of the remaining code remains the same.

I have encountered the following problem:

In the function m2l_kernel (line 322) a variable "jnk" is declared on line 429.
jnk=(j+n)*(j+n)+j+n;

jnk is used to access an array of size numExpansion2 (10**2) but jnk oftentimes exceeds 100 which results in an illegal memory access.

Is this a bug or am I doing something wrong?

[rioyokota]

I really don't remember, but I guess numExpansion2 should be 4 * numExpansion * numExpansion.
What exactly are you trying to do in PyTorch with FMM?

[zjin-lcf]

Thank you for fixing the issues.

N = 100000 (iteration = 8)
level : 3
Cuda failure gpukernel_p3.cu:1419: 'an illegal memory access was encountered'

N = 100000 (iteration = 8)
level : 3
Cuda failure gpukernel_p4.cu:733: 'an illegal memory access was encountered'