下面做点啥

1，周三和老板聊过了，因此要开始新的工作。先把他们的工作弄明白！！很重要。要有信心并且脑子灵活

2，筹划写写短程相互作用的论文了。

3，给老张打个电话。

4，把 analytical diff 的程序改进下：在估计平方误差的时候请忽略高频震荡，因为对结果影响不本质。

5，想想从解析表达式能得到有什么有趣的性质，对什么样的体系进行模拟，能说明什么问题。。

把误差估计的文档写了

14页公式，写的手疼。。

终于完成 analytical differetiation 的验证。。

程序一千三百多行。。。。

对测试的结果观察有问题。。

算的好慢

Analytical differentiation...

通过平均误差测试

在重新推方差误差公式。。

Analytical differentiation...

验证程序写成了一坨 shi 。。。

周日写程序就是不中。。。

SPME ik

测试通过 SPME ik，非均匀体系。

下面去做 analytical 的误差估计吧。

bug 坑死爹。。。

上午实现 SPME 的 ik 部分，把误差力的平均写好后脑抽认为应该先去验证下 Ewald sum 的这部分结果。于是一个下午都在 debug 中。。。。怎么调怎么发现我的估计是所有粒子都是正电荷的结果，但是显然不对啊，因为一半正一半负的。。。最后发现是 Ewald sum 的程序的电荷赋值有个 bug。。坑爹啊，一个下午啊。。。这从侧面说明误差估计的有多准

一个月过去了干成了哪些

0，done

1，写出了 Ewald sum 和 SPME ik 的误差估计，计算量的阶为 N log N。昨天数值验证了 Ewald sum。现在在写 ik 的验证程序。analytical 完全类似，但是我想先验证了 ik 再来做 analytical。

2，done

3，done

看来在多一吃帝国清心寡欲效率还挺高。。。