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