%-----------------------------------------------------------------------
% Script File: matBench
%-----------------------------------------------------------------------
disp('    n         Dot         Saxpy       MatVec        Outer        Direct ')
disp('------------------------------------------------------------------------')
for L = [100,200,400]
   A = randn(L,L);
   B = randn(L,L);
   tic; C = MatMatDot(A,B); f1 = toc;
   tic; C = MatMatSax(A,B); f2 = toc;
   tic; C = MatMatVec(A,B); f3 = toc;
   tic; C = MatMatOuter(A,B); f4 = toc;
   tic; C = A*B; f5 = toc;
disp(sprintf('%6.0f %12.8f %12.8f %12.8f %12.8f %12.8f',L,f1,f2,f3,f4,f5))
end
