市來健吾の日記

プログラマ、(元)物理屋(ナノテク、流体)

周期系 FMM


  • 新 adap-FMM の debug

    • n2-2-2 に集中

      • 部分毎に押えて行く

      • NaN の原因は発見 (vd_tmp の初期化)

      • 問題は real にある

      • さらに絞る、問題は cp_tns_2_tns_ にある

      • さらに絞る、発見

      • オチは、cp_2_tns で、初期化が入っていた…

    • n3-3-3 に移行、別の bug

      • non-FMM と違う

      • 色々つついた結果、 calc_down_nnn... の children の扱いに問題があった (new だけでなく)

  • やっと debug 終了、登録

  • 展開精度の検証に移る!

    • code を書く

      • とりあえず、書き上げる

      • 結果は、すばらしくはないが、こんなものかもしれない、という程度のもの

      • 球面上の一様分布は、gsl の routine を使う

        • なんでも、Knuth の本に解説してある scheme らしい