1 C FFTC.FT FFT-V4A 4/22/76
3 SUBROUTINE FFTC (A,N,M,SGNEX)
5 C-------FFT VON N=2**M KOMPLEXEN DATEN
6 C A = IN- UND OUTPUTARRAY, LAENGE N=2**M, KOMPLEX
7 C-------SGNEXP = VORZEICHEN DES EXPONENTEN DER E-FUNKTIONEN
10 DIMENSION U(2),W(2),T(2),V(2)
42 T(1)=A(1,IP)*U(1)-A(2,IP)*U(2)
43 T(2)=A(1,IP)*U(2)+A(2,IP)*U(1)
49 V(1)=U(1)+U(1)*W(1)-U(2)*W(2)
50 V(2)=U(2)+U(1)*W(2)+U(2)*W(1)
53 COR=.5-.5*(U(1)*U(1)+U(2)*U(2))