Eigen 3 Square Array, Разделенный на столбец Array

Question

Я хочу разделить каждую строку квадратного массива (b) путем транспонирования массива столбцов (a).

Тестовый код также находится на https://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

Утверждение говорит вам, что вы используете матрицу (она же 2D массив), тогда как ожидается вектор времени компиляции (он же 1D массив). Таким образом, решение состоит в том, чтобы определить a как 1D массив во время компиляции:

ArrayXf a(3);

Вы также можете увидеть эту операцию в терминах стандартной линейной алгебры:

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