inline void multiply1(
    Eigen::MatrixXd &m,
    std::vector<double> &vector,
    std::vector<double> &result
)
    {
        int size = (int)vector.size();
        result.resize(size);
        Eigen::Map<Eigen::VectorXd> Evector(vector.data(), size);

        Eigen::VectorXd a;
        a.noalias() = m * Evector;

        Eigen::VectorXd::Map(&result[0], a.size()) = a;
    }