計算

EAM ポテンシャルの分子動力学計算にCG法使うため二階微分を計算、コードに入れる。
http://d.hatena.ne.jp/ita/19000102
 \nabla_i \nabla_i E という対角ブロックだけでとりあえず近似。本当は粒子数*3 次元の行列の逆行列計算が必要だけど面倒なんで。

しまった明日はサンクスギビングで休日、翌日からボストンのMRSに行くため一週間つぶれる。何もできん。なんだらかんだらたまってるのに。