c++ - 使用 runge_kutta4 的高维数组

Question

我想使用 boost::numeric::odeint::runge_kutta4 求解耦合微分方程组。它是一个 3D 晶格系统，所以我很自然(也很方便)使用 3D 阵列。有没有办法让 runge_kutta4 使用用户定义的数据结构或 boost multi_array ？

Answer 1

原则上这是可能的。 odeint 提供了一种使用自定义数据结构的机制——代数和运算。看看here .您要么使用现有的代数之一，并尝试调整您的数据结构以使用该代数。或者您实现自己的代数并用它实例化 Runge Kutta 步进器。

您可能还想看看 Eigen、Mtl4、boost.ublas 或 Armadillo 等库。它们可能具有高阶张量的数据类型。例如，Eigen 与 odeint 配合得很好。