Merge pull request #3 from loiseaujc/intent_specs

Intent specs

Author: loiseaujc

src/dcost.f90

@@ -2,7 +2,7 @@ subroutine dcost(n, x, wsave)
          use fftpack_kind, only: rk
          implicit none
          integer, intent(in) :: n
-         real(rk), intent(in) :: wsave(*)
+         real(rk), intent(inout) :: wsave(*)
          real(rk), intent(inout) :: x(*)
          real(rk) :: c1, t1, t2, tx2, x1h, x1p3, &
                      xi, xim2

src/dfftb.f90

@@ -3,7 +3,7 @@ subroutine dfftb(n, r, wsave)
          implicit none
          integer, intent(in) :: n
          real(dp), intent(inout) :: r(*)
-         real(dp), intent(in) :: wsave(*)
+         real(dp), intent(inout) :: wsave(*)
          if (n == 1) return
          call rfftb1(n, r, wsave, wsave(n + 1), wsave(2*n + 1))
       end subroutine dfftb

src/dfftf.f90

@@ -3,7 +3,7 @@ subroutine dfftf(n, r, wsave)
          implicit none
          integer, intent(in) :: n
          real(dp), intent(inout) :: r(*)
-         real(dp), intent(in) :: wsave(*)
+         real(dp), intent(inout) :: wsave(*)
          if (n == 1) return
          call rfftf1(n, r, wsave, wsave(n + 1), wsave(2*n + 1))
       end subroutine dfftf

src/dzfftb.f90

@@ -3,7 +3,7 @@ subroutine dzfftb(n, r, azero, a, b, wsave)
          implicit none
          integer, intent(in) :: n
          real(dp), intent(out) :: r(*)
-         real(dp), intent(in) :: wsave(*)
+         real(dp), intent(inout) :: wsave(*)
          real(dp), intent(in) :: azero, a(*), b(*)
          integer :: i, ns2
          if (n < 2) then

src/zfftb.f90

@@ -3,7 +3,7 @@ subroutine zfftb(n, c, wsave)
          implicit none
          integer, intent(in) :: n
          real(dp), intent(inout) :: c(*)
-         real(dp), intent(in) :: wsave(*)
+         real(dp), intent(inout) :: wsave(*)
          integer :: iw1, iw2
          if (n == 1) return
          iw1 = n + n + 1

src/zfftf.f90

@@ -3,7 +3,7 @@ subroutine zfftf(n, c, wsave)
          implicit none
          integer, intent(in) :: n
          real(dp), intent(inout) :: c(*)
-         real(dp), intent(in) :: wsave(*)
+         real(dp), intent(inout) :: wsave(*)
          integer :: iw1, iw2
          if (n == 1) return
          iw1 = n + n + 1