У меня две собственные матрицы, и я хотел бы их конкатенировать, как в matlab cat(0, A, B)
Есть ли что-то эквивалентное в собственном?
Спасибо.
У меня две собственные матрицы, и я хотел бы их конкатенировать, как в matlab cat(0, A, B)
Есть ли что-то эквивалентное в собственном?
Спасибо.
Для этого вы можете использовать синтаксис запятой инициализатора.
По горизонтали:
MatrixXd C(A.rows(), A.cols()+B.cols());
C << A, B;
Вертикально:
// eigen uses provided dimensions in declaration to determine
// concatenation direction
MatrixXd D(A.rows()+B.rows(), A.cols()); // <-- D(A.rows() + B.rows(), ...)
D << A, B; // <-- syntax is the same for vertical and horizontal concatenation
Для удобства чтения можно отформатировать вертикальные конкатенации с пробелами:
D << A,
B; // <-- But this is for readability only.
У меня был немного другой вариант использования: для вертикальной компоновки std :: vector из собственных матриц. Вот как я реализовал функцию более общего назначения. Дайте мне знать, если это может быть улучшено:
// matrix_eig = Eigen::MatrixXf in RowMajor format
matrix_eig VStack(const std::vector<matrix_eig> &mat_vec) {
assert(!mat_vec.empty());
long num_cols = mat_vec[0].cols();
size_t num_rows = 0;
for (size_t mat_idx = 0; mat_idx < mat_vec.size(); ++mat_idx) {
assert(mat_vec[mat_idx].cols() == num_cols);
num_rows += mat_vec[mat_idx].rows();
}
matrix_eig vstacked_mat(num_rows, num_cols);
size_t row_offset = 0;
for (size_t mat_idx = 0; mat_idx < mat_vec.size(); ++mat_idx) {
long cur_rows = mat_vec[mat_idx].rows();
vstacked_mat.middleRows(row_offset, cur_rows) = mat_vec[mat_idx];
row_offset += cur_rows;
}
return vstacked_mat;
}
Я бы использовал индексирование Eigen-блоков таким же образом, как этот пост (который объединяет существующую матрицу).
Блочная индексация позволяет избежать двусмысленности направления в принятом подходе и имеет довольно компактный синтаксис:
MatrixXd C(A.rows(), A.cols()+B.cols());
C.leftCols(A.cols()) = A;
C.rightCols(B.cols()) = B;