Eigen::Matrix 上的 odeint 在integrate_const 上输出NaN

编程问答 2022-05-01

问题描述

我有以下状态： using state_type= Eigen::Matrix<double,42,1>; 使用定义的 abs 和无限范数，我可以使用以下配置运行它：

typedef runge_kutta_dopri5<state_type,double,state_type,vector_space_algebra> stepper_type;
    auto t0= 0.00,t1= 42.00;
    size_t steps = integrate_adaptive(make_controlled<stepper_type>(1e-5,1e-5),odeState,x,t0,t1,1.00);

但是，当调用类似的东西

size_t steps = integrate_const( stepper_type(),odeSt,0.01,write_state);

我在第一步就已经知道 NaN。我想知道可能是什么原因。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

eigenodeint