DIMENSION XR(128),XI(128) DATA XR/128*0./,XI/128*0./ TPI=8*ATAN(1.) N=128 XN=N KK=6 XINC=TPI/XN C GENERATE THE TEST TIME SEQUENCE CONSISTING OF A COSINE FUNDEMENTAL C AND A 3RD HARMONIC COSINE, AND A DC TERM. 701 FORMAT('0THIS PROGRAM GENERATE THE TEST TIME SEQUENCE CONSISTING O 1F A COSINE'/' FUNDEMENTAL, A 3RD HARMONIC COSINE, AND A DC TERM.'//) TYPE 701 DO 1 I=1,N ARG=(I-1)*XINC 1 XR(I)=COS(ARG)+1.0 DO 2 I=1,N ARG=(I-1)*3*XINC 2 XR(I)=XR(I)-COS(ARG)/3. C TRANSFORM AND UNSCRAMBLE CALL FFTBA2(XR,XI,N,1,1) C PRINT OUT SOME OF THE VALUES. TYPE901,KK 901 FORMAT('0THE FIRST AND LAST ',I4,' VALUES OF THE TRANSFORMED ARRAY 1:') CALL YPRNT(XR,XI,N,KK) CALL FOLDFT(XR,XI,N,-1) TYPE 902,KK 902 FORMAT('0THE FIRST AND LAST ',I4,' VALUES OF THE TRANSFORMED ARRAY ' 1/' AFTER FOLDING USING SUBROUTINE FOLDFT WITH ITYPE = -1 :') CALL YPRNT(XR,XI,N,KK) CALL EXIT END SUBROUTINE YPRNT(XR,XI,N,KK) DIMENSION XR(N),XI(N) 888 FORMAT('+ ') 889 FORMAT('+ --------NOTE--------') TYPE 900 900 FORMAT(9X,'INDEX',5X,'VECTOR XR',5X,'VECTOR XI') DO 801 I=1,KK TYPE 901,I,XR(I),XI(I) IF(I.LT.3.OR.I.EQ.4)TYPE 889 IF(.NOT.LK)TYPE 888 801 CONTINUE KS=N-KK+1 DO 802 I=KS,N TYPE 901,I,XR(I),XI(I) IF(I.EQ.N.OR.I.EQ.N-2)TYPE 889 901 FORMAT(8X,I5,1X,1P2E14.4,$) TYPE 888 802 CONTINUE RETURN END