Speed test octave for windows

[Tatsuro MATSUOKA]

Report on cygwin build on 2.9.12

At present Cxsparse could not be built for cygwin using
gcc 3.3.3-3.
So at present, all other requireid components, hdf5, glpk,
... were
sucessfully included.

I carried out speed test. (Dell inspriron I6000, Celeron M
1.3GHz, 512MB)
****lutest

2.9.12 cygwin( with atlas, with sse2
Simple left division
ans =  10.036
LU decomposition
ans =  0.23000

2.9.12 MSVC by Michael
Simple left division
ans =  14.774
LU decomposition
ans =  0.20000

2.9.9+ MSVC by Michael
Simple left division
ans =  8.8586
LU decomposition
ans =  0.20239

2.1.73 Cygwin distrbution
Simple left division
ans = 14.563
LU decomposition
ans = 15.666

2.1.50 binary package with old cygwin
Simple left division
ans = 37.914
LU decomposition
ans = 37.910

*****oregonator test
2.9.12 cygwin( with atlas, with sse2
ans =  10.988

2.9.12 MSVC by Michael
ans =  2.8236

2.9.9+ MSVC by Michael
ans =  2.3504

2.1.73 Cygwin distrbution
ans = 7.8460

2.1.50 binary package with old cygwin
ans = 1.7000

I could not build octave 2.9.12 gcc 3.4.4-3 so I used gcc
3.3.3-3.
I hope the slowness from sjlj excwption on cygwin  will be
solved by
the cywin distributer.

T. Matsuoka



********source lutest.m
%
Num=1000; ItNum=10;
id=1:Num;
rand("seed",1);
A=rand(Num)-0.5;
rand("seed",2);
B=rand(Num,ItNum)-0.5;
[L U P]=lu(A);
%
disp("Simple left division");
ts=time();
for k=1:ItNum
  b=B(:,k);
  x=A\b;
end
time()-ts
x1=x;
%
ts=time();
disp("LU decomposition");
ts=time();
for k=1:ItNum
  b=B(:,k);
  c=P*b; y=L\c; x=U\y;
end
time()-ts
x2=x;
plot(id,x1,"o1",id,x2,"+2");

********source testOregoB.m
clear;
## The `oregonator'.
##
## Reference:
##
##   Oscillations in chemical systems.  IV.  Limit cycle
behavior in a
##   model of a real chemical reaction. Richard J. Field
and Richard
##   M. Noyes, The Journal of Chemical Physics, Volume 60
Number 5,
##   March 1974.function dx = oregonator_m (x, t)

function dx = oregonator_m (x, t)

  dx = zeros (3, 1);

  dx(1) = 77.27*(x(2) - x(1)*x(2) + x(1) -
8.375e-06*x(1)^2);
  dx(2) = (x(3) - x(1)*x(2) - x(2)) / 77.27;
  dx(3) = 0.161*(x(1) - x(3));

end

% The test of `oregonator'.
x0 = [ 4; 1.1; 4 ];
%t = [0, logspace (-1, log10(303), 150), logspace
(log10(304),
log10(500), 150)];
%t=linspace(0,500,1000);
t=0:0.5:500;
ts=time();
y = lsode ('oregonator_m', x0, t);
time()-ts
plot (t',y(:,1));




--------------------------------------
Easy + Joy + Powerful = Yahoo! Bookmarks x Toolbar
http://pr.mail.yahoo.co.jp/toolbar/