5.7.2 大浦FFT

Cで書かれた 汎用 FFT (高速 フーリエ/コサイン/サイン 変換) パッケージ by 大浦拓哉 は、比較的簡単に Java に書き直せる。

コンパイル＆テスト

oyabun% tar xzf fft.tgz
oyabun% cd fft
oyabun% cd sample1
oyabun% make
oyabun% ./test4g
data length n=? (must be 2^m)
1024
cdft err= 5.55112e-16 
rdft err= 4.44089e-16 
ddct err= 7.77156e-16 
ddst err= 6.66134e-16 
dfct err= 4.44089e-16 
dfst err= 7.77156e-16 
 test4g_h, test8g, test8g_h, testsg, testsg_h についても同様
oyabun% cd ../sample2
oyabun% make
oyabun% ./pi_fft4g
  ... (略) ...

2007年度卒研で一木君は fft4g.c を Java に書き換えて使用した。