a few cleanup fixes to previous commit

Author: hsnyder

src/fftpack.f90

@@ -30,7 +30,7 @@ end subroutine zffti
 
         !> Version: experimental
         !>
-        !> Forward transform of a double complex periodic sequence.
+        !> Forward transform of a complex periodic sequence.
         !> ([Specification](../page/specs/fftpack.html#zfftf))
         pure subroutine zfftf(n, c, wsave)
             import rk
@@ -62,7 +62,7 @@ end subroutine dffti
 
         !> Version: experimental
         !>
-        !> Forward transform of a double real periodic sequence.
+        !> Forward transform of a real periodic sequence.
         !> ([Specification](../page/specs/fftpack.html#dfftf))
         pure subroutine dfftf(n, r, wsave)
             import rk
@@ -94,7 +94,7 @@ end subroutine dzffti
 
         !> Version: experimental
         !>
-        !> Simplified forward transform of a double real periodic sequence.
+        !> Simplified forward transform of a real periodic sequence.
         !> ([Specification](../page/specs/fftpack.html#dzfftf))
         pure subroutine dzfftf(n, r, azero, a, b, wsave)
             import rk
@@ -154,7 +154,7 @@ end subroutine dcosqb
 
     !> Version: experimental
     !>
-    !> Forward transform of a double complex periodic sequence.
+    !> Forward transform of a complex periodic sequence.
     !> ([Specifiction](../page/specs/fftpack.html#fft))
     interface fft
         pure module function fft_rk(x, n) result(result)
@@ -166,7 +166,7 @@ end function fft_rk
 
     !> Version: experimental
     !>
-    !> Backward transform of a double complex periodic sequence.
+    !> Backward transform of a complex periodic sequence.
     !> ([Specifiction](../page/specs/fftpack.html#ifft))
     interface ifft
         pure module function ifft_rk(x, n) result(result)
@@ -178,7 +178,7 @@ end function ifft_rk
 
     !> Version: experimental
     !>
-    !> Forward transform of a double real periodic sequence.
+    !> Forward transform of a real periodic sequence.
     !> ([Specifiction](../page/specs/fftpack.html#rfft))
     interface rfft
         pure module function rfft_rk(x, n) result(result)
@@ -190,7 +190,7 @@ end function rfft_rk
 
     !> Version: experimental
     !>
-    !> Backward transform of a double real periodic sequence.
+    !> Backward transform of a real periodic sequence.
     !> ([Specifiction](../page/specs/fftpack.html#irfft))
     interface irfft
         pure module function irfft_rk(x, n) result(result)

src/fftpack_fft.f90

@@ -2,7 +2,7 @@
 
 contains
 
-    !> Forward transform of a double complex periodic sequence.
+    !> Forward transform of a complex periodic sequence.
     pure module function fft_rk(x, n) result(result)
         complex(kind=rk), intent(in) :: x(:)
         integer, intent(in), optional :: n

src/fftpack_ifft.f90

@@ -2,7 +2,7 @@
 
 contains
 
-    !> Backward transform of a double complex periodic sequence.
+    !> Backward transform of a complex periodic sequence.
     pure module function ifft_rk(x, n) result(result)
         complex(kind=rk), intent(in) :: x(:)
         integer, intent(in), optional :: n

src/fftpack_irfft.f90

@@ -2,7 +2,7 @@
 
 contains
 
-    !> Backward transform of a double real periodic sequence.
+    !> Backward transform of a real periodic sequence.
     pure module function irfft_rk(x, n) result(result)
         real(kind=rk), intent(in) :: x(:)
         integer, intent(in), optional :: n

src/fftpack_rfft.f90

@@ -2,7 +2,7 @@
 
 contains
 
-    !> Forward transform of a double real periodic sequence.
+    !> Forward transform of a real periodic sequence.
     pure module function rfft_rk(x, n) result(result)
         real(kind=rk), intent(in) :: x(:)
         integer, intent(in), optional :: n

test/tstfft.f

@@ -53,11 +53,12 @@ PROGRAM TSTFFT
 
 
 C
-      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
+      USE fftpack_kind
+      IMPLICIT REAL(RK) (A-H,O-Z)
       DIMENSION       ND(10)     ,X(200)     ,Y(200)     ,W(2000)    ,
      1                A(100)     ,B(100)     ,AH(100)    ,BH(100)    ,
      2                XH(200)    ,CX(200)    ,CY(200)
-      DOUBLE COMPLEX         CX         ,CY
+      COMPLEX(RK)         CX         ,CY
       DATA ND(1),ND(2),ND(3),ND(4),ND(5),ND(6),ND(7)/120,54,49,32,4,3,2/
       SQRT2 = SQRT(2.0D0)
       NNS = 7