Initial commit

3 files changed, 154 insertions, 0 deletions

diff --git a/math/ieee/chap1/time/time12.f b/math/ieee/chap1/time/time12.f
new file mode 100644
index 00000000..774c2cc9
--- /dev/null
+++ b/math/ieee/chap1/time/time12.f
@@ -0,0 +1,53 @@
+c-----------------------------------------------------------------------
+c
+c-----------------------------------------------------------------------
+c
+      parameter  (SIZE = 1024, ILOOP = 100)
+      complex  a, w
+      real     breal(SIZE), bimag(SIZE), qbreal(SIZE), qbimag(SIZE)
+c
+      ioutd = i1mach(2)
+      nn = SIZE
+      tpi = 8.*atan(1.)
+      tpion = tpi/float(nn)
+      w = cmplx(cos(tpion),-sin(tpion))
+c
+c generate a**k as test function
+c result to b(i) for modification by dft and idft subroutines and
+c a copy to qb(i) to compare final result with for error difference.
+c
+      a = (.9,.3)
+      breal(1) = 1.0
+      bimag(1) = 0.0
+      qbreal(1) = 1.0
+      qbimag(1) = 0.0
+      do 10 k=2,nn
+        w = a**(k-1)
+        breal(k) = real(w)
+        bimag(k) = aimag(w)
+        qbreal(k) = breal(k)
+        qbimag(k) = bimag(k)
+  10  continue
+c
+c now compute dft, idft, dft, idft, ...
+c first dft is computed specially, in case subroutine needs to be started.
+c
+      call fft842(0, SIZE, breal, bimag)
+      do 25 icount = 1, ILOOP
+         call fft842(1, SIZE, breal, bimag)
+         call fft842(0, SIZE, breal, bimag)
+   25 continue
+      call fft842(1, SIZE, breal, bimag)
+c
+c calculate rms error between b(i) and qb(i).
+c
+      err = 0.
+      do 30 i=1,nn
+        err = err + (breal(i)-qbreal(i))**2
+     *            + (bimag(i)-qbimag(i))**2
+  30  continue
+      err = sqrt(err / float(nn))
+      write (ioutd,9994) ILOOP, err
+ 9994 format(' rms error, after ', i5, ' loops = ', e15.8)
+      stop
+      end
diff --git a/math/ieee/chap1/time/time17.f b/math/ieee/chap1/time/time17.f
new file mode 100644
index 00000000..adcce4c2
--- /dev/null
+++ b/math/ieee/chap1/time/time17.f
@@ -0,0 +1,53 @@
+c-----------------------------------------------------------------------
+c
+c-----------------------------------------------------------------------
+c
+      parameter  (SIZE = 1008, ILOOP = 100)
+      complex  a, w
+      real     breal(SIZE), bimag(SIZE), qbreal(SIZE), qbimag(SIZE)
+c
+      ioutd = i1mach(2)
+      nn = SIZE
+      tpi = 8.*atan(1.)
+      tpion = tpi/float(nn)
+      w = cmplx(cos(tpion),-sin(tpion))
+c
+c generate a**k as test function
+c result to b(i) for modification by dft and idft subroutines and
+c a copy to qb(i) to compare final result with for error difference.
+c
+      a = (.9,.3)
+      breal(1) = 1.0
+      bimag(1) = 0.0
+      qbreal(1) = 1.0
+      qbimag(1) = 0.0
+      do 10 k=2,nn
+        w = a**(k-1)
+        breal(k) = real(w)
+        bimag(k) = aimag(w)
+        qbreal(k) = breal(k)
+        qbimag(k) = bimag(k)
+  10  continue
+c
+c now compute dft, idft, dft, idft, ...
+c first dft is computed specially, in case subroutine needs to be started.
+c
+      call wfta(breal, bimag, SIZE, 0, 0, ierr)
+      do 25 icount = 1, ILOOP
+         call wfta(breal, bimag, SIZE, 1, 1, ierr)
+         call wfta(breal, bimag, SIZE, 0, 1, ierr)
+   25 continue
+      call wfta(breal, bimag, SIZE, 1, 1, ierr)
+c
+c calculate rms error between b(i) and qb(i).
+c
+      err = 0.
+      do 30 i=1,nn
+        err = err + (breal(i)-qbreal(i))**2
+     *            + (bimag(i)-qbimag(i))**2
+  30  continue
+      err = sqrt(err / float(nn))
+      write (ioutd,9994) ILOOP, err
+ 9994 format(' rms error, after ', i5, ' loops = ', e15.8)
+      stop
+      end
diff --git a/math/ieee/chap1/time/time18.f b/math/ieee/chap1/time/time18.f
new file mode 100644
index 00000000..87680554
--- /dev/null
+++ b/math/ieee/chap1/time/time18.f
@@ -0,0 +1,48 @@
+c-----------------------------------------------------------------------
+c
+c-----------------------------------------------------------------------
+c
+      parameter  (SIZE = 1024, ILOOP = 100)
+      complex  w, b(SIZE), qb(SIZE), a
+      common   /aa/ b
+c
+      ioutd = i1mach(2)
+      nn = SIZE
+      tpi = 8.*atan(1.)
+      tpion = tpi/float(nn)
+      w = cmplx(cos(tpion),-sin(tpion))
+c
+c generate a**k as test function
+c result to b(i) for modification by dft and idft subroutines and
+c a copy to qb(i) to compare final result with for error difference.
+c
+      a = (.9,.3)
+      b(1) = (1.,0.)
+      qb(1) = b(1)
+      do 10 k=2,nn
+        b(k) = a**(k-1)
+        qb(k) = b(k)
+  10  continue
+c
+c now compute dft, idft, dft, idft, ...
+c first dft is computed specially, in case subroutine needs to be started.
+c
+      call radix4(5, 1, -1)
+      do 25 icount = 1, ILOOP
+         call radix4(5, 0,  1)
+         call radix4(5, 0, -1)
+   25 continue
+      call radix4(5, 0,  1)
+c
+c calculate rms error between b(i) and qb(i).
+c
+      err = 0.
+      do 30 i=1,nn
+        de = cabs(qb(i)-b(i))
+        err = err + de**2
+  30  continue
+      err = sqrt(err / float(nn))
+      write (ioutd,9994) ILOOP, err
+ 9994 format(' rms error, after ', i5, ' loops = ', e15.8)
+      stop
+      end