我在Trilinos LOCA(c ++)中一直是伪弧长连续,非线性求解器的用户多年,但我最近遇到了python包petsc4py。这简化了我的编程,以减少速度。
petc4py唯一缺失的类是伪arclength延续求解器。 PETSC有一个伪时间步进求解器,但这并不像转折点那样有用。
Trilinos LOCA有一个PETSC包装器,但我一直在寻找一个解决方案,在python(来自petsc4py)而不是c或c ++中指定问题。
因此,我正在寻找一个用petsc4py编写的类,或者可以从petsc4py调用,实现伪arclength延续。
非常感谢你。