Question

system of equations

您好。我想在时间上从零到10 ^ 16和初始条件x（0）= 10 ^ 8和y（0）= 0.5进化这些方程。由于方程式对分母中x的依赖性，我认为使用odeint和runge_kutta_dopri5是一个很好的选择，因为自适应步控制。问题是我不知道如何在实践中做到这一点因为我对c ++和odeint没什么经验。我搜索了很多关于使用odeint但是对我没有帮助的例子。另外我想在x达到零时停止计算我看到了这个https://stackoverflow.com/questions/33334073/stop-integration-in-odeint-with-stiff-ode

基于我迄今为止写的这个例子没有运气

#include <iostream>
#include <vector>
#include <cmath>
#include <boost/array.hpp>
#include <boost/numeric/odeint.hpp>

using namespace std;
using namespace boost::numeric::odeint;

const double b = 43.0e17;

typedef boost::array< double , 2 > state_type;

void binary(const state_type &x , state_type &dxdt , double t )
{
 dxdt[0] = -b*(64.0/5)*(1 + (73.0/24)*pow(x[1],2) 
   +  37.0/96)*pow(x[1],4) )/pow(x[0],3)*pow(1-pow(x[1],2),7.0/2);

 dxdt[1] = -b*(304.0/96)*x[1]*(1 + (121.0/304)*pow(x[1],2))
 /pow(x[0],4)*pow((1 - pow(x[1],2)),5.0/2);

 }

 void write_binary( const state_type &x , const double t )
{
cout << t << '\t' << x[0] << '\t' << x[1] << '\t' << x[2] << endl;
}
 //I dont know what this does but the examples used it
struct streaming_observer
{
 std::ostream& m_out;

streaming_observer( std::ostream &out ) : m_out( out ) { }

template< class State , class Time >
void operator()( const State &x , Time t ) const
{
    m_out << t;
    for( size_t i=0 ; i<x.size() ; ++i ) m_out << "\t" << x[i] ;
    m_out << "\n";
 }
};
 //This was a first try with a given stepper but i want to replace it
 int main(int argc, char **argv)
  {
state_type x = { 20.871e8 , 0.5  }; // initial conditions
integrate( binary , x , 0.0 , 1000.0 , 0.1 , write_binary );
}

当我编译运行时，我收到了此错误

内部程序错误 - 在const T＆amp;中断言（i

我怎样才能完成这项工作？

Answer 1

write_binary函数写入数组边界并导致断言。 x[2]无效。

使用dopri5和odeint boost库

1 个答案: