[BACK]Return to fftws CVS log [TXT] [DIR] Up to [CSRG BSD Unix] / 43BSD / contrib / apl / lib

Annotation of 43BSD/contrib/apl/lib/fftws, revision 1.1

1.1     ! root        1: n��*.KB�       H fftx{ fft z ;n;n2;mk;bitrev;i2;arg;cs;t;i1
        !             2:  n{1YRz
        !             3:  x{z
        !             4:  n2{n%2
        !             5:  bitrev { 2B O-((2O*n)R2)NIn
        !             6: loop: mk{(n2R1),n2R0
        !             7:  i2 { n2+i1 { (nRmk)/In
        !             8:  arg { bitrev[Di1%n2]%n
        !             9:  cs { O\ 2 1 J.O 2XO arg
        !            10:  t { x[i2; 0 1 1 0]Xcs,cs
        !            11:  t { (t[;,0]+t[;,1]),t[;,2]-t[;,3]
        !            12:  x[i2;] { x[i1;]-t
        !            13:  x[i1;] { x[i1;]+t
        !            14:  }(1$n2 { n2%2)/loop
        !            15:  x { x[bitrev;]
        !            16: &aX&�@&n&&&&�A&n2&&@&mk$&&�@&i2$&&�@@A&i1$&&A&arg$&&�?&bitrev4&&A�@@A&cs8&�@���˙���@&t8&A&i&&D    rampr { ramp steps
        !            17: r { (O|Isteps +1),1UIsteps
        !            18: &rmp�&�A@AA�@�@A@A  timetime prog; FtiFm
        !            19: FtiFm { NB21
        !            20: E prog
        !            21: L { ''
        !            22: L { 'execution time is',(,B (.01XS((BN21)-FtiFm)%.6)), ' secs.'
        !            23: ''
        !            24:               magphaseresult { magphase y ;atans;phase;real0
        !            25: real0 { y[;0]=0
        !            26: atans { `3Oy[;1]%real0 +y[;0]
        !            27: phase { (~real0 )X((Oy[;0]<0)XXy[;1])+ atans
        !            28: phase { phase +0.5XO real0 XXy[;1]
        !            29: result { O\(O|Ry)R((+/yXy)*0.5),phase
        !            30: 
unix.superglobalmegacorp.com

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.