软件进展（二零零八年三月）

三月一日二日 重新写了算力的部分，以求效率上获得本质改进。但是发现力算的不对，找了半天错误，发现其实误差是在允许范围之内的，崩溃。。。
三月八日 对于 \psi 和 \phi 中无穷求和做了精度上的优化，发现对于K为32及以上时没有什么效果。把 c2c 的 Fourier transform 修改为 r2c and c2r 的 Fourier transform。貌似 Fourier transform 部分应该有50％的时间节省，但是在笔记本上无法测试。同时空间上50％的节省是肯定的。
三月九日 完成了 B-Spline 的程序，精度提高一个数量级，效率方面有待测试。把原先的 ElectrostaticInteraction.cpp 分成了几个.cpp。作为副产品的是递归计算 Cardinal B-Splines 的程序。
三月二十四日 大概推导了一下实空间部分的误差，但是实际证明这个估计给的太过夸张，大概大了4个数量级。
三月二十五－二十八日 找到了对 Ewald sum 和 PME 误差了理论分析的文章。发现能量随着 \beta 的变化有微小的变化，对于 -80 左右的能量，变化可以达到 0.4 左右。终于发现了实空间能量部分有一个 bug，更正后满足了误差的理论估计。理论估计的妙处在假设 cut off 之外电荷均匀分布，这个假设是合理的，因为对于液体来说 rdf 是收敛到 1 的。
三月三十一日 掌握了 PME 误差分析的技术，准备做 B-spline 的误差分析，但是发现不太知道 B-spline 的逼近结果。。。