help vectorize code for GPU

Dear All,

I have a for loop code snippet which in essence, involves 2D matrix or elementwise multiplication, and 2D matrix addition, for 2D sheets along the 3rd dimension in 3D double arrays. Here is what I mean:

for i=1:M
A(:,:,i)=B(:,:,i)* C(:,:,i);
D(:,:,i)=E(:,:,i).*F(:,:,i);
G(:,:,i)=H(:,:,i)+I(:,:,i);
end

Besides using parfor, is there a way to vectorize this code for GPU use to speed it (A-I are huge matrices).

Thank you, as always, Octavian

Answer 1

[Matt J]

I assume all the arrays shown are gpuArrays. You can use pagefun for creating A

    A=pagefun(@mtimes, B,C);

There is also mtimesx ( Download ) on the File Exchange which optimizes this kind of thing for CPU variables.

The remaining operations are already internally vectorized, both for ordinary CPU array variables and gpuArrays. You can just do

   D=E.*F;
   G=H+I;