RFFTFM.inc

1 C***********************************************************************
2       SUBROUTINE RFFTFM (N,INC,JUMP,LOT, R, WA,IFAC, WSAVE)
3       DIMENSION       R(*)       ,WSAVE(*)   ,WA(N)    ,IFAC(*)
4 cch?     REAL*8 CTIME1,CTIME2
5 cch?     CALL XCLOCK(CTIME1,5)
6 C
7       IF (N .EQ. 1) RETURN
8 C
9       IF (JUMP .NE. 1) THEN
10         CALL RFTF2M (N,INC,JUMP,LOT, R, WA,IFAC, WSAVE)
11       ELSE
12         CALL RFTF3M (N,INC,     LOT, R, WA,IFAC, WSAVE)
13       ENDIF
14 cch?     CALL XCLOCK(CTIME2,5)
15 cch?     CTIME1=CTIME2-CTIME1
16 cch?     CTIME2=1000.0*CTIME1/FLOAT(N*LOT)
17 cch?     WRITE(6,500) INC,JUMP,N,LOT, CTIME1,CTIME2
18 cch? 500 FORMAT(1H0,' RFFTBM: INC,JUMP,N,LOT =',4I6,'  TOTAL CPU =',F13.6
19 cch?    *      ,' (SEC) ;  UNIT CPU =',F13.6,' (MSEC)')
20       RETURN
21       END SUBROUTINE RFFTFM
22