print
print
print
print"=========================================================="
print
[v, d] = version;
print " Simple Benchmark test of O-Matrix -- DOUBLE Precision"
print " Version : ", v,d
print " Dell Precision 620, Pentium III Xeon, 800MHz, 640MB ram"
print " Windows 2000  SP4"
print " BLAS Kernel : Intel mkl_p3v7.dll ??"
print " Written by Derek O'Connor, Dec 2004"
print " derekroconnor@eircom.net"
print
print " Run Date : ", clock
print"=========================================================="
print

format double "f10.3";
format int "5";

n = atod(input("Enter matrix size (n): "))
A = double(rand(n,n));

print"    Matrix Size : ",   n, type(A),"precision";
print
begin

# Warmup : loads dlls (I hope)

A = A*A;

totalt = 0.0;
tic
	A*A;
      t = toc
	print "Multiply   ",t, "secs  ", 2*n*n*n/(t*1.0e6), "MFlops (Double Precision)"
	totalt = totalt + t;
tic
	[L,U,P] = lu(A);	
	print "LUP Decomp ",toc, "secs"
	totalt = totalt + toc;
tic
	inv(A);	
	print "Invert     ", toc, "secs"
	totalt = totalt + toc;
tic
	[u,s,v] = svd(A);	
	print "SVD        ", toc, "secs"
	totalt = totalt + toc;
tic
	[Q,R] = qr(A);	
	print "QR fact    ", toc, "secs"
	totalt = totalt + toc;
tic
	[e,d] = eigen(A);	
	print "Eigen      ", toc, "secs"
	totalt = totalt + toc;
tic
	s = det(A);	
	print "Determ     ", toc, "secs"
	totalt = totalt + toc;
tic
	r = rank(A);	
	print "Rank       ", toc, "secs"
	totalt = totalt + toc;
tic
	c = cond(A);	
	print "Condition  ", toc, "secs"
	totalt = totalt + toc;
      print "---------------------------"
	print "Total Time ", totalt, "secs"

end
print
# print profile
print
print "========== Finished SimpleBenchD ======================="