Boost C++ Libraries: Ticket #10680: Cannot use Eigen vector as state type in integrate?

karsten — Thu, 29 Jan 2015 21:41:10 GMT

There are two problems with your code. First, eigen is not support by default, you need to include the eigen adaption of odeint. And secondly, integrate can not determine the value from the eigen vector by itself. If you change the example to

#include "boost/numeric/odeint.hpp"
#include "boost/numeric/odeint/external/eigen/eigen.hpp"
#include "Eigen/Core"
using namespace boost::numeric::odeint;
int main()
{
    Eigen::Vector3d state(0, 0, 0);
    integrate< double >(
        [](const auto& state, auto& change, const double /*time*/){
            change = Eigen::Vector3d(1, 0, 0);
        },
        state,
        0.0,
        1.0,
        0.1
        );
        return 0;
}

everything should work.

status changed; resolution set

karsten — Thu, 29 Jan 2015 21:41:28 GMT

status new → closed
resolution → wontfix

anonymous — Thu, 29 Jan 2015 21:53:20 GMT

#include "boost/numeric/odeint/external/eigen/eigen.hpp"

I don't have that file in boost-1.57.0, only eigen_algebra.hpp and eigen_resize.hpp. I just downloaded the bz2 tarball again from sourceforge and still nothing.

status changed; resolution deleted

anonymous — Thu, 29 Jan 2015 21:53:32 GMT

status closed → reopened
resolution wontfix

karsten — Thu, 29 Jan 2015 22:12:58 GMT

Ok, I see. Can you use the odeint version from github: https://github.com/headmyshoulder/odeint-v2. It will be included in 1.58. Be sure to set the include path to odeint-v2/include before the boost include path.

anonymous — Thu, 29 Jan 2015 22:25:59 GMT

Thanks, it worked once I added -std=c++1y compilation args.

status changed; resolution set

mariomulansky — Wed, 10 May 2017 04:18:26 GMT

status reopened → closed
resolution → fixed