SUBROUTINE VRADB3 (MP,IDO,L1,CC,CH,MDIMC,WA1,WA2)
C
C VRFFTPK, VERSION 1, AUGUST 1985
C
DIMENSION CC(MDIMC,IDO,3,L1) ,CH(MDIMC,IDO,L1,3),
1 WA1(IDO) ,WA2(IDO)
ARG=2.*PIMACH(1.0)/3.
TAUR=COS(ARG)
TAUI=SIN(ARG)
DO 101 K=1,L1
DO 1001 M=1,MP
CH(M,1,K,1) = CC(M,1,1,K)+2.*CC(M,IDO,2,K)
CH(M,1,K,2) = CC(M,1,1,K)+(2.*TAUR)*CC(M,IDO,2,K)
1 -(2.*TAUI)*CC(M,1,3,K)
CH(M,1,K,3) = CC(M,1,1,K)+(2.*TAUR)*CC(M,IDO,2,K)
1 +2.*TAUI*CC(M,1,3,K)
1001 CONTINUE
101 CONTINUE
IF (IDO .EQ. 1) RETURN
IDP2 = IDO+2
DO 103 K=1,L1
DO 102 I=3,IDO,2
IC = IDP2-I
DO 1002 M=1,MP
CH(M,I-1,K,1) = CC(M,I-1,1,K)+(CC(M,I-1,3,K)+CC(M,IC-1,2,K))
CH(M,I,K,1) = CC(M,I,1,K)+(CC(M,I,3,K)-CC(M,IC,2,K))
CH(M,I-1,K,2) = WA1(I-2)*
1 ((CC(M,I-1,1,K)+TAUR*(CC(M,I-1,3,K)+CC(M,IC-1,2,K)))-
* (TAUI*(CC(M,I,3,K)+CC(M,IC,2,K))))
2 -WA1(I-1)*
3 ((CC(M,I,1,K)+TAUR*(CC(M,I,3,K)-CC(M,IC,2,K)))+
* (TAUI*(CC(M,I-1,3,K)-CC(M,IC-1,2,K))))
CH(M,I,K,2) = WA1(I-2)*
4 ((CC(M,I,1,K)+TAUR*(CC(M,I,3,K)-CC(M,IC,2,K)))+
8 (TAUI*(CC(M,I-1,3,K)-CC(M,IC-1,2,K))))
5 +WA1(I-1)*
6 ((CC(M,I-1,1,K)+TAUR*(CC(M,I-1,3,K)+CC(M,IC-1,2,K)))-
8 (TAUI*(CC(M,I,3,K)+CC(M,IC,2,K))))
CH(M,I-1,K,3) = WA2(I-2)*
7 ((CC(M,I-1,1,K)+TAUR*(CC(M,I-1,3,K)+CC(M,IC-1,2,K)))+
8 (TAUI*(CC(M,I,3,K)+CC(M,IC,2,K))))
8 -WA2(I-1)*
9 ((CC(M,I,1,K)+TAUR*(CC(M,I,3,K)-CC(M,IC,2,K)))-
8 (TAUI*(CC(M,I-1,3,K)-CC(M,IC-1,2,K))))
CH(M,I,K,3) = WA2(I-2)*
1 ((CC(M,I,1,K)+TAUR*(CC(M,I,3,K)-CC(M,IC,2,K)))-
8 (TAUI*(CC(M,I-1,3,K)-CC(M,IC-1,2,K))))
2 +WA2(I-1)*
3 ((CC(M,I-1,1,K)+TAUR*(CC(M,I-1,3,K)+CC(M,IC-1,2,K)))+
8 (TAUI*(CC(M,I,3,K)+CC(M,IC,2,K))))
1002 CONTINUE
102 CONTINUE
103 CONTINUE
RETURN
END
.
