Open aw02m opened 1 year ago
分岐計算プログラムまで作りました. 設計としてはnewton法に相当する関数newton(func, x0, args=())を用意してあり,funcは目的関数のベクトルFおよびそのヤコビアンJを返すことを満たす関数を与えます. これによりfixやG,PD,NS用の関数を用意してそれを与えるだけで良く,newton法本体の部分はscipy.optimize.root()のように使い回しが良い形式になります.(使用者側はbif.py内のobjective_func()を変更する,すなわちNewton法の目的関数を準備するだけとなります.) まだ数値のハードコーディングをしている部分がありますが,とりあえず共有ということでお願いします. 設計に意見等あればぜひよろしくお願いします.
SymPy は使わない感じですかね?あるいは学習を意識して最初は SymPy にしてないとか?
Sympyは導入予定ですが,とりあえずbifが動いたのでpushしたという感じですね.
分岐計算プログラムまで作りました. 設計としてはnewton法に相当する関数newton(func, x0, args=())を用意してあり,funcは目的関数のベクトルFおよびそのヤコビアンJを返すことを満たす関数を与えます. これによりfixやG,PD,NS用の関数を用意してそれを与えるだけで良く,newton法本体の部分はscipy.optimize.root()のように使い回しが良い形式になります.(使用者側はbif.py内のobjective_func()を変更する,すなわちNewton法の目的関数を準備するだけとなります.) まだ数値のハードコーディングをしている部分がありますが,とりあえず共有ということでお願いします. 設計に意見等あればぜひよろしくお願いします.