我目前正在尝试使用odeint和Eigen3来集成一个nBody系统(目标是提供用于行星形成的高级例程的库,例如混合变量辛或MVS的混合变量)。在尝试使用不同的步进器时,我发现当使用普通步进器时,state_type std::vector<Eigen::Vector3d>
工作正常,但是使用受控步进器(例如burlisch_stoer
)编译失败,第一条错误信息是:
/usr/include/boost/numeric/odeint/stepper/controlled_runge_kutta.hpp:87:40: error: cannot convert ‘boost::numeric::odeint::norm_result_type<std::vector<Eigen::Matrix<double, 3, 1> >, void>::type {aka Eigen::Matrix<double, 3, 1>}’ to ‘boost::numeric::odeint::default_error_checker<double, boost::numeric::odeint::range_algebra, boost::numeric::odeint::default_operations>::value_type {aka double}’ in return
return algebra.norm_inf( x_err );
这是否意味着norm_result_type被错误推断?这个规范究竟做了什么?它应该是x_err中找到的最高value_type吗?
和第二个:
/usr/include/boost/numeric/odeint/algebra/range_algebra.hpp:132:89: error: call of overloaded ‘Matrix(int)’ is ambiguous
static_cast< typename norm_result_type<S>::type >( 0 ) );
我是否必须提供自己的代数才能以这种方式使用它?我宁愿不切换到std::vector<double>
或Eigen :: VectorNd,因为分组的坐标对于ODE右侧的可读性非常有利。
这是我使用的代码的简化示例。
#include "boost/numeric/odeint.hpp"
#include "boost/numeric/odeint/external/eigen/eigen.hpp"
#include "Eigen/Core"
using namespace boost::numeric::odeint;
typedef std::vector<Eigen::Vector3d> state_type;
struct f
{
void operator ()(const state_type& state, state_type& change, const double /*time*/)
{
};
};
int main()
{
// Using this compiles
typedef euler <state_type> stepper_euler;
// Using this does not compile
typedef bulirsch_stoer <state_type> stepper_burlisch;
state_type x;
integrate_const(
stepper_burlisch(),
f(),
x,
0.0,
1.0,
0.1
);
return 0;
}
Headmyshoulders解决方案有效。我创建了一个继承自range_algebra
:
class custom_algebra : public boost::numeric::odeint::range_algebra {
public:
template< typename S >
static double norm_inf( const S &s )
{
double norm = 0;
for (const auto& inner : s){
const double tmp = inner.maxCoeff();
if (tmp > norm)
norm = tmp;
}
return norm;
}
} ;
我认为应该可以使用range_algebra
和vector_space_algebra
来创建这样的代数,但我还没有尝试过。
答案 0 :(得分:1)
我担心您的州类型不容易被整合。问题是,这是一个类似范围的状态类型(vector<>
)与类似向量空间的状态类型(Vector3d
)混合在一起。您需要range_algebra
迭代向量。但norm_inf
中的range_algebra
在您的情况下不起作用。
您可以做的是复制range_algebra
并保留所有for_eachX
方法,只重新实施norm_inf
以适合您的州类型。那不应该是困难的。然后你只需通过
typedef bulirsch_stoer< state_type , double , state_type , double , your_algebra > stepper_burlisch