Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

Initialize an Eigen::MatrixXd from a 2d std::vector

Tags:

c++

eigen

This should hopefully be pretty simple but i cannot find a way to do it in the Eigen documentation.

Say i have a 2D vector, ie

std::vector<std::vector<double> > data

Assume it is filled with 10 x 4 data set.

How can I use this data to fill out an Eigen::MatrixXd mat.

The obvious way is to use a for loop like this:

#Pseudo code
Eigen::MatrixXd mat(10, 4);
for i : 1 -> 10
   mat(i, 0) = data[i][0];
   mat(i, 1) = data[i][1];
   ...
 end

But there should be a better way that is native to Eigen?

like image 346
Fantastic Mr Fox Avatar asked Sep 16 '13 23:09

Fantastic Mr Fox


1 Answers

Sure thing. You can't do the entire matrix at once, because vector<vector> stores single rows in contiguous memory, but successive rows may not be contiguous. But you don't need to assign all elements of a row:

std::vector<std::vector<double> > data;
MatrixXd mat(10, 4);
for (int i = 0; i < 10; i++)
  mat.row(i) = VectorXd::Map(&data[i][0],data[i].size());
like image 63
us2012 Avatar answered Nov 02 '22 22:11

us2012