Cumulative 2x2 Complex Matrix Multiplication

This takes a 2x2xm matrix and calculates the cumulative matrix multiplication along the 3rd diim.

This function calculates the cumulative matrix multiplication of input complex matrices A(2,2,m)
This function should be compiled using the Maltab Mex function.
m = 1000; % number of matrices
A = rand(2,2,m) + 1i * rand(2,2,m); % Example matrix;
B = cM(A); % B is a 2x2 matrix

This function is up to 180 times faster that the equivalent for loop.

B = eye(2);
for ii = 1 : 1000
B = B * A(:,:,ii);
end