我想用python中的newton_krylov解算器来求解一个非线性代数方程组。我还有一个函数可以给出系统的雅可比矩阵;它看起来像这样 def jac(x, param1, param2, param3): dim, =x.shape # this is thereturn jacob 此函数在输入时接受一个矢量x,以及与从param4中sys一个零件相同的参数,并返回另一个矢
我正在尝试使用MatCont包来解决一个非常大的非线性问题。由于大量的维度和非线性性质,我相信将我的系统的雅可比提供给MatCont算法将极大地提高速度。然而,我不能让它认识到它有一个雅可比可以使用!作为一个最小的工作示例,我修改了帮助文档中的圆查找器,使其包含一个Jacobian:
function out