home *** CD-ROM | disk | FTP | other *** search
/ CP/M / CPM_CDROM.iso / simtel / cpmug / cpmug091.ark / TWIDDLE.FOR < prev    next >
Encoding:
Text File  |  1984-04-29  |  1.4 KB  |  46 lines
  1. C ***********************************************************
  2. C
  3. C THIS PROGRAM GENERATES TWIDDLE FACTORS FOR A
  4. C DATA A STATEMENT OF AN FFT PROGRAM.
  5. C FILENAME IS "TWIDDLE.LIB"
  6. C
  7. C Written by Victor De Pinto
  8. C Revised July 1981
  9. C
  10.       PROGRAM JENIFER
  11.       DOUBLE PRECISION PI, A, A1, A2, A3, A4
  12.       DATA PI, LFEED / 3.141592653590D+00, X'0A' /
  13.       CALL OPEN (3, 'TWIDDLE LIB', 0)
  14.       WRITE (3,10)
  15. 10    FORMAT ()
  16.       DO 1 I=1,512,5
  17.       I1=I+1
  18.       I2=I+2
  19.       I3=I+3
  20.       I4=I+4
  21.       A =DCOS( (I  )*(-PI/512.0D+00) )
  22.       A1=DCOS( (I+1)*(-PI/512.0D+00) )
  23.       A2=DCOS( (I+2)*(-PI/512.0D+00) )
  24.       A3=DCOS( (I+3)*(-PI/512.0D+00) )
  25.       A4=DCOS( (I+4)*(-PI/512.0D+00) )
  26.       WRITE (3,2) LFEED,I,I1,I2,I3,I4
  27. 1     WRITE (3,3) LFEED,A,A1,A2,A3,A4
  28. 2     FORMAT (A1, '     X ', 5( ',W(1,', I3, ')' ), '/')
  29. 3     FORMAT (A1, '     X ', 4(F10.7, ','), F10.7, '/')
  30.       DO 6 I=1,512,5
  31.       I1=I+1
  32.       I2=I+2
  33.       I3=I+3
  34.       I4=I+4
  35.       A =DSIN( (I  )*(-PI/512.0D+00) )
  36.       A1=DSIN( (I+1)*(-PI/512.0D+00) )
  37.       A2=DSIN( (I+2)*(-PI/512.0D+00) )
  38.       A3=DSIN( (I+3)*(-PI/512.0D+00) )
  39.       A4=DSIN( (I+4)*(-PI/512.0D+00) )
  40.       WRITE (3,7) LFEED,I,I1,I2,I3,I4
  41. 6     WRITE (3,8) LFEED,A,A1,A2,A3,A4
  42. 7     FORMAT (A1, '     X ', 5( ',W(2,', I3, ')' ), '/')
  43. 8     FORMAT (A1, '     X ', 4(F10.7, ','), F10.7, '/')
  44.       ENDFILE 3
  45.       END
  46.