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 ;
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.
Subscribe to:
Post Comments (Atom)
203 comments:
«Oldest ‹Older 201 – 203 of 203ツイッターで出 合 いを求めるのです。気の合う異性と交流して楽しいひと時をお過ごしください。登録無料で使えるので気軽さは100点満点!
夏休みももう終わりに近づきこのまま一人は寂しいのちがう?そこでモテる度チェッカーを使い自分がなぜモテないか診断してもらいましょう。10~20代の女性アドバイザーが多数在籍してるので、アドバイスを参考に夏終盤を楽しもう
モバゲーで専用でぁいコミュニティOPEN!でぁいを求める男性と女性に素敵な時間をご提供します。
Post a Comment