Я только начал использовать Boost Odeint для интеграции системы ODE. Для удобства я хотел бы использовать его с Armadillo, так как обе современные библиотеки C++ с удобным API. Однако, если я укажу arma::vec
как тип состояния, сразу получается, что на первом этапе интеграции integrate_adaptive()
изменяет размер вектора состояния на 0x1
. Я отправляю тривиальный пример здесь:Armadillo конфликтует с Boost Odeint: Odeint изменяет вектор состояния на ноль во время интеграции
#include <iostream>
#include <armadillo>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace arma;
using namespace boost::numeric::odeint;
typedef vec state_type;
class harm_osc
{
private:
mat my_A;
public:
harm_osc(double gam)
{
my_A = { {0.0, 1.0}, {-gam*gam, 0.0} };
}
harm_osc()
{
my_A = { {0.0, 1.0}, {-1.0, 0.0} };
}
void operator() (const vec& x, vec& dxdt, const double t)
{
cout << "size of x: " << size(x) << endl;
cout << "size of dxdt: " << size(dxdt) << endl;
dxdt = my_A*x;
}
};
class observer
{
private:
mat& my_states;
vec& my_times;
public:
observer(mat& states, vec& times):
my_states(states),
my_times(times)
{}
void operator() (const vec& x, double t)
{
my_states.insert_rows(my_states.n_rows, x);
my_times.insert_rows(my_times.n_rows, t);
cout << t << '\t';
for(auto elem : x)
cout << elem << '\t';
cout << endl;
}
};
typedef runge_kutta_cash_karp54<state_type> error_stepper_type;
typedef controlled_runge_kutta<error_stepper_type> controlled_stepper_type;
int main()
{
state_type x = {0.0, 1.0};
vec t;
mat x_full;
integrate_adaptive(make_controlled<error_stepper_type>(1e-5, 1e-5), harm_osc(1.0), x, 0.0, 200.0, 0.01, observer(x_full, t));
}
Если я указать arma::vec::fixed<2>
вместо arma::vec
как state_type
, это просто демо работает правильно. Моя проблема в том, что в текущем проекте, над которым я работаю, я не знаю размер моего вектора состояния во время компиляции, поэтому я не могу исправить его с помощью параметра шаблона, как упоминалось ранее.
Есть ли решение использовать Armadillo с Boost Odeint без фиксации размера вектора состояния во время компиляции? В других частях моего проекта я могу правильно использовать Армадилло. Я хотел бы использовать его во всем моем проекте. Прямо сейчас, при интеграции ODE, я использую Boost Ublas для определения системы ODE.