cfft1b cfft1f cfft1i cfftmb cfftmf cfftmi

cfftmf() void cfftmf ( int  lot, int  jump, int  n, int  inc, doublecomplex  c[], int  lc, double  wsave[], int  lwsave, double  work[], int  lwork, int *  info  ) 1次元複素フーリエ変換 (複数データ列) 目的 cfftmfは複素配列中の複数の周期数列の1次元フーリエ変換を計算する. この変換はフーリエ変換あるいはフーリエ解析と呼ばれ, 数列を物理空間からスペクトル空間に変換する. c[l*jump+k] = (1/n)Σc[l*jump+j]exp(-2πijk/n) (Σは j = 0 〜 n-1) (l = 0 〜 lot-1, k = 0 〜 n-1) (iは虚数単位) この変換は正規化されており, cfftmbに続くcfftmfの呼び出し(あるいはその逆)により, アルゴリズム上の制約, 丸め誤差などを除き, 元の配列を復元する. 引数 [in] lot 入力データ列の数. (lot >= 1) [in] jump c[]内でのデータ列の最初の要素と次のデータ列の最初の要素の間隔. (jump >= 1) [in] n 各入力データ列の長さ. (n >= 1) nが小さな素数の積で表されると効率が良い. [in] inc データ列の連続する要素の配列c[]内でのインデックスの間隔. (inc >= 1) [in,out] c[] 配列 c[lc] [in] 入力データ列. [out] フーリエ変換されたデータ列. [in] lc 配列 c[] の大きさ. (lc >= (lot - 1)*jump + inc*(n - 1) + 1) [in] wsave[] 配列 wsave[lwsave] 作業データ. 入力データ列の長さnごとに, cfftmfあるいはcfftmbを最初に呼び出す前にcfftmiにより初期化しておかなければならない. [in] lwsave 配列 wsave[] の大きさ. (lesave >= 2*n + ln(n)/ln(2) + 4) [out] work[] 配列 work[lwork] 作業領域. [in] lwork 配列 work[] の大きさ. (lwork >= 2*lot*n) [out] info = 0: 正常終了 = -1: 入力パラメータ lot の誤り (lot < 1) = -2: 入力パラメータ jump の誤り (jump < 1) = -3: 入力パラメータ n の誤り (n < 1) = -4: 入力パラメータ inc の誤り (inc < 1) = -6: 入力パラメータ lc の誤り (lcが小さい) = -8: 入力パラメータ lwsave の誤り (lwsaveが小さい) = -10: 入力パラメータ lwork の誤り (lworkが小さい) 出典 FFTPACK 5.1