# HG changeset patch # User jwe # Date 1130884882 0 # Node ID 480035e1c2b03e9b997c7ee3cbd5d74619658fe4 # Parent ce71ed4c6da60acfe80b9d8335a52791dfa2fdb8 [project @ 2005-11-01 22:41:22 by jwe] diff -r ce71ed4c6da6 -r 480035e1c2b0 libcruft/ChangeLog --- a/libcruft/ChangeLog Tue Nov 01 22:15:49 2005 +0000 +++ b/libcruft/ChangeLog Tue Nov 01 22:41:22 2005 +0000 @@ -1,3 +1,7 @@ +2005-11-01 John W. Eaton + + * ranlib/ignbin.f, ranlib/ignpoi.f: Avoid arithmetic IF statements. + 2005-10-17 John W. Eaton * lapack/Makefile.in (dlamc1.o pic/dlamc1.o): diff -r ce71ed4c6da6 -r 480035e1c2b0 libcruft/ranlib/ignbin.f --- a/libcruft/ranlib/ignbin.f Tue Nov 01 22:15:49 2005 +0000 +++ b/libcruft/ranlib/ignbin.f Tue Nov 01 22:41:22 2005 +0000 @@ -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 diff -r ce71ed4c6da6 -r 480035e1c2b0 libcruft/ranlib/ignpoi.f --- a/libcruft/ranlib/ignpoi.f Tue Nov 01 22:15:49 2005 +0000 +++ b/libcruft/ranlib/ignpoi.f Tue Nov 01 22:41:22 2005 +0000 @@ -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