diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
| commit | fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch) | |
| tree | bdda434976bc09c864f2e4fa6f16ba1952b1e555 /math/ieee/chap1/rad4sb.f | |
| download | iraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz | |
Initial commit
Diffstat (limited to 'math/ieee/chap1/rad4sb.f')
| -rw-r--r-- | math/ieee/chap1/rad4sb.f | 38 |
1 files changed, 38 insertions, 0 deletions
diff --git a/math/ieee/chap1/rad4sb.f b/math/ieee/chap1/rad4sb.f new file mode 100644 index 00000000..dfebf0cc --- /dev/null +++ b/math/ieee/chap1/rad4sb.f @@ -0,0 +1,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 |