Mercurial > hg > octave-lyh

changeset 5528:480035e1c2b0

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2005-11-01 22:41:22 by jwe]
author jwe
date Tue, 01 Nov 2005 22:41:22 +0000
parents ce71ed4c6da6
children 29d855ffcfad
files libcruft/ChangeLog libcruft/ranlib/ignbin.f libcruft/ranlib/ignpoi.f
diffstat 3 files changed, 15 insertions(+), 5 deletions(-) [+]
line wrap: on
line diff
--- a/libcruft/ChangeLog
+++ b/libcruft/ChangeLog
@@ -1,3 +1,7 @@
+2005-11-01  John W. Eaton  <jwe@octave.org>
+
+	* ranlib/ignbin.f, ranlib/ignpoi.f: Avoid arithmetic IF statements.
+
 2005-10-17  John W. Eaton  <jwe@octave.org>
 
 	* lapack/Makefile.in (dlamc1.o pic/dlamc1.o):
--- a/libcruft/ranlib/ignbin.f
+++ b/libcruft/ranlib/ignbin.f
@@ -178,7 +178,8 @@
 C     .. Executable Statements ..
       IF (pp.NE.psave) GO TO 10
       IF (n.NE.nsave) GO TO 20
-      IF (xnp-30.) 150,30,30
+      IF (xnp-30.0.LT.0.0) GOTO 150
+      GO TO 30
 C
 C*****SETUP, PERFORM ONLY WHEN PARAMETERS CHANGE
 C
@@ -259,7 +260,9 @@
       f = 1.0
       r = p/q
       g = (n+1)*r
-      IF (m-ix) 80,120,100
+      IF (m-ix.LT.0) GO TO 80
+      IF (m-ix.EQ.0) GO TO 120
+      GO TO 100
    80 mp = m + 1
       DO 90 i = mp,ix
           f = f* (g/i-r)
@@ -270,7 +273,8 @@
       DO 110 i = ix1,m
           f = f/ (g/i-r)
   110 CONTINUE
-  120 IF (v-f) 170,170,30
+  120 IF (v-f.LE.0) GOTO 170
+      GO TO 30
 C
 C     SQUEEZING USING UPPER AND LOWER BOUNDS ON ALOG(F(X))
 C
@@ -296,7 +300,8 @@
      +    140./f2)/f2)/f2)/f2)/f1/166320.+ (13860.- (462.- (132.- (99.-
      +    140./z2)/z2)/z2)/z2)/z/166320.+ (13860.- (462.- (132.- (99.-
      +    140./x2)/x2)/x2)/x2)/x1/166320.+ (13860.- (462.- (132.- (99.-
-     +    140./w2)/w2)/w2)/w2)/w/166320.)) 170,170,30
+     +    140./w2)/w2)/w2)/w2)/w/166320. .LE. 0.)) GO TO 170
+      GO TO 30
 C
 C     INVERSE CDF LOGIC FOR MEAN LESS THAN 30
 C
--- a/libcruft/ranlib/ignpoi.f
+++ b/libcruft/ranlib/ignpoi.f
@@ -225,7 +225,8 @@
       xx = x*x
       fx = -0.5*xx
       fy = omega* (((c3*xx+c2)*xx+c1)*xx+c0)
-      IF (kflag) 40,40,60
+      IF (kflag.LE.0) GO TO 40
+      GO TO 60
 C
 C     C A S E  B. (START NEW TABLE AND CALCULATE P0 IF NECESSARY)
 C