我使用boost / odeint来解决一阶颂歌。
raman_system.h
中系统的定义如下:
class raman_signal_system {
double alpha_s;
double dB_to_neper;
double g_pump_b;
double P_p;
public:
raman_signal_system()
: alpha_s(0.00021181),
dB_to_neper(0.23025850929940458),
g_pump_b(0.00044446),
P_p(0.494) {}
void operator() (const double P_s, double& dPsdz, const double z)
{
dPsdz = -alpha_s * dB_to_neper * P_s + P_s * P_p * g_pump_b;
}
};
我在raman.cpp中的main
方法如下:
/**
* @author Harun Kara
* @date 23.01.2017
* @brief Implementaion of the Raman amplification for fiber communications systems
*/
#include "raman_system.h"
#include <fstream>
#include <iostream>
#include <vector>
#include<boost/numeric/odeint.hpp>
int main(int argc, char *argv[]) {
namespace odeint = boost::numeric::odeint;
odeint::runge_kutta4< double > rk4;
raman_signal_system sys();
double P_s = 0.00079433;
double z = 0.0; // initial (forward/signal)
double l = 80e3; // Fiberlenght in m
size_t steps = 200;
const double dz = l/steps; // stepsize
std::vector<double> P_s_vec;
std::vector<double> places;
P_s_vec.push_back(P_s);
places.push_back(z);
for( size_t i=0; i<steps; ++i ) {
rk4.do_step( sys, P_s, z, dz );
z += dz;
P_s_vec.push_back(P_s);
places.push_back(z);
}
return 0;
}
如果我尝试构建,我会收到以下错误:
/home/kara/raman/src/raman.cpp:43: required from here
/home/kara/raman/include/boost/numeric/odeint/stepper/base/explicit_stepper_base.hpp:229: error: too many arguments to function
sys( x , m_dxdt.m_v ,t );
我无法找到导致此错误的原因。请有人帮助我,找出我的代码有什么问题?
使用Boost 1.60.0,Qt,CMake
答案 0 :(得分:1)
这一行
raman_signal_system sys();
被编译器识别为本地函数声明,因此错误消息将太多参数传递给该函数。
要使用默认构造函数初始化变量,只需编写
即可raman_signal_system sys; // <<<< No parenthesis
正确声明sys
变量后,重载的调用operator()
应该可以正常工作:
sys( x , m_dxdt.m_v ,t );