#include <thrust/device_vector.h>
#include <thrust/host_vector.h>
#include <boost/numeric/odeint/util/copy.hpp>