python - 如何找到 scipy.integrate.ode 的默认 atol 和 rtol？

Question

如何在 scipy.integrate.ode 中找到默认参数？具体来说，atol 和 rtol in integrator dopri5？可以使用 set_integrator 方法设置参数，但如果我不设置参数，我看不到任何查询它的方法以了解它一直在使用什么。

Answer 1

据我所知，找出答案的唯一方法是调查 the source .在那里，从第 966 行开始(截至目前)，您会发现:

class dopri5(IntegratorBase):

    […]

    def __init__(self,
                 rtol=1e-6, atol=1e-12,
                 nsteps=500,
                 max_step=0.0,
                 […]
                 ):

因此，dopri 的默认绝对容差为 10⁻¹²，其默认相对容差为 10⁻⁶。 (是的，这应该记录在案。)