This https://stackoverflow.com/questions/33856910/passing-runtime-parameters-to-odeint-integrator/答案很有帮助,但我想知道如何将不同类型的多个参数传递给 ODE 模型,也许是在结构中。对于我的直接用例,我需要能够通过一个std::array<double, 6>
, two std::vector<std::vector<double>>
和两个两个double
总共要传递四个参数的标量。在链接的示例中以及谐波振荡器.cpp https://github.com/headmyshoulder/odeint-v2/blob/master/examples/harmonic_oscillator.cpp,只有一个double
传递的参数。谢谢。
下面是我需要传递到 ODE 力模型并在速率方程中使用的结构示例。
struct T
{
std::array<double, 6> IC;
double S;
double M;
std::vector<std::vector<double>> C;
std::vector<std::vector<double>> WT;
};
我相信我已经想出了一个可行的结构解决方案,但不确定它是否有任何变量/内存范围禁忌。这是一个例子:
#include <vector>
#include <boost/numeric/odeint.hpp>
// define structure
struct T
{
std::array<double, 6> IC;
double S;
};
// force model
class harm_osc
{
struct T T1;
public:
harm_osc(struct T G) : T1(G) {}
void operator() ( const std::vector< double > &x , std::vector< double > &dxdt , const double /* t */ )
{
dxdt[0] = x[1];
dxdt[1] = -x[0] - T1.IC[0]*x[1] + T1.S;
}
};
// print integrated state solution
void write_solution( const std::vector< double > &x , const double t )
{
printf("%-6.2f %-6.2f %-6.2f\n", t, x[0], x[1]);
}
// problem setup
int main()
{
std::vector< double > x(2);
x[0] = 1.0;
x[1] = 0.0;
struct T T2;
T2.IC = {0.15, 0.15, 0.15, 0.15, 0.15, 0.15};
T2.S = 0.0;
harm_osc ho(T2);
boost::numeric::odeint::integrate(ho, x, 0.0, 10.0, 0.1, write_solution);
}
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)