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