# Example SGEMM for N = [2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048, 4096, 8192] A = single( rand(N,N) ); B = single( rand(N,N) ); start = clock(); C = A*B; elapsedTime = etime(clock(), start); gFlops = 2*N*N*N / (elapsedTime * 1e+9); disp(sprintf("N = %4d, elapsed Time = %9.6f, GFlops = %9.6f ", ... N, elapsedTime, gFlops) ); endfor