C TEST OF 1024 COMPLEX POINT RADIX 4 FFT C LINKING:RX10B=RX10B/F,RDX405,RX4BUT DIMENSION C(2048) DIMENSION A(2048) COMMON /XX/A PI=3.141592 RMAX=1024. N=1024 XN=N DO 1 I=1,N XI=I-1 C(2*I)=0 1 C(2*I-1)=RMAX*SIN(2.*PI*XI/XN) 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