RFTF3M.inc

1 C***********************************************************************
2       SUBROUTINE RFTF3M (N,INC,LOT, C, WA,IFAC, CH)
3       DIMENSION       C(INC,N)   ,CH(INC,N)  ,WA(N)    ,IFAC(*)
4         CALL RFTF9M (N,INC,LOT, NA,C, WA, IFAC, CH)
5 C
6            CF = 1.0/FLOAT(N)
7            N4=(N/4)*4
8       IF (NA .EQ. 1) THEN
9         IF (N4 .GE. 4) THEN
10           DO 117 K=1,N4,4
11           DO 117 L=1,LOT
12             C(L,K  ) = CF*C(L,K  )
13             C(L,K+1) = CF*C(L,K+1)
14             C(L,K+2) = CF*C(L,K+2)
15             C(L,K+3) = CF*C(L,K+3)
16   117     CONTINUE
17         ENDIF
18         IF (N4 .NE. N) THEN
19           DO 118 K=N4+1,N
20           DO 118 L=1,LOT
21             C(L,K) = CF*C(L,K)
22   118     CONTINUE
23         ENDIF
24       ELSE
25         IF (N4 .GE. 4) THEN
26           DO 119 K=1,N4,4
27           DO 119 L=1,LOT
28             C(L,K  ) = CF*CH(L,K  )
29             C(L,K+1) = CF*CH(L,K+1)
30             C(L,K+2) = CF*CH(L,K+2)
31             C(L,K+3) = CF*CH(L,K+3)
32   119     CONTINUE
33         ENDIF
34         IF (N4 .NE. N) THEN
35           DO 120 K=N4+1,N
36           DO 120 L=1,LOT
37             C(L,K) = CF*CH(L,K)
38   120     CONTINUE
39         ENDIF
40       ENDIF
41       RETURN
42       END SUBROUTINE RFTF3M
43