the kind used by fftpack can now be chosen at build time by modifying the defintion of 'rk' in 'src/rk.f'

Author: hsnyder

src/cfftb1.f

@@ -1,5 +1,6 @@
       SUBROUTINE CFFTB1 (N,C,CH,WA,IFAC)
-      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
+      USE fftpack_kind
+      IMPLICIT REAL(RK) (A-H,O-Z)
       DIMENSION       CH(*)      ,C(*)       ,WA(*)      ,IFAC(*)
       NF = IFAC(2)
       NA = 0

src/cfftf1.f

@@ -1,5 +1,6 @@
       SUBROUTINE CFFTF1 (N,C,CH,WA,IFAC)
-      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
+      USE fftpack_kind
+      IMPLICIT REAL(RK) (A-H,O-Z)
       DIMENSION       CH(*)      ,C(*)       ,WA(*)      ,IFAC(*)
       NF = IFAC(2)
       NA = 0

src/cffti1.f

@@ -1,5 +1,6 @@
       SUBROUTINE CFFTI1 (N,WA,IFAC)
-      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
+      USE fftpack_kind
+      IMPLICIT REAL(RK) (A-H,O-Z)
       DIMENSION       WA(*)      ,IFAC(*)    ,NTRYH(4)
       DATA NTRYH(1),NTRYH(2),NTRYH(3),NTRYH(4)/3,4,2,5/
       NL = N

src/cosqb1.f

@@ -1,5 +1,6 @@
       SUBROUTINE COSQB1 (N,X,W,XH)
-      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
+      USE fftpack_kind
+      IMPLICIT REAL(RK) (A-H,O-Z)
       DIMENSION       X(1)       ,W(1)       ,XH(1)
       NS2 = (N+1)/2
       NP2 = N+2

src/cosqf1.f

@@ -1,5 +1,6 @@
       SUBROUTINE COSQF1 (N,X,W,XH)
-      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
+      USE fftpack_kind
+      IMPLICIT REAL(RK) (A-H,O-Z)
       DIMENSION       X(1)       ,W(1)       ,XH(1)
       NS2 = (N+1)/2
       NP2 = N+2

src/dcosqb.f

@@ -1,5 +1,6 @@
       SUBROUTINE DCOSQB (N,X,WSAVE)
-      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
+      USE fftpack_kind
+      IMPLICIT REAL(RK) (A-H,O-Z)
       DIMENSION       X(*)       ,WSAVE(*)
       DATA TSQRT2 /2.82842712474619009760D0/
       IF (N-2) 101,102,103

src/dcosqf.f

@@ -1,5 +1,6 @@
       SUBROUTINE DCOSQF (N,X,WSAVE)
-      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
+      USE fftpack_kind
+      IMPLICIT REAL(RK) (A-H,O-Z)
       DIMENSION       X(*)       ,WSAVE(*)
       DATA SQRT2 /1.41421356237309504880D0/
       IF (N-2) 102,101,103

src/dcosqi.f

@@ -1,5 +1,6 @@
       SUBROUTINE DCOSQI (N,WSAVE)
-      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
+      USE fftpack_kind
+      IMPLICIT REAL(RK) (A-H,O-Z)
       DIMENSION       WSAVE(1)
       DATA PIH /1.57079632679489661923D0/
       DT = PIH/FLOAT(N)

src/dcost.f

@@ -1,5 +1,6 @@
       SUBROUTINE DCOST (N,X,WSAVE)
-      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
+      USE fftpack_kind
+      IMPLICIT REAL(RK) (A-H,O-Z)
       DIMENSION       X(*)       ,WSAVE(*)
       NM1 = N-1
       NP1 = N+1

src/dcosti.f

@@ -1,5 +1,6 @@
       SUBROUTINE DCOSTI (N,WSAVE)
-      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
+      USE fftpack_kind
+      IMPLICIT REAL(RK) (A-H,O-Z)
       DIMENSION       WSAVE(1)
       DATA PI /3.14159265358979323846D0/
       IF (N .LE. 3) RETURN