MATLAB Answers

Paul
1

GPU, Interpn and Meshgrid

Asked by Paul
on 20 Aug 2014
Latest activity Answered by Jill Reese
on 27 Aug 2014
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 Comments

Sign in to comment.

2 Answers

Matt J
Answer by Matt J
on 20 Aug 2014
Edited by Matt J
on 20 Aug 2014

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 Comments

Show 1 older comment
Matt J
on 20 Aug 2014
Are you using version R2014? In earlier versions, interpn does not support gpuArray, only interp1 and interp2. It might default back to interp2 in that case, though I admit that would be strange.
Have you tried instead,
DEV_gpu=interpn(lgrid(:),pgrid(:).',ngrid(:),Xgrid(:).',...
target_gpu,xi(:,1),xi(:,2),xi(:,3),xi(:,4),'linear');
If so, what happens?
Paul
on 21 Aug 2014
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
on 21 Aug 2014
I'd report it as a bug, then.

Sign in to comment.


Answer by Jill Reese
on 27 Aug 2014

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?

  0 Comments

Sign in to comment.