two - loop with GPU

Question

Anh 2013-4-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/69698-two-loop-with-gpu

I just have GPU on my desktop with matlab R2012b I plotted a Gaussian profile using two - loop on CPUas below:

for n=1:1:N
       for m=1:1:N
            x(n)=n*ds-L/2;
            y(m)=m*ds-L/2;
            G(n,m)=exp(-((x(n)^2)+(y(m)^2)))/(w^2);
       end
end

I read somewhere that GPU does not support index like that, I even tried but it not work either. Is there any way possible to do the same thing on GPU?

Thank so much

2 个评论
显示无隐藏无

James Lebak 2013-4-5

It seems like it should be possible to do this on the GPU. Can you tell us what variables were on the GPU, and what happened when it failed to work? Did you receive an error message, or did MATLAB compute an incorrect answer?

Anh 2013-4-15

All variables were on GPU (N, ds, N). I got the message: Undefined function 'colon' for input arguments of type 'gpuArray' Even one "For loop", if N is on GPU, it not work either. I know there is one way we can deal with "For loop" of an array by using gpuArray.colon.(...), but I have loop inside the loop to create a matrix G(n,m), so I don't know how to solve this problem.

请先登录，再进行评论。

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

Answer 1

James Lebak 2013-4-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/69698-two-loop-with-gpu#answer_83239

在 MATLAB Online 中打开

Yes, 1:gpuArray(N) for scalar N doesn't work in MATLAB R2013a. Try the following:

x=gpuArray.colon(1,N)*ds-L/2;
y=gpuArray.colon(1,N)*ds-L/2;
G = exp(-((x.*x).'*gpuArray.ones(1,N)+gpuArray.ones(N,1)*(y.*y)))./(w^2);

This vectorizes the calculation of x and y and calculates all the elements in the matrix at once, which should be much more efficient on the GPU than using a loop.