blob: dfebf0ccad56402b22ea5791cc56edff44305924 (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
|
c
c-----------------------------------------------------------------------
c subroutine: rad4sb
c used by subroutine radix4. never directly accessed by user.
c-----------------------------------------------------------------------
c
subroutine rad4sb(ntype)
c
c input: ntype = type of butterfly invoked
c output: parameters used by subroutine radix4
c
dimension ix(2996)
common /xx/ix
common ntypl,kkp,index,ixc
if(ntype.eq.ntypl) go to 7
ix(ixc)=0
ix(ixc+1)=ntype
ixc=ixc+2
if(ntype.ne.4) go to 4
indexp=(index-1)*9
ix(ixc)=kkp+1
ix(ixc+1)=indexp+1
ixc=ixc+2
go to 6
4 ix(ixc)=kkp+1
ixc=ixc+1
6 ntypl=ntype
return
7 if(ntype.ne.4) go to 8
indexp=(index-1)*9
ix(ixc)=kkp+1
ix(ixc+1)=indexp+1
ixc=ixc+2
return
8 ix(ixc)=kkp+1
ixc=ixc+1
return
end
|
