GPU, Interpn and Meshgrid

Question

Paul 2014-8-20

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/151721-gpu-interpn-and-meshgrid

回答： Jill Reese 2014-8-27

Hi,

I am trying to do some interpolation on a 4 dimensional grid, taking advantage of GPU computing. However, when creating my grid with ndgrid, matlab says that only grid generated with meshgrid are valid. This is true only for the GPU version, there is no problem when I don't use a GPU array.

The problem is that meshgrid cannot be 4 dimensional (maximum is 3). Does it mean than interpn cannot be use with GPU arrays? There is probably a way to avoid this (otherwise a gpu version of interpn is imposible to use), but I cannot find it.

Advice would be welcome!

Thanks

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Matt J 2014-8-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/151721-gpu-interpn-and-meshgrid#answer_149320

编辑：Matt J 2014-8-20

I suspect you are inadvertently calling interp2, rather than interpn. interp2 does require that the data be as if produced by meshgrid. In any case, you should show the code you have run so that we can reproduce it.

Ultimately, though, there is no reason to use either ndgrid or meshgrid when working with interpn or griddedInterpolant. It consumes extra memory unnecessarily and is, I believe, purely for backward compatibility. You can just call interpn with the original 1D grid vectors used as arguments to ndgrid. This is true in the CPU version and should be true on the GPU as well.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Paul 2014-8-21

Yes I am using R2014a, such that interpn is gpu-compatible.

I tried using the grids vector as you recommended, but I got the same error message.

Matt J 2014-8-21

I'd report it as a bug, then.

请先登录，再进行评论。

Answer 2

Jill Reese 2014-8-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/151721-gpu-interpn-and-meshgrid#answer_149988

If you look at the help for gpuArray/interpn in R2014a it shows an example of 4D interpolation. My theory is that your inputs to NDGRID don't satisfy some of the restrictions imposed on the grid by the GPU version of interpn, namely that:

The grid arguments must be finite vectors or N-D arrays with increasing and non-repeating elements in corresponding dimension.
When the grid arguments are N-D arrays, they must align with cartesian axes, as if they were produced by ndgrid.

Can you share what your inputs to ndgrid look like?