1 subroutine zfftf1 (n,c,ch,wa,ifac)
2 implicit double precision (a-h,o-z)
14 if (ip .ne. 4) go
to 103
17 if (na .ne. 0) go
to 101
18 call
zpassf4(idot,l1,c,ch,wa(iw),wa(ix2),wa(ix3))
20 101 call
zpassf4(idot,l1,ch,c,wa(iw),wa(ix2),wa(ix3))
23 103
if (ip .ne. 2) go
to 106
24 if (na .ne. 0) go
to 104
25 call
zpassf2(idot,l1,c,ch,wa(iw))
27 104 call
zpassf2(idot,l1,ch,c,wa(iw))
30 106
if (ip .ne. 3) go
to 109
32 if (na .ne. 0) go
to 107
33 call
zpassf3(idot,l1,c,ch,wa(iw),wa(ix2))
35 107 call
zpassf3(idot,l1,ch,c,wa(iw),wa(ix2))
38 109
if (ip .ne. 5) go
to 112
42 if (na .ne. 0) go
to 110
43 call
zpassf5(idot,l1,c,ch,wa(iw),wa(ix2),wa(ix3),wa(ix4))
45 110 call
zpassf5(idot,l1,ch,c,wa(iw),wa(ix2),wa(ix3),wa(ix4))
48 112
if (na .ne. 0) go
to 113
49 call
zpassf(nac,idot,ip,l1,idl1,c,c,c,ch,ch,wa(iw))
51 113 call
zpassf(nac,idot,ip,l1,idl1,ch,ch,ch,c,c,wa(iw))
52 114
if (nac .ne. 0) na = 1-na
subroutine zpassf5(ido, l1, cc, ch, wa1, wa2, wa3, wa4)
subroutine zpassf3(ido, l1, cc, ch, wa1, wa2)
subroutine zpassf4(ido, l1, cc, ch, wa1, wa2, wa3)
subroutine zfftf1(n, c, ch, wa, ifac)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
subroutine zpassf(nac, ido, ip, l1, idl1, cc, c1, c2, ch, ch2, wa)
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)
subroutine zpassf2(ido, l1, cc, ch, wa1)