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

Question

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

Answer 1

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

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