与solve_ivp集成的速度比c语言慢++

Question

使用solve_ivp from整合集成一个刚性系统大约需要3分钟才能完成。积分器是：

sol = solve_ivp(DC_model,[t0,tf],y0,method='LSODA')

用c++和Boost库解决同一个系统需要2.12742秒。这是一个巨大的区别。 有什么方法可以提高python脚本的时间？在

完整代码：

^{pr2}$