C TEST OF 1024 POINT COMPLEX FLOATING POINT FFT C LINKING: RX10=RX10/F,RDX405,RX4BUT DIMENSION C(2048) DIMENSION A(2048) COMMON /XX/A PI=3.141592 RMAX=2.**13-1. N=1024 XN=N C WAVEFORM USED IS FOURIER SERIES KERNEL, A "PERIODIC" C SINX/X. DFT IS A CONSTANT FOR THIS SIGNAL. C FOR DETAILS, SEE PAPAOULIS, "THE FOURIER INTEGRAL C AND IT'S APPLICATIONS", PP.44-45, MCGRAW-HILL 1962. XN2=XN/2.+.5 XP1=XN+1. XX=RMAX/XP1 C(1)=XX*XP1 C(2)=0 DO 1 I=2,N XI=I-1 C(2*I)=0 1 C(2*I-1)=XX*SIN(2.*XN2*PI*XI/XN)/SIN(PI*XI/XN) C WRITE(7,25) 25 FORMAT(1H ,'INPUT N FOR N**2 ITERATIONS, FORMAT I4') READ(7,26) NS 26 FORMAT(I4) XS=NS XS=XS*XS WRITE(7,27) XS 27 FORMAT(1H ,'NUMBER OF TIMES=', 2X,E14.6) C DO 10 JJ=1,NS DO 10 KK=1,NS DO 2 J=1,N A(2*J-1)=C(2*J-1) 2 A(2*J)=C(2*J) CALL RDX405 10 CONTINUE DO 3 J=1,N WRITE(7,5) A(2*J-1),A(2*J) 5 FORMAT(1H ,2(F14.6,5X)) 3 CONTINUE STOP END