C ++ Eigen ::传递给Class后无法访问矩阵

Question

我试图通过引用将Eigen :: Matrix传递给类并访问其中的元素。

当我尝试在main函数和类中再次访问矩阵时，如果矩阵大小很大，它将失败并产生分段错误。我使用调试器检查矩阵的值，发现矩阵的元素不可访问。

以下是我的代码：

main.cpp中

Eigen::MatrixXf A = Eigen::MatrixXf::Random(3,640); //this is OK
//Eigen::MatrixXf A = Eigen::MatrixXf::Random(3,640*480); //but not this

std::cout << "in main A col " << A.cols() << " row " << A.rows()
          << "\nA\n" << A.col(100) << "\n\n";

model.testLoadMat(A);

model.testReadMat();

std::cout << "in main testMat col " << model.testMat->cols() << " row " << model.testMat->rows()
          << "\ntestMat\n" << model.testMat->col(100) << "\n\n"; //fails here if A is large

return 0;

model.h

class model
{
    public:
        const Eigen::Matrix<float,3,Eigen::Dynamic> *testMat;
        void testLoadMat(const Eigen::Matrix<float,3,Eigen::Dynamic> &tMat);
        void testReadMat();
}

model.cpp

void model::testLoadMat(const Eigen::Matrix<float,3,Eigen::Dynamic> &tMat)
{
    testMat = &tMat;

    std::cout << "in testLoadMat col " << testMat->cols() << " row " << testMat->rows()
           << "\ntestMat\n" << testMat->col(100) << "\n\n";
}

void model::testReadMat()
{
    std::cout << "in testReadMat col " << testMat->cols() << " row " << testMat->rows()
           << "\ntestMat\n" << testMat->col(100) << "\n\n"; //fails here if A is large
}

他们仍然提供正确的行数和列数，但我无法访问其中的元素。

为什么它不起作用，我该如何解决？否则，有没有更好的方法来做到这一点？

另外，我注意到执行testLoadMat()的时间随矩阵大小而增加。这是否意味着我正在做错误的“通过引用传递”？

Answer 1

在A中main的类型是Eigen::MatrixXf，或Matrix<float, Eigen::Dynamic, Eigen::Dynamic>。你的testLoadMat函数想要引用Eigen::Matrix<float,3,Eigen::Dynamic>对象。由于这些是不同的类型，因此会创建一个正确类型的临时Matrix，它是对传递给testLoadMat的临时对象的引用。然后testLoadMat存储一个指向此临时Matrix的指针，当函数返回并且临时Matrix被销毁时，该指针变为无效。当您调用testReadMat时，您取消引用此无效指针，从而导致未定义的行为。

修复方法是更改​​testLoadMat函数的参数类型以匹配您传入的类型。