переданные параметры для повышения odeint в переполнении стека

Question

переданные параметры для повышения odeint в переполнении стека

это Ответ полезен, но я хотел бы знать, как передать несколько параметров разных типов в модель ODE, возможно, в виде структуры. Для моего непосредственного случая использования, я должен быть в состоянии передать один std::array<double, 6>, два std::vector<std::vector<double>> и два два double скаляры для четырех параметров. В связанном примере, а также в 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;
};

3

boost c++integrator odeint parameters

Решение

Другие решения

Других решений пока нет …

Источник

Accepted Answer

Я полагаю, что придумала структурное решение, которое работает, но я не уверена, имеет ли оно какую-либо переменную / область памяти. Вот пример:

#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);

}

1