Cumulative Matrix Multiplication

This function takes a matrix A(n,n,m) and calculates the cumulative matrix product B(n,n)

This C function does a cumulative matrix multiplication of square 2D complex matrices. The matrices need to be put in the function as a 3D matrix of size n x n x m where n is the rank of the square matrices and m is the number of matrices.
I could make a version for real matrices too. If you need this message me.
This function should be compiled using the Mex function of Matlab.

The code is faster when the rank of the matrix is hard-coded in the C function. I have done this for ranks 2, 3 and 4.

% Example of use in Matlab
n = 3 % rank of the square matrices;
m = 1000 % number of matrices;
B = rand(n,n,m); % example input matrices
A = cMn(B, n); % A is a n x n matrix