1 SUBROUTINE RADF3 (IDO,L1,CC,CH,WA1,WA2) 2 IMPLICIT DOUBLE PRECISION (A-H,O-Z) 3 DIMENSION CH(IDO,3,L1) ,CC(IDO,L1,3) , 4 1 WA1(1) ,WA2(1) 5C *** TAUI IS -SQRT(3)/2 *** 6 DATA TAUR,TAUI /-0.5D0,0.86602540378443864676D0/ 7 DO 101 K=1,L1 8 CR2 = CC(1,K,2)+CC(1,K,3) 9 CH(1,1,K) = CC(1,K,1)+CR2 10 CH(1,3,K) = TAUI*(CC(1,K,3)-CC(1,K,2)) 11 CH(IDO,2,K) = CC(1,K,1)+TAUR*CR2 12 101 CONTINUE 13 IF (IDO .EQ. 1) RETURN 14 IDP2 = IDO+2 15 DO 103 K=1,L1 16 DO 102 I=3,IDO,2 17 IC = IDP2-I 18 DR2 = WA1(I-2)*CC(I-1,K,2)+WA1(I-1)*CC(I,K,2) 19 DI2 = WA1(I-2)*CC(I,K,2)-WA1(I-1)*CC(I-1,K,2) 20 DR3 = WA2(I-2)*CC(I-1,K,3)+WA2(I-1)*CC(I,K,3) 21 DI3 = WA2(I-2)*CC(I,K,3)-WA2(I-1)*CC(I-1,K,3) 22 CR2 = DR2+DR3 23 CI2 = DI2+DI3 24 CH(I-1,1,K) = CC(I-1,K,1)+CR2 25 CH(I,1,K) = CC(I,K,1)+CI2 26 TR2 = CC(I-1,K,1)+TAUR*CR2 27 TI2 = CC(I,K,1)+TAUR*CI2 28 TR3 = TAUI*(DI2-DI3) 29 TI3 = TAUI*(DR3-DR2) 30 CH(I-1,3,K) = TR2+TR3 31 CH(IC-1,2,K) = TR2-TR3 32 CH(I,3,K) = TI2+TI3 33 CH(IC,2,K) = TI3-TI2 34 102 CONTINUE 35 103 CONTINUE 36 RETURN 37 END 38