Mercurial > hg > octave-terminal
changeset 5528:480035e1c2b0
[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