市來健吾の日記

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

周期系 FMM


  • 最単純 FMM 実装、その5

    • bug 一つ、 周期では level 0 まで計算しないとダメだった

    • 非周期と周期の code を排他的にしていたが、 周期を非周期の上に乗っけることにする

    • 全体の test code を書く

    • 間違い発見

      • real 部も格子和をとっていた

      • recip 部に self の寄与を考えていなかった

      • これで order は合う程度まで来た