boost - boost::numeric::odeint 用于刚性系统，如果我们无法得到解析的雅可比矩阵怎么办？

Question

我正在使用 boost::numeric::odeint rosenbrock4 积分器。问题是我无法得到一个分析的雅可比矩阵，有什么解决方案可以解决这个问题吗？

Answer 1

您可以使用数值方案来计算近似雅可比。但这必须由您自己或其他一些库来实现。odeint 不提供此功能。

出于好奇：您找不到雅可比行列式的系统到底是什么？如果 rhs 不可微分，您可能会遇到其他数字问题。