1 C FFTR.FT FFT-V4A 4/22/76
3 SUBROUTINE FFTR (A,N,DIR,SGNEX)
5 C-------ZUSATZPROGRAMM ZUR FFT FUER 2**N REELLE DATEN
6 C A = IN- UND OUTPUTARRAY, LAENGE N+1, KOMPLEX (A,B)
7 C A(1),...,A(N) ENTHALTEN DIE 2*N REELLEN DATEN
9 C CALL REALTR MIT SGNEXP=+1.
10 C RUECKWAERTS CALL REALTR
11 C CALL FFT MIT SGNEXP=-1.
12 C-------(DAZU IN EINER RICHTUNG MIT 1/N NORMIEREN)
17 SD=3.14159265/FLOAT(N)
22 IF (DIR.LT.0.) GOTO 10
41 C THE FOLLOWING THREE STATMENTS COMPENSATE FOR TRUNCATION
42 C ERROR. IF ROUNDED ARITHMETIC IS USED, SUBSTITUTE
44 CN=.5-.5*(AA*AA+SN*SN)