Mercurial > hg > octave-max

comparison libcruft/ranlib/genchi.f @ 10103:0e71ead7359d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use CALL XSTOPX instead of STOP in Fortran ranlib routines New call allows Octave's error handler to intercept otherwise fatal errors in Fortran code
author Rik <rdrider0-list@yahoo.com>
date Wed, 13 Jan 2010 22:22:46 -0800
parents df7c57a6639d
children
comparison
equal deleted inserted replaced
10102:8b4e3388a254 10103:0e71ead7359d
37 C .. 37 C ..
38 C .. Executable Statements .. 38 C .. Executable Statements ..
39 IF (.NOT. (df.LE.0.0)) GO TO 10 39 IF (.NOT. (df.LE.0.0)) GO TO 10
40 WRITE (*,*) 'DF <= 0 in GENCHI - ABORT' 40 WRITE (*,*) 'DF <= 0 in GENCHI - ABORT'
41 WRITE (*,*) 'Value of DF: ',df 41 WRITE (*,*) 'Value of DF: ',df
42 STOP 'DF <= 0 in GENCHI - ABORT' 42 CALL XSTOPX ('DF <= 0 in GENCHI - ABORT')
43 43
44 C JJV changed this to call sgamma directly 44 C JJV changed this to call sgamma directly
45 C 10 genchi = 2.0*gengam(1.0,df/2.0) 45 C 10 genchi = 2.0*gengam(1.0,df/2.0)
46 10 genchi = 2.0*sgamma(df/2.0) 46 10 genchi = 2.0*sgamma(df/2.0)
47 RETURN 47 RETURN