share/fftpack5/fortran/tables.f

   1 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
   2 C
   3 C   FFTPACK 5.0
   4 C
   5 C   Authors:  Paul N. Swarztrauber and Richard A. Valent
   6 C
   7 C   $Id$
   8 C
   9 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  10
  11       SUBROUTINE TABLES (IDO,IP,WA)
  12       REAL  WA(IDO,IP-1,2)
  13 C
  14       TPI = 8.*ATAN(1.)
  15       ARGZ = TPI/REAL(IP)
  16       ARG1 = TPI/REAL(IDO*IP)
  17       DO 110 J=2,IP
  18          ARG2 = REAL(J-1)*ARG1
  19          DO 100 I=1,IDO
  20             ARG3 = REAL(I-1)*ARG2
  21             WA(I,J-1,1) = COS(ARG3)
  22             WA(I,J-1,2) = SIN(ARG3)
  23   100    CONTINUE
  24          IF (IP .LE. 5) GO TO 110
  25          ARG4 = REAL(J-1)*ARGZ
  26          WA(1,J-1,1) = COS(ARG4)
  27          WA(1,J-1,2) = SIN(ARG4)
  28   110 CONTINUE
  29       RETURN
  30       END