周期系 FMM

• 
  library から sparse matrix の教科書を借りて来る
  • 
    saad (GMRES の開発者) の本もあったので借りておく
  • 
    Yousef Saad, "Iterative Methods for Sparse Linearsystems" (1996) ISBN: 053494776-x (amazon)
• 
  FFT による recip 部の高速化、その２
  • 
    補間に GSL を使う
    • 
      今までは悪しき Numerical Recipe を使っていた
    • 
      exp-periodic-fft.h: 1.2 2002/06/24 19:56:01 ichiki Exp
    • 
      exp-periodic-fft.c: 1.2 2002/06/24 19:55:33 ichiki Exp
  • 
    比較
    

    scheme	time [m sec]	comment
    GSL(linear)	229.046	精度悪い
    GSL(akima_periodic)	675.900	精度悪い
    NR(polint)	1067.065
    GSL(cspline)	1147.931
    GSL(cspline_periodic)	1500.142	精度良い

• 
  xi 設定の tuning
  • 
    全て 0 次関数の実測で限界を決める
    • 
      exp-periodic.c: 3.9 2002/06/24 21:51:15 ichiki Exp
  • 
    実際の lattice summation の打切り誤差の挙動
    • 
      parameter が打切り精度になっていない… こういうもの、とする