Mercurial > hg > octave-terminal
changeset 3118:74cc8e2fe2c0
[project @ 1997-11-29 19:47:53 by jwe]
author | jwe |
---|---|
date | Sat, 29 Nov 1997 19:47:53 +0000 |
parents | f735f3ea1ee7 |
children | f3e1da120048 |
files | libcruft/specfun/ribesl.f libcruft/specfun/rjbesl.f libcruft/specfun/rybesl.f |
diffstat | 3 files changed, 18 insertions(+), 19 deletions(-) [+] |
line wrap: on
line diff
--- a/libcruft/specfun/ribesl.f +++ b/libcruft/specfun/ribesl.f @@ -180,15 +180,18 @@ PARAMETER (ONE = 1.0D0, TWO = 2.0D0, ZERO = 0.0D0) PARAMETER (HALF = 0.5D0, CONST = 1.585D0) C------------------------------------------------------------------- + DATA FIRST /.TRUE./ +C------------------------------------------------------------------- +C Statement functions for conversion +C------------------------------------------------------------------- + CONV(N) = DBLE(N) + FUNC(X) = DGAMMA(X) +C------------------------------------------------------------------- C Machine-dependent parameters C------------------------------------------------------------------- - DATA NSIG,XLARGE,EXPARG /16,1.0D4,709.0D0/ - DATA ENTEN,ENSIG,RTNSIG/1.0D308,1.0D16,1.0D-4/ - DATA ENMTEN/8.9D-308/ - DATA FIRST /.TRUE./ IF (FIRST) THEN NSIG = NINT (-LOG (D1MACH (1))) - ENTEN = 1.0D1 ** (INT (LOG10 (D1MACH (2))) + ENTEN = 1.0D1 ** (INT (LOG10 (D1MACH (2)))) ENSIG = 1.0D1 ** NSIG RTNSIG = 1.0D1 ** (-NINT (NSIG / 4.0)) ENMTEN = 4.0D0 * D1MACH (1) @@ -197,11 +200,6 @@ FIRST = .FALSE. ENDIF C------------------------------------------------------------------- -C Statement functions for conversion -C------------------------------------------------------------------- - CONV(N) = DBLE(N) - FUNC(X) = DGAMMA(X) -C------------------------------------------------------------------- C Check for X, NB, OR IZE out of range. C------------------------------------------------------------------- IF ((NB.GT.0) .AND. (X .GE. ZERO) .AND.
--- a/libcruft/specfun/rjbesl.f +++ b/libcruft/specfun/rjbesl.f @@ -155,7 +155,7 @@ C TWOPI1 + TWOPI2 = 2 * PI to extra precision. C--------------------------------------------------------------------- PARAMETER (PI2 = 0.636619772367581343075535D0) - PARAMETER (TWOPI1 = 6.28125D0, + PARAMETER (TWOPI1 = 6.28125D0) PARAMETER (TWOPI2 = 1.935307179586476925286767D-3) PARAMETER (ZERO = 0.0D0, EIGHTH = 0.125D0, HALF = 0.5D0) PARAMETER (ONE = 1.0D0, TWO = 2.0D0, THREE = 3.0D0) @@ -171,12 +171,18 @@ 4 5.109094217170944D19,1.12400072777760768D21, 5 2.585201673888497664D22,6.2044840173323943936D23/ C--------------------------------------------------------------------- + DATA FIRST /.TRUE./ +C--------------------------------------------------------------------- +C Statement functions for conversion and the gamma function. +C--------------------------------------------------------------------- + CONV(I) = DBLE(I) + FUNC(X) = DGAMMA(X) +C--------------------------------------------------------------------- C Machine-dependent parameters C--------------------------------------------------------------------- - DATA FIRST /.TRUE./ IF (FIRST) THEN NSIG = NINT (-LOG (D1MACH (1))) - ENTEN = 1.0D1 ** (INT (LOG10 (D1MACH (2))) + ENTEN = 1.0D1 ** (INT (LOG10 (D1MACH (2)))) ENSIG = 1.0D1 ** NSIG RTNSIG = 1.0D1 ** (-NINT (NSIG / 4.0)) ENMTEN = 4.0D0 * D1MACH (1) @@ -184,11 +190,6 @@ FIRST = .FALSE. ENDIF C--------------------------------------------------------------------- -C Statement functions for conversion and the gamma function. -C--------------------------------------------------------------------- - CONV(I) = DBLE(I) - FUNC(X) = DGAMMA(X) -C--------------------------------------------------------------------- C Check for out of range arguments. C--------------------------------------------------------------------- MAGX = INT(X)
--- a/libcruft/specfun/rybesl.f +++ b/libcruft/specfun/rybesl.f @@ -158,7 +158,7 @@ PARAMETER (PIBY2 = 1.5707963267948966192D0) PARAMETER (PI = 3.1415926535897932385D0) PARAMETER (SQ2BPI = 7.9788456080286535588D-1) - PARAMETER (PIM5 = /7.0796326794896619231D-1) + PARAMETER (PIM5 = 7.0796326794896619231D-1) PARAMETER (ONBPI = 3.1830988618379067154D-1) C---------------------------------------------------------------------- C Coefficients for Chebyshev polynomial expansion of