Ошибка собственного утверждения при преобразовании из MatrixXf в ArrayXf - PullRequest
Новая
       98

Ошибка собственного утверждения при преобразовании из MatrixXf в ArrayXf

0 голосов
/ 28 марта 2020

Я написал простую программу, которая добавляет скаляр к матрице Random. 

#include <Eigen/Dense>

int main() {
    Eigen::MatrixXf mat = Eigen::MatrixXf::Random(100, 100);
    Eigen::ArrayXf arr = mat.array() + 1;
}

Программа компилируется без ошибок. Однако при выполнении этой строки: 

Eigen::ArrayXf arr = mat.array() + 1;

я получаю следующую ошибку: 

Assertion failed: (other.rows() == 1 || other.cols() == 1), function resizeLike, file /usr/local/Cellar/eigen/3.3.7/include/eigen3/Eigen/src/Core/PlainObjectBase.h, line 374.

Я читаю официальную документацию (https://eigen.tuxfamily.org/dox/group__TutorialArrayClass.html), и я надеваю не понимаю, что я делаю не так.

Ответы [ 3 ]

0 голосов
/ 29 марта 2020

Решил это. Видимо, мне пришлось хранить результат сложения в MatrixXf, а не в ArrayXf. 

#include <Eigen/Dense>

int main() {
    Eigen::MatrixXf mat = Eigen::MatrixXf::Random(100, 100);
    Eigen::MatrixXf arr = mat.array() + 1;
}
0 голосов
/ 29 марта 2020

Соответствующий Array объект MatrixXf (он же Matrix<float, Dynamic, Dynamic>) будет ArrayXXf (он же Array<float, Dynamic, Dynamic>). ArrayXf - это typedef для Array<float, Dynamic, 1>, похожий на VectorXf (он же Matrix<float, Dynamic, 1>).

Следующее должно работать как положено: 

#include <Eigen/Core>

int main() {
    Eigen::MatrixXf mat = Eigen::MatrixXf::Random(100, 100);
    Eigen::ArrayXXf arr = mat.array() + 1;
}
0 голосов
/ 29 марта 2020

Прошло много времени с тех пор, как я что-то сделал с Эйгеном, но я верю, что вы не можете смешивать Массивы и Матрицы подобным образом.

Однако создается впечатление, что вы генерируете матрицу 100x100 со случайными числами, и я не уверен, почему она срабатывает. В настоящий момент у меня нет готового Eigen, но изменение его на использование массивов фиксированного размера, похоже, работает на Godbolt. 

// Type your code here, or load an example.
#include <iostream>
#include <Eigen/Dense>

int main() 
{
    Eigen::Matrix<float, 100, 100> mat = Eigen::Matrix<float, 100, 100>::Random();
    Eigen::Array<float, 100, 100> arr = mat.array() + 1;

    std::cout << "Mat 0,0: " << mat(0,0) << "\nArr 0,0: " << arr(0,0) << "\n";
    std::cout << "Mat 2,0: " << mat(2,0) << "\nArr 2,0: " << arr(2,0);

    return 0;
}

ASM generation compiler returned: 0

Execution build compiler returned: 0

Program returned: 0

Mat 0,0: 0.680375
Arr 0,0: 1.68038
Mat 2,0: 0.566198
Arr 2,0: 1.5662

https://godbolt.org/z/9xu5oW

...