RFTB9M.inc
References to this file elsewhere.
1 C***********************************************************************
2 SUBROUTINE RFTB9M (N,INC,LOT, NA,C, WA,IFAC, CH)
3 DIMENSION C(INC,N) ,CH(INC,N) ,WA(N) ,IFAC(*)
4 NF = IFAC(2)
5 NA = 1
6 L1 = 1
7 IW = 1
8 DO 116 K1=1,NF
9 IP = IFAC(K1+2)
10 L2 = IP*L1
11 IDO = N/L2
12 IDL1 = IDO*L1
13 NA = 1-NA
14 IF (IP .EQ. 4) THEN
15 IX2 = IW +IDO
16 IX3 = IX2+IDO
17 IF (NA .EQ. 0) THEN
18 CALL RADB4M (INC,LOT,IDO,L1,C,CH,WA(IW),WA(IX2),WA(IX3))
19 ELSE
20 CALL RADB4M (INC,LOT,IDO,L1,CH,C,WA(IW),WA(IX2),WA(IX3))
21 ENDIF
22 ELSE IF (IP .EQ. 2) THEN
23 IF (NA .EQ. 0) THEN
24 CALL RADB2M (INC,LOT,IDO,L1,C,CH,WA(IW))
25 ELSE
26 CALL RADB2M (INC,LOT,IDO,L1,CH,C,WA(IW))
27 ENDIF
28 ELSE IF (IP .EQ. 3) THEN
29 IX2 = IW +IDO
30 IF (NA .EQ. 0) THEN
31 CALL RADB3M (INC,LOT,IDO,L1,C,CH,WA(IW),WA(IX2))
32 ELSE
33 CALL RADB3M (INC,LOT,IDO,L1,CH,C,WA(IW),WA(IX2))
34 ENDIF
35 ELSE IF (IP .EQ. 5) THEN
36 IX2 = IW +IDO
37 IX3 = IX2+IDO
38 IX4 = IX3+IDO
39 IF (NA .EQ. 0) THEN
40 CALL RADB5M (INC,LOT,IDO,L1,C,CH
41 * ,WA(IW),WA(IX2),WA(IX3),WA(IX4))
42 ELSE
43 CALL RADB5M (INC,LOT,IDO,L1,CH,C
44 * ,WA(IW),WA(IX2),WA(IX3),WA(IX4))
45 ENDIF
46 ELSE
47 112 IF (NA .EQ. 0) THEN
48 CALL RADBGM (INC,LOT,IDO,IP,L1,IDL1,C,C,C,CH,CH,WA(IW))
49 NA = 1
50 ELSE
51 CALL RADBGM (INC,LOT,IDO,IP,L1,IDL1,CH,CH,CH,C,C,WA(IW))
52 NA = 0
53 ENDIF
54 IF (IDO .EQ. 1) NA = 1-NA
55 ENDIF
56 L1 = L2
57 IW = IW+(IP-1)*IDO
58 116 CONTINUE
59 RETURN
60 END SUBROUTINE RFTB9M
61