传递参数以增强集成符
pass parameter to boost integrator
我已经阅读了这些答案1,2关于将参数传递到odeint(boost)。我在这里尝试进行相同的过程。但是代码给出了错误的答案。这是我的尝试:这是Lorenz振荡器Lyapunov指数的增强计算的简化代码。
#include <iostream>
#include <boost/array.hpp>
#include <boost/numeric/odeint.hpp>
#include "gram_schmidt.hpp"
using namespace std;
using namespace boost::numeric::odeint;
typedef vector<double> dim1;
const double sigma = 10.0;
const double R = 28.0;
const double b = 8.0 / 3.0;
const size_t n = 3;
const size_t num_of_lyap = 3;
const size_t N = n + n * num_of_lyap;
// system with out perturbation
struct lorenz
{
void operator()( const dim1 &x , dim1 &dxdt , double t ) const
{
dxdt[0] = sigma * ( x[1] - x[0] );
dxdt[1] = R * x[0] - x[1] - x[0] * x[2];
dxdt[2] = -b * x[2] + x[0] * x[1];
}
};
// system with perturbation
void lorenz_with_lyap( const dim1 &x , dim1 &dxdt , double t )
{
lorenz()( x , dxdt , t );
for( int l=0 ; l<num_of_lyap ; ++l )
{
const double *pert = &x[3 + l * 3];
double *dpert = &dxdt[3 + l * 3];
dpert[0] = - sigma * pert[0] + 10.0 * pert[1];
dpert[1] = ( R - x[2] ) * pert[0] - pert[1] - x[0] * pert[2];
dpert[2] = x[1] * pert[0] + x[0] * pert[1] - b * pert[2];
}
}
//------------------------------------------------------------------------
int main( int argc , char **argv )
{
const double dt = 0.01;
dim1 x(N);
x[0] = 10.; x[1] = 10.0; x[2] = 5.0; // initial condition
dim1 lyap(3);
runge_kutta4<dim1, double, dim1, double, range_algebra> rk4;
// perform 10000 transient steps
integrate_n_steps(rk4, lorenz(), x, 0.0, dt, 10000);
fill( x.begin()+n , x.end() , 0.0 );
for( size_t i=0 ; i<num_of_lyap ; ++i ) x[n+n*i+i] = 1.0;
fill( lyap.begin() , lyap.end() , 0.0 );
double t = 0.0;
size_t count = 0;
while( true )
{
t = integrate_n_steps( rk4 , lorenz_with_lyap , x , t , dt , 100 );
gram_schmidt< num_of_lyap >( x , lyap , n );
++count;
if( !(count % 1000) )
{
cout << t;
for( size_t i=0 ; i<num_of_lyap ; ++i ) cout << "t" << lyap[i] / t ;
cout << endl;
}
}
return 0;
}
首先我定义了这些类
class lorenz
{
public:
void operator()...
和
class lorenz_with_lyap
{
public:
void operator()(const dim1 &x, dim1 &dxdt, double t) const
{
dxdt[0] = sigma * (x[1] - x[0]);
dxdt[1] = R * x[0] - x[1] - x[0] * x[2];
dxdt[2] = -b * x[2] + x[0] * x[1];
for (int l = 0; l < num_of_lyap; ++l)
{
const double *pert = &x[3 + l * 3];
double *dpert = &dxdt[3 + l * 3];
dpert[0] = -sigma * pert[0] + sigma * pert[1];
dpert[1] = (R - x[2]) * pert[0] - pert[1] - x[0] * pert[2];
dpert[2] = x[1] * pert[0] + x[0] * pert[1] - b * pert[2];
}
}
};
然后我可以更改它以定义一个参数:
class lorenz
{
double sigma;
public:
lorenz(double si) : sigma(si) { }
void operator()...
但是定义lorenz_with_lyap
类会更改结果,我不知道我在哪里做错了什么。
你有什么想法吗?(现在这是一项家庭作业,我将使用另一种方法检查此代码的结果)。谢谢。
挣扎后,我终于对此进行了更正,如果有人可能有兴趣,请在此处发布:
class lorenz
{
public:
void operator()(const dim1 &x, dim1 &dxdt, double t) const
{
dxdt[0] = sigma * (x[1] - x[0]);
dxdt[1] = R * x[0] - x[1] - x[0] * x[2];
dxdt[2] = -b * x[2] + x[0] * x[1];
}
};
class lorenz_with_lyap
{
const double sig;
public:
lorenz_with_lyap(double si) : sig(si) { }
void operator()(const dim1 &x, dim1 &dxdt, double t) const
{
dxdt[0] = sig * (x[1] - x[0]);
dxdt[1] = R * x[0] - x[1] - x[0] * x[2];
dxdt[2] = -b * x[2] + x[0] * x[1];
for (int l = 0; l < num_of_lyap; ++l)
{
const double *pert = &x[3 + l * 3];
double *dpert = &dxdt[3 + l * 3];
dpert[0] = -sig * pert[0] + 10.0 * pert[1];
dpert[1] = (R - x[2]) * pert[0] - pert[1] - x[0] * pert[2];
dpert[2] = x[1] * pert[0] + x[0] * pert[1] - b * pert[2];
}
}
};
//------------------------------------------------------------------------
// in main
lorenz_with_lyap lo(10.0);
t = integrate_n_steps(rk4, lo, x, t, dt, 100);
...
相关文章:
- 将公共但非静态的成员函数与ALGLIB集成
- 无法将结构注册为增强几何体3D点
- 将IBM Rhapsody模型集成到VS 2019中
- 从R调用C++函数并对其进行集成时出错
- 增强基于 XML class_id的反序列化
- 增强精神解析器规则以检测语句中的特殊结尾
- C++ 使用增强正则表达式库时断言崩溃
- 增强精神 X3:错误:在"..."中没有名为"大小"的类型
- 如何集成 HID USB 控制器?
- 如何包装(撰写)增强 hana 地图并访问括号运算符(运算符 [])?
- 有没有办法在C ++中检测汉字?(使用增强)
- 如何将增强MSM与增强信号结合使用?
- 在VS2019项目中集成ImageMagick:x64-windows-static library
- 将Qt集成到现有的VS项目中以取代WinAPI
- 将 Crashpad 与 Windows Qt 应用程序集成
- Python3.6 模板中的 CGAL C++ 集成错误
- 如何使用增强::几何计算多边形的旋转固体体积?
- 复制部分文件的C++/增强方式
- 增强 ASIO 和串行端口异步读取
- 传递参数以增强集成符