RADB3M.inc

1 C***********************************************************************
2       SUBROUTINE RADB3M (INC,LOT,IDO,L1,CC,CH,WA1,WA2)
3       DIMENSION  CC(INC,IDO,3,L1)  ,CH(INC,IDO,L1,3)  ,WA1(*)  ,WA2(*)
4       DATA TAUR,TAUI /-.5,.8660254037844390/
5 C
6       DO 101 K=1,L1
7 *VOPTION NOFVAL
8         DO 101 L=1,LOT
9           TR2         = CC(L,IDO,2,K)+CC(L,IDO,2,K)
10           CR2         = CC(L,  1,1,K)+TAUR*TR2
11           CI3         = TAUI*(CC(L,1,3,K)+CC(L,1,3,K))
12           CH(L,1,K,1) = CC(L,  1,1,K)+TR2
13           CH(L,1,K,2) = CR2-CI3
14           CH(L,1,K,3) = CR2+CI3
15   101   CONTINUE
16       IF (IDO .GT. 1) THEN
17           IDP2 = IDO+2
18         DO 103 K=1,L1
19 *VOPTION NOFVAL
20         DO 103 I=3,IDO,2
21             IC = IDP2-I
22 *VOPTION NOFVAL
23           DO 102 L=1,LOT
24             TR2           = CC(L,I-1,3,K)+CC(L,IC-1,2,K)
25             TI2           = CC(L,I  ,3,K)-CC(L,IC  ,2,K)
26             CR2           = CC(L,I-1,1,K)+TAUR*TR2
27             CI2           = CC(L,I  ,1,K)+TAUR*TI2
28             CR3           = TAUI*(CC(L,I-1,3,K)-CC(L,IC-1,2,K))
29             CI3           = TAUI*(CC(L,I  ,3,K)+CC(L,IC  ,2,K))
30             CH(L,I-1,K,1) = CC(L,I-1,1,K)+TR2
31             CH(L,I  ,K,1) = CC(L,I  ,1,K)+TI2
32             DR2           = CR2-CI3
33             DI2           = CI2+CR3
34             DR3           = CR2+CI3
35             DI3           = CI2-CR3
36             CH(L,I-1,K,2) = WA1(I-2)*DR2-WA1(I-1)*DI2
37             CH(L,I  ,K,2) = WA1(I-2)*DI2+WA1(I-1)*DR2
38             CH(L,I-1,K,3) = WA2(I-2)*DR3-WA2(I-1)*DI3
39             CH(L,I  ,K,3) = WA2(I-2)*DI3+WA2(I-1)*DR3
40   102     CONTINUE
41   103   CONTINUE
42       END IF
43       RETURN
44       END SUBROUTINE RADB3M
45