fix import/export of ierode common on Windows
[scilab.git] / scilab / modules / differential_equations / src / fortran / rchek.f
index 6257230..c9d762f 100644 (file)
@@ -1,5 +1,7 @@
       subroutine rchek (job, g, neq, y, yh, nyh, g0, g1, gx, jroot, irt)
 clll. optimize
+      include 'stack.h'
+      
       external g
       integer job, neq, nyh, jroot, irt
       double precision y, yh, g0, g1, gx
@@ -21,8 +23,6 @@ clll. optimize
      5   maxord, maxcor, msbp, mxncf, n, nq, nst, nfe, nje, nqu
       common /lsr001/ rownr3(2), t0, tlast, toutc,
      1   iownd3(3), iownr3(2), irfnd, itaskc, ngc, nge      
-      integer         iero
-      common /ierode/ iero
 
 
 c!purpose
@@ -74,7 +74,7 @@ c evaluate g at initial t, and check for zero values. ------------------
  100  continue
       t0 = tn
       call g (neq, t0, y, ngc, g0)
-      if(iero.gt.0) return
+      if(ierror.gt.0) return
       nge = 1
       zroot = .false.
       do 110 i = 1,ngc
@@ -87,7 +87,7 @@ c g has a zero at t.  look at g at t + (small increment). --------------
       do 120 i = 1,n
  120    y(i) = y(i) + temp2*yh(i,2)
       call g (neq, t0, y, ngc, g0)
-      if(iero.gt.0) return
+      if(ierror.gt.0) return
       nge = nge + 1
       zroot = .false.
       do 130 i = 1,ngc
@@ -106,7 +106,7 @@ c
 c if a root was found on the previous step, evaluate g0 = g(t0). -------
       call intdy (t0, 0, yh, nyh, y, iflag)
       call g (neq, t0, y, ngc, g0)
-      if(iero.gt.0) return
+      if(ierror.gt.0) return
       nge = nge + 1
       zroot = .false.
       do 210 i = 1,ngc
@@ -126,7 +126,7 @@ c g has a zero at t0.  look at g at t + (small increment). -------------
       go to 240
  230  call intdy (t0, 0, yh, nyh, y, iflag)
  240  call g (neq, t0, y, ngc, g0)
-      if(iero.gt.0) return
+      if(ierror.gt.0) return
       nge = nge + 1
       zroot = .false.
       do 250 i = 1,ngc
@@ -157,7 +157,7 @@ c set t1 to tn or toutc, whichever comes first, and get g at t1. -------
       do 320 i = 1,n
  320    y(i) = yh(i,1)
  330  call g (neq, t1, y, ngc, g1)
-      if(iero.gt.0) return
+      if(ierror.gt.0) return
       nge = nge + 1
 c call roots to search for root in interval from t0 to t1. -------------
       jflag = 0
@@ -166,7 +166,7 @@ c call roots to search for root in interval from t0 to t1. -------------
       if (jflag .gt. 1) go to 360
       call intdy (x, 0, yh, nyh, y, iflag)
       call g (neq, x, y, ngc, gx)
-      if(iero.gt.0) return
+      if(ierror.gt.0) return
       nge = nge + 1
       go to 350
  360  t0 = x