Multiply slices of a 4D array with vectors without nested loops? Is there a function comparable to pagefun for the CPU?

Question

Petter Stefansson 2016-9-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/303201-multiply-slices-of-a-4d-array-with-vectors-without-nested-loops-is-there-a-function-comparable-to-p

编辑： James Tursa 2016-9-15

I want to multiply each 2D matrix slice of a 4D matrix, X, with a number of vectors, stored in Beta. If I was doing it on the GPU a function like pagefun would be perfect for this since it passes on slices of high dimensional arrays so they can be done in parallel. But since that function does not exist on the CPU I’m really struggling to find an efficient way of writing this. Right now I have this horrible nested for-loop:

for a = 1 : 10
    for b = 1 : 245
        Y(:,a,b) = X(:,:,a,b) * Beta(:,b,a);
    end
end
% size(X) =  19000 x 13 x 10 x 245
% size(Beta) = 13 x 245 x 10

Which is as ugly as it is slow. Can anyone give be some tips on how this can be done in a better and faster way without nested loops? (The Y variable is of course pre-allocated before I run this).