[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

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.