The primary Mingw Octave 2.9.13 binary test

[Tatsuro MATSUOKA]

I have build an primary binary Octave 2.9.13 on Mingw.
I preformed a simple oregoanator test.
The snapshot is

http://www.geocities.jp/tmacchant/Files/mingwsnap070805.png



The code was cited in the end og this mail.

The result was
octave.exe:1> testOregoB2
ans =  1.7188

The same test was performed with MSVC 2.9.13.
octave:1> testOregoB2
ans =  1.7188

My Cygwin binary was
octave:1> testOregoB2
ans =  2.0010         

The results of my cygwin binary is  a little bit slow.
All test were carried out on Dell Inspiron I6000. (Celeron
M 1.3 GHz, 512MB Memoey)

Build Environmnts
mingw 5.1.3, gcc-3.4.5
msys-1.0.10+msysDTK-1.0.1.exe
Flex, bison, and readline were obtained from GNUwin32.
glob-1.0.tar.bz2 was obtained from the octaveworkshop
package.

Only one modification was carried out in config.h after
./configure to avoid defenition confilct.

983: #if !defined(HAVE_SIGSET_T)
984: //typedef int sigset_t;
985: #endif

Excuse me I do not remember what was the conflict.

Note that any other modifications were not performed.
Another note is that install with 'make install-strip' was
failed, so I used make install. 

I think that the build prcess itself in mingw environment
is not now a hard prcess owing to efforts many
contributers.
The biggest problem is the communication with gnuplot.
The low lelvel ploting "___gnuplot_plot__" has no
probelem.
Oridinary plot did not work.
Perhaps Michael Goffioul has been paid much efforts to
solve the plotting problem.

Perhaps I have to modify the gnuplot code.  But it is too
hard for me now.

T. Matsuoka

****code for testOregoB2
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=cputime();
y = lsode ('oregonator_m', x0, t);
cputime()-ts
% plot (t',y(:,1));
ty1=[t' y(:,1)]; ty2=[t' y(:,2)]; ty3=[t' y(:,3)];
__gnuplot_plot__ ty1, ty2, ty3;

***********************


  *******************************************
  Tatsuro MATSUOKA 
  Email : tmacchant_at_yahoo_dot_co_dot_jp
  (_at_ --> @; _dot_ --> . )
  *******************************************


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