■
テクノロジー : 日経電子版
今日も足を取られっぱなし。コードが動かん。いままでのコードをちょこっといじっただけのコードだと killed の表示が出て動かん。しかたないので動くところから、じわじわとスクラッチして書き足す。複素フーリェのラッパーが変。実FFTで比較するとライブラリのラッパーが10数倍早いので、ラッパーのデバッグをする方が、最終的に結果が速く出ると考えて進めているのだが。
実2変数バージョンの複素フーリェ変換を使ったら Numerical Recipe 版と同じウェーブレット係数を出した。(精度の問題ではなくて、違う答えを出していたのだ。)これでいこう。それにしても実フーリェも複素フーリェも動作が変だ。同じようにマニュアルにしたがってコーディングしたのに。
夜中になってバグ出し完了。新規作成の ∇×(U×B) と以前バグ出ししたコードを2回呼び出す (B・∇)U-(U・∇)B が同じ結果を出したから多分だいじょうぶ。作成したいコードとバグ出しクロスチェック用コードの双方にバグが見つかった。多分、両方ともに全滅というバグ混入はもう無いだろう。
[典侍のいいまつがい]「おたませ!(お待たせ)」「みきしばい(紙芝居)」