特征垂直堆叠行到矩阵中

Question

我想创建一个大小2N x 9矩阵，其中N是通过垂直堆叠矩阵2N1x9动态值。

这是我尝试做的事情。

using CoefficientMatrix = Eigen::Matrix<T, Eigen::Dynamic, 9>;
using CoefficientRow = Eigen::Matrix<T, 1, 9>;
CoefficientMatrix A(2*N, 9);
for (int i = 0; i < N; i++) {
CoefficientRow ax;
CoefficientRow ay;
// fill in ax and ay
A << ax, ay;
}

但是，我收到以下运行时错误。

Assertion failed: (((m_row+m_currentBlockRows) == m_xpr.rows() || m_xpr.cols() == 0) && m_col == m_xpr.cols() && "Too few coefficients passed to comma initializer (operator<<)"), function finished, file /usr/local/include/eigen3/Eigen/src/Core/CommaInitializer.h, line 120.

我尝试解析断言语法，但我不确定这些内部变量名称在我的代码中指的是什么(Eigen 的新功能)。

感谢您的任何帮助。

Answer 1

TLDR：写这样的东西：

CoefficientMatrix A(2*N, 9);
for (int i = 0; i < N; i++) {
CoefficientRow ax;
CoefficientRow ay;
// fill in ax and ay
A.row(2*i)   = ax;
A.row(2*i+1) = ay;
}

---

您错误的原因是(如 Avi 解释的那样)operator<<旨在一次填充整个矩阵。事实上，调用operator<<(Array &A, Array const &b)会将b分配给A的左上角，并返回一个代理对象，该对象包含对A的引用，并跟踪已经分配给A的条目数量(存储在m_row、m_currentBlockRowsm_col中)并重载operator,将下一个表达式分配给a的相应位置并相应地增加位置。最后，当该代理对象被销毁时(通常发生在";")，析构函数检查A的所有条目是否已填充(如果没有，则引发失败的断言)。

如果您更喜欢使用<< ,语法，您还可以编写：

A.middleRows<2>(2*i) << ax, ay;

启用优化后，应生成与上述简单实现相同的代码(因此请选择您更容易阅读的代码)。

---

注意：从技术上讲，你可以在循环中使用CommaInitializer，方法是在循环外构造它，将其分配给一个变量，然后只在循环中使用,运算符。我故意不会提供有关如何执行此操作的更多详细信息......

Answer 2

operator <<的重载一次填充整个矩阵。如果维度不匹配，则会收到遇到的运行时错误。

这样看。假设代码可以按预期工作(operator <<插入一行？一次插入两行？)。如何连续两次调用operator <<？现在每个矩阵都必须跟踪调用operator <<的次数吗？或者相反，如果插入了部分行，应该如何处理？