next up previous
Next: 4 波動方程式の差分近似 Up: メモ Previous: 2 MitsuiWorld の改良

3 Java で FFT (1)

どうするのが一番良いかイマイチ分からないが、 とりあえず Ooura's Mathematical Software Packages にある 「汎用 FFT (高速 フーリエ/コサイン/サイン 変換) パッケージ」 を使ってみる。

これはデータ数が $ 2$ の羃でないとダメという制限はあるが、 性能の高さには定評がある。

oyabun% tar xzf fft.tgz
oyabun% cd fft
oyabun% ls
oyabun% cd sample1
oyabun% make
mathpc% ./test4g
data length n=? (must be 2^m)
1024
cdft err= 4.24682e-16 
rdft err= 4.24682e-16 
ddct err= 4.95724e-16 
ddst err= 4.98625e-16 
dfct err= 3.86626e-16 
dfst err= 5.69098e-16 
oyabun% cd ../sample2 
oyabun% make
oyabun% ./pi_fft4g
oyabun% less pi.dat

例えば test4g.c を Java プログラムに書き換えると良い。

I君はその作業をスタートしたようだが、上手く行ったかな?

FFTPACK を Java にするのは無理か? あるじゃん! jfftpack.tgz -- でも上手くコンパイル出来ない。なぜだろう。


next up previous
Next: 4 波動方程式の差分近似 Up: メモ Previous: 2 MitsuiWorld の改良
Masashi Katsurada
平成20年2月11日