Этот ответ полезен, но я хотел бы знать, как передать несколько параметры разных типов модели 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;
};
