GPU, Interpn and Meshgrid

17 次查看(过去 30 天)
Paul
Paul 2014-8-20
回答: 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

回答(2 个)

Matt J
Matt J 2014-8-20
编辑: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 个评论
Paul
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
Matt J 2014-8-21
I'd report it as a bug, then.

请先登录,再进行评论。


Jill Reese
Jill Reese 2014-8-27
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?

类别

Help CenterFile Exchange 中查找有关 GPU Computing 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by