How to execute several backslash operations in parallel on GPU?

Is is possible several backslash operations in parallel on GPU?

I converted the MATLAB code to CUDA code.

I expected the three backslash operations would be executed in parallel, but It seemed that the code was executed in serial.

function [z]=backslash(A, a, B, b, C, c)
%#codegen
coder.gpu.kernelfun();
z = [A\a;B\b;C\c]
end

I generated a mex file with codegen command.

cfg=coder.gpuConfig('mex')
codegen -config cfg -args {A, b, B, b, C, c} backslash

Answer 1

[Walter Roberson]

blkdiag(A, B, C) \ [a; b; c]

Answer 2

[Joss Knight]

Try using gpuArray/pagefun before you resort to code generation.

[Walter Roberson]

pagefun will, though, require that all the systems are the same size.

The solution I posted requires that all of the right hand sides have the same number of columns, but does not require that the systems be the same size. I figured that requiring them to have the same number of columns was reasonable considering that if it were not the case then the user's function would have failed in the vertcat stage.