% % Do benchmark % function [name, time] = bm_test(f,rep) % Actual test functions global t; start = cputime; for i = 1:rep if (f==1) name='Matrix inversion (LINPACK/LAPACK)'; rand('seed', 0); bm_x=inv(rand(256)); elseif (f==2) name='Schur decomposition (LAPACK)'; rand('seed', 0); bm_x=schur(rand(128)); elseif (f==3) name='Linear equation (LINPACK/LAPACK)'; rand('seed', 0); bm_A=rand(256); bm_y=ones(256,1); bm_x=bm_A\bm_y; elseif (f==4) name='Linear equation N RHS (LINPACK/LAPACK)'; rand('seed', 0); bm_A=rand(256); bm_y=ones(256,4); bm_x=bm_A\bm_y; elseif (f==5) name='Matrix determinant (LINPACK/LAPACK)'; rand('seed', 0); bm_x=det(rand(256)); elseif (f==6) name='Matrix Multiply (BLAS)'; rand('seed', 0); bm_x=rand(256)*rand(256); elseif (f==7) name='Fourier transforms (FFTPACK)'; rand('seed', 0); bm_x=ifft2(fft2(rand(256))); elseif (f==8) name='for loop'; for i=1:6000;bm_x=i^2;end end end time = (cputime - start)/rep; return endfunction