私は何年もの間、Trilinos LOCA (c++) で疑似弧長連続、非線形ソルバーのユーザーでしたが、最近、python パッケージ petsc4py に出くわしました。これにより、速度がわずかに低下するためのプログラミングが大幅に簡素化されました。
petsc4py で唯一欠けているクラスは、疑似アークレングス継続ソルバーです。PETSC には疑似時間ステッパー ソルバーがありますが、ターニング ポイントがアルゴリズムを分解するので、これはほとんど役に立ちません。
Trilinos LOCA には PETSC のラッパーがありますが、c や c++ ではなく Python (petsc4py から) で問題が指定されているソリューションを探していました。
したがって、petsc4py で記述された、または petsc4py から呼び出し可能で、疑似弧長継続を実装するクラスを探しています。
どうもありがとうございました。