RADF2M.inc

1 C***********************************************************************
2       SUBROUTINE RADF2M (INC,LOT,IDO,L1,CC,CH,WA1)
3       DIMENSION       CH(INC,IDO,2,L1)  ,CC(INC,IDO,L1,2)  ,WA1(*)
4       DO 101 K=1,L1
5         DO 101 L=1,LOT
6           CH(L,  1,1,K) = CC(L,1,K,1)+CC(L,1,K,2)
7           CH(L,IDO,2,K) = CC(L,1,K,1)-CC(L,1,K,2)
8   101   CONTINUE
9       IF (MOD(IDO,2) .EQ. 0) THEN
10         DO 102 K=1,L1
11           DO 102 L=1,LOT
12             CH(L,  1,2,K) = -CC(L,IDO,K,2)
13             CH(L,IDO,1,K) =  CC(L,IDO,K,1)
14   102     CONTINUE
15       END IF
16       IF (IDO .GT. 2) THEN
17           IDP2 = IDO+2
18         DO 104 K=1,L1
19 *VOPTION NOFVAL
20         DO 104 I=3,IDO,2
21             IC = IDP2-I
22 *VOPTION NOFVAL
23           DO 103 L=1,LOT
24             TR2 = WA1(I-2)*CC(L,I-1,K,2)+WA1(I-1)*CC(L,I  ,K,2)
25             TI2 = WA1(I-2)*CC(L,I  ,K,2)-WA1(I-1)*CC(L,I-1,K,2)
26             CH(L,I -1,1,K) = CC(L,I-1,K,1)+TR2
27             CH(L,I   ,1,K) = TI2          +CC(L,I,K,1)
28             CH(L,IC-1,2,K) = CC(L,I-1,K,1)-TR2
29             CH(L,IC  ,2,K) = TI2          -CC(L,I,K,1)
30   103     CONTINUE
31   104   CONTINUE
32       END IF
33       RETURN
34       END SUBROUTINE RADF2M
35