RFTB1M.inc

1 C***********************************************************************
2       SUBROUTINE RFTB1M (N,INC,LOT, C, WA,IFAC, CH)
3       DIMENSION       C(INC,N)   ,CH(INC,N)  ,WA(N)    ,IFAC(*)
4         CALL RFTB9M (N,INC,LOT, NA,C, WA, IFAC, CH)
5       IF (NA .EQ. 1) RETURN
6            N4 = (N/4)*4
7       IF (N4 .GE. 4) THEN
8          DO 117 K=1,N4,4
9          DO 117 L=1,LOT
10            C(L,K  ) = CH(L,K  )
11            C(L,K+1) = CH(L,K+1)
12            C(L,K+2) = CH(L,K+2)
13            C(L,K+3) = CH(L,K+3)
14   117    CONTINUE
15       ENDIF
16       IF (N4 .NE. N) THEN
17          DO 118 K=N4+1,N
18          DO 118 L=1,LOT
19            C(L,K) = CH(L,K)
20   118    CONTINUE
21       ENDIF
22       RETURN
23       END SUBROUTINE RFTB1M
24