多粒子系の Conjugate Gradient 緩和

CG法は関数f(x)の極値を求める方法。x -= f'(x)/f''(x) とすればいい。多次元なら二回微分行列の逆をかける。多粒子系だと逆行列とか非対角要素の計算面倒だなと思ってたけど、各粒子について他を全部固定して一つずつ順番に更新すれば問題なさそうだ。気が付くのが遅い。

小ネタ

ボストンの公園で機雷が飾ってあった。隣にいたおじさんにこれ何？って聞いたら It's a mine. と教えてくれた。そこですかさず Oh It's yours! と返せばよかったと後で後悔。