Should pass steppers by reference in integrate functions

[anikin@…]

The integrate family of functions takes steppers by value, but steppers like the most useful runge_kutta_dopri5 may contain state which is not the best type to copy around. The worst happens if you try to use SSE-aligned vector for state type, e.g. Eigen::Array. Passing aligned struct by value cannot preserve the alignment and the program crashes. Which makes it impossible to use SSE vectorization with odeint. I think you need to change to passing steppers by reference or add odeint::unwrap_reference< Stepper > wherever steppers are passed, so that boost::ref could be used with them.

[karsten]

I think using odeint::unwrap_reference< Stepper > is the best solution since it does not break any code and one can pass steppers by reference.

[karsten]

It should already work with std::ref or boost::ref