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