Sunday, January 7, 2007

FFT

Wrote FFT in Factor. Developing is very nice. Just the final code looks a bit ugly because of all the stack operations. I guess that's the price you have to pay for a stack based programming language.


IN: fft

USING: arrays sequences math kernel ;

: omega ( n -- n )
recip -2 pi i * * * exp ;

: n^v ( n v -- w )
[ ^ ] map-with ;

: v^n ( v n -- w )
swap n^v ;

: even ( seq -- seq )
2 group 0 <column> ;

: odd ( seq -- seq )
2 group 1 <column> ;

: two ( seq -- seq )
1/2 v*n dup append ;

: twiddle ( seq -- seq )
dup length dup omega v^n v* ;

: fft ( seq -- seq )
dup length 1 =
[ dup odd fft two twiddle swap even fft two v+ ] unless ;

PROVIDE: demos/fft ;

203 comments:

«Oldest   ‹Older   201 – 203 of 203
ツイッター said...

ツイッターで出 合 いを求めるのです。気の合う異性と交流して楽しいひと時をお過ごしください。登録無料で使えるので気軽さは100点満点!

モテる度チェッカー said...

夏休みももう終わりに近づきこのまま一人は寂しいのちがう?そこでモテる度チェッカーを使い自分がなぜモテないか診断してもらいましょう。10~20代の女性アドバイザーが多数在籍してるので、アドバイスを参考に夏終盤を楽しもう

モバゲー said...

モバゲーで専用でぁいコミュニティOPEN!でぁいを求める男性と女性に素敵な時間をご提供します。

«Oldest ‹Older   201 – 203 of 203   Newer› Newest»