的数据为什么当我用Radau方法运行scipy.integrate.solve_ivp时，雅可比矩阵求值的次数不等于函数求值的次数？具体来说，在我的代码中，它只对雅可比矩阵求值一次或两次，尽管雅可比矩阵在每一步都在变化。
toe950271#
雅可比矩阵的计算和因式分解是一个昂贵的操作。为了使常微分方程积分器正确有效地工作，一个近似雅可比矩阵就足够了。求解器在更新雅可比矩阵之前尽可能长时间地计算一次雅可比矩阵。因此，如果雅可比矩阵在解段所在的常微分方程域上变化很小，那么一个雅可比矩阵的计算就足够了。还应注意，对于单个Radau 5步骤，可以有3、6、9或12或甚至3的更高倍数的函数求值次数。
1条答案
按热度按时间toe950271#
雅可比矩阵的计算和因式分解是一个昂贵的操作。为了使常微分方程积分器正确有效地工作，一个近似雅可比矩阵就足够了。求解器在更新雅可比矩阵之前尽可能长时间地计算一次雅可比矩阵。因此，如果雅可比矩阵在解段所在的常微分方程域上变化很小，那么一个雅可比矩阵的计算就足够了。
还应注意，对于单个Radau 5步骤，可以有3、6、9或12或甚至3的更高倍数的函数求值次数。