| |
| | [No title] (Site not responding. Last check: 2007-11-06) |
 | | 1) return c call srftb1 (n,r,wsave,wsave(n+1),wsave(2*n+1)) c return end subroutine srfftf (n,r,wsave) real r(1), wsave(1) c if (n.eq. |
 | | 1) return c call srfti1 (n,wsave(n+1),wsave(2*n+1)) c return end subroutine srftb1 (n,c,ch,wa,ifac) real c(1), ch(1), wa(1) integer ifac(*) c nf = ifac(2) na = 0 l1 = 1 iw = 1 do 116 k1=1,nf ip = ifac(k1+2) l2 = ip*l1 ido = n/l2 idl1 = ido*l1 if (ip.ne. |
 | | 1) go to 104 ifac(1) = n ifac(2) = nf c argh = tpi/float(n) is = 0 nfm1 = nf-1 l1 = 1 if (nfm1.eq. |
| www.astro.ku.dk /~aake/ComputerPhysics/13_StarFormation/fftpack.f (1318 words) |
|