FFTWをFortranで使う

今仕事で書いてるプログラムの仕様…というかいわゆる数式を眺めていて「一見数値積分しろといわんばかりではあるけど、どう考えてもフーリエ変換だよね？」ということに気付いて慌ててFFTWを準備。

インストールはlennyなのでいいとして、考えるとFortranでFFTW使ったことないかもと。当然ドキュメント読めばいいのだけど、誰かF90辺りでサンプルコード書いてないかなーと思ったらあった。

心から感謝。