Я рву на себе волосы, пытаясь получить доступ к столбцу данных в Eigen::Tensor<double, 3>
как Eigen::VectorXd
.
Нарезка, согласно этому ответу, отлично работает, чтобы получить колонка хочу. Но я не могу назначить это вектору.
Что я имею:
Eigen::Tensor<double, 3> my_tens(2, 3, 4);
my_tens.setRandom();
Eigen::array<Eigen::Index, 3> dims = my_tens.dimensions();
Eigen::array<Eigen::Index, 3> offsets = {0, 1, 0};
Eigen::array<Eigen::Index, 3> extents = {dims[0], 1, 1};
// This works perfectly, and is exactly the column I want:
std::cout << my_tens.slice(offsets, extents);
// ... and I want it as a VectorXd, so:
Eigen::VectorXd my_vec(dims[0]);
Следующие вещи, которые я пробовал, все FAIL:
// Direct assignment (won't compile, no viable overloaded '=')
my_vec = my_tens.slice(offsets, extents);
// Initialisation (won't compile, no viable overloaded '<<')
my_vec << my_tens.slice(offsets, extents);
// Same problem with reshaping:
g_a = signature_a.g.slice(offsets, extents).reshape(Eigen::array<Eigen::Index, 2>{dims[0], 1});
// Converting the base (won't compile, no member 'matrix')
my_vec << my_tens.slice(offsets, extents).matrix();
Я также пытался сопоставить как в этом ответе, но это тоже не работает (EDIT: я думал, что это связано с порядком хранения, но на самом деле это связано с неправильным смещением, см. мой ответ):
// This produces a part of a row of the tensor, not a column. Gah!
auto page_offset = offsets[2] * dims[0] * dims[1];
auto col_offset = offsets[1] * dims[0];
auto bytes_offset = sizeof(double) * (page_offset + col_offset)
Eigen::Map<Eigen::VectorXd> my_mapped_vec(my_tens.data() + bytes_offset, dims[0]);
Это действительно должно быть так сложно или я упускаю что-то простое? Спасибо за любую помощь!