fprintf(2,'Starting SubReference Benchmark\n');
P=128;
L=[1000 1000 P];
l=[1000 1000];

M=reshape(1:prod(L),L);
fprintf(2,'Test Matrix Dimensions:\t'); disp(L); fflush(1);

tic();
s = sum(M(:));
t1=toc;
fprintf(2,'Flattened Summation time:\t%g\n',t1);

tic();
s = 0;
for i=1:P
    S=M(:,:,i);
    s = s + sum(S(:));
end
t2=toc();
fprintf(2,'Sliced Summation time:\t%g\n',t2);
fprintf(2,'Efficiency is:\t%g%%\n',100*t1/t2);


fprintf(2,'Starting SubAssignment Benchmark');

m=reshape(1:prod(l),l);
fprintf(2,'Test Matrix Dimensions:\t'); disp(l); fflush(1);

tic();
for i=1:P
    M(:,:,i)=m;
end
t3=toc();
fprintf(2,'Sliced setting time:\t%g\n',t3);