Mercurial > octave-nkf

comparison liboctave/cruft/ranlib/ignpoi.f @ 19575:04d6d592c614

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
avoid valgrind warning: conditional move depends on uninitialized variable * ignpoif: Explicitly initialize and save local array PP.
author John W. Eaton <jwe@octave.org>
date Wed, 07 Jan 2015 17:26:50 -0500
parents 648dabbb4c6b
children
comparison
equal deleted inserted replaced
19574:146f1daa805b 19575:04d6d592c614
97 INTRINSIC abs,alog,exp,float,ifix,max0,min0,sign,sqrt 97 INTRINSIC abs,alog,exp,float,ifix,max0,min0,sign,sqrt
98 C .. 98 C ..
99 C JJV added this for case: mu unchanged 99 C JJV added this for case: mu unchanged
100 C .. Save statement .. 100 C .. Save statement ..
101 SAVE s, d, l, ll, omega, c3, c2, c1, c0, c, m, p, q, p0, 101 SAVE s, d, l, ll, omega, c3, c2, c1, c0, c, m, p, q, p0,
102 + a0, a1, a2, a3, a4, a5, a6, a7, fact, muprev, muold 102 + a0, a1, a2, a3, a4, a5, a6, a7, fact, pp, muprev, muold
103 C .. 103 C ..
104 C JJV end addition - I am including vars in Data statements 104 C JJV end addition - I am including vars in Data statements
105 C .. Data statements .. 105 C .. Data statements ..
106 C JJV changed initial values of MUPREV and MUOLD to -1.0E37 106 C JJV changed initial values of MUPREV and MUOLD to -1.0E37
107 C JJV if no one calls IGNPOI with MU = -1.0E37 the first time, 107 C JJV if no one calls IGNPOI with MU = -1.0E37 the first time,
108 C JJV the code shouldn't break 108 C JJV the code shouldn't break
109 DATA muprev,muold/-1.0E37,-1.0E37/ 109 DATA muprev,muold/-1.0E37,-1.0E37/
110 DATA a0,a1,a2,a3,a4,a5,a6,a7/-.5,.3333333,-.2500068,.2000118, 110 DATA a0,a1,a2,a3,a4,a5,a6,a7/-.5,.3333333,-.2500068,.2000118,
111 + -.1661269,.1421878,-.1384794,.1250060/ 111 + -.1661269,.1421878,-.1384794,.1250060/
112 DATA fact/1.,1.,2.,6.,24.,120.,720.,5040.,40320.,362880./ 112 DATA fact/1.,1.,2.,6.,24.,120.,720.,5040.,40320.,362880./
113 DATA pp/35*0.0/
113 C .. 114 C ..
114 C .. Executable Statements .. 115 C .. Executable Statements ..
115 116
116 IF (mu.EQ.muprev) GO TO 10 117 IF (mu.EQ.muprev) GO TO 10
117 IF (mu.LT.10.0) GO TO 120 118 IF (mu.LT.10.0) GO TO 120