Eigen3 colwise与另一个矩阵交叉乘积

Question

我正在尝试将以下代码编写为一个特征表达式：

Eigen::Matrix<double,3,2> a, b, c(3,2);
a << 1,2,3,4,5,6;
b << 7,8,9,10,11,12;
for(int i=0; i<2; ++i)
  c.col(i) = a.col(i).cross(b.col(i));

即。我想写点像

c = a.colwise().cross(b.colwise())

但这不能用Eigen 3.2.7编译。我发现我可以写

c = a.colwise().cross(b.col(0))

但这不是我想要的......

非常感谢您的帮助。在我的FEM代码中，a和b的列数不是2，而是大约1e6，效率非常重要。

提前多多感谢， 拉斐尔