594びぼーろく; 只今絶賛計算中

重い腰を上げてmodified gradient flow method を使ったpythonコードを書き始めた。

この2,3ヶ月かけて色々やってみて結構面倒なことに気づいた。

しかも適当にやると二階微分が壊れる。一般にflow timeの間隔は空間の分割の2乗より小さく取る必要があるけど、どうやらものによって違うらしい。例えばx^2f''みたいなパターンの場合は2乗ではうまく行かない。加えて、与えられているパラメータより空間分割は十分に小さく取る必要があるようだ。

結論：面倒くさい。