mttestb.m

%MTTESTB runs iterative timing tests for Matrix Transpose algorithms.

% $Source: /home/stpierre/cvsroot/matlab/simulink/matrix/test/mttestb.m,v $
% $Revision: 1.2 $
% $Date: 2009-07-19 22:03:27 $

% Copyright (c) 2000-2009, Jay A. St. Pierre.  All rights reserved.

stop_time=2e6;
A=[1 2 3; 4 5 6; 7 8 9];
iterations=6;

save mttest_vars stop_time A iterations

mttest1
mttest2
mttest3
disp(' ')

clear, pack
disp('mttest1:')
disp('========')
load mttest_vars
for k=1:iterations
  ts(k)=timesim('mttest1');
  disp(['Time: ', num2str(ts(k))])
end
simAt=rm2mat(simout);
disp(['Min/Mean/Max Time: ', ...
      num2str(min(ts)),  '/', ...
      num2str(mean(ts)), '/', ...
      num2str(max(ts))])
disp(['Standard Deviation for Time: ', num2str(std(ts))])
disp(['Total Time: ', num2str(sum(ts))])
if simAt == A.'
  disp('simAt == A.''  ***PASSED***')
else
  disp('simAt ~= A.''  ***FAILED***')
  simAt
  Ax
end
disp(' ')

clear, pack
disp('mttest2:')
disp('========')
load mttest_vars
for k=1:iterations
  ts(k)=timesim('mttest2');
  disp(['Time: ', num2str(ts(k))])
end
simAt=rm2mat(simout);
disp(['Min/Mean/Max Time: ', ...
      num2str(min(ts)),  '/', ...
      num2str(mean(ts)), '/', ...
      num2str(max(ts))])
disp(['Standard Deviation for Time: ', num2str(std(ts))])
disp(['Total Time: ', num2str(sum(ts))])
if simAt == A.'
  disp('simAt == A.''  ***PASSED***')
else
  disp('simAt ~= A.''  ***FAILED***')
  simAt
  Ax
end
disp(' ')

clear, pack
disp('mttest3:')
disp('========')
load mttest_vars
for k=1:iterations
  ts(k)=timesim('mttest3');
  disp(['Time: ', num2str(ts(k))])
end
simAt=rm2mat(simout);
disp(['Min/Mean/Max Time: ', ...
      num2str(min(ts)),  '/', ...
      num2str(mean(ts)), '/', ...
      num2str(max(ts))])
disp(['Standard Deviation for Time: ', num2str(std(ts))])
disp(['Total Time: ', num2str(sum(ts))])
if simAt == A.'
  disp('simAt == A.''  ***PASSED***')
else
  disp('simAt ~= A.''  ***FAILED***')
  simAt
  Ax
end
disp(' ')