本征 3 方阵除以列阵

Question

我想将方形数组 (b) 的每一行除以列数组 (a) 的转置。

测试代码也在godbolt.org，里面安装了eigen 3库

预期结果：

// before (a)
2
2
2
// before (b)
2 2 2
4 4 4
6 6 6
// after
1 1 1
2 2 2
3 3 3

测试：

#include <Eigen/Eigen>
#include <iostream>
using namespace Eigen;

int main() {
    ArrayXXf a(3, 1);
    a << 2, 2, 2;
    ArrayXXf b(3, 3);
    b << 2, 2, 2, 4, 4, 4, 6, 6, 6;
    std::cout << a << "\n";
    std::cout << b << "\n";
    b.rowwise() /= a.transpose();
    std::cout << b << "\n";
}

编译错误

...
error: static assertion failed: YOU_TRIED_CALLING_A_VECTOR_METHOD_ON_A_MATRIX
...

相关问题：

接受的答案是对转换为特征数组的特征矩阵使用行向。我只是不明白为什么它在我的情况下不起作用...

Answer 1

该断言告诉您您正在使用矩阵（又名二维数组），而预期使用编译时向量（又名一维数组）。所以解决方案是将 a 定义为编译时一维数组：

ArrayXf a(3);

你也可以从标准线性代数的角度来看这个运算：

VectorXd a(3);
MatrixXd b(3,3);
b = b * a.asDiagonal().inverse();

本征 3 方阵除以列阵

Eigen 3 Square Array Divided by Column Array

c++

eigen

eigen3