%-----------------------------------------------------------------------
% script: matflops
%-----------------------------------------------------------------------
disp('    n         Dot         Saxpy       MatVec        Outer        Direct ')
disp('------------------------------------------------------------------------')
for L = [100,200,400]
   A = randn(L,L);
   B = randn(L,L);
   flops(0); C = MatMatDot(A,B); f1 = flops;
   flops(0); C = MatMatSax(A,B); f2 = flops;
   flops(0); C = MatMatVec(A,B); f3 = flops;
   flops(0); C = MatMatOuter(A,B); f4 = flops;
   flops(0); C = A*B; f5 = flops;
disp(sprintf('%6.0f %12.0f %12.0f %12.0f %12.0f %12.0f',L,f1,f2,f3,f4,f5))
end
