Question

我遇到了一个错误，使用Eigen和SparseMatrix和Visual Studio Pro 2013我不明白

连续修剪产品在DEBUG模式下完美运行，但它们很长。我想切换到RELEASE，而不更改代码中的任何内容，但程序正在崩溃这种模式。

我遇到访问冲突错误c0000005（法国VS，对不起）：

Message de résultat  :  Code de l'exception : C0000005
StackTrace de résultat  :   
à Eigen::internal::sparse_sparse_product_with_pruning_impl<Eigen::SparseMatrix<short,0,int>,Eigen::SparseMatrix<short,0,int>,Eigen::SparseMatrix<short,0,int> >() dans [..]\eigen\src\sparsecore\sparsesparseproductwithpruning.h:ligne 61
    à Eigen::internal::sparse_sparse_product_with_pruning_selector<Eigen::SparseMatrix<short,1,int>,Eigen::SparseMatrix<short,1,int>,Eigen::SparseMatrix<short,0,int>,1,1,0>::run() dans [..]\eigen\src\sparsecore\sparsesparseproductwithpruning.h:ligne 143
    à Eigen::internal::unary_evaluator<Eigen::SparseView<Eigen::Product<Eigen::SparseMatrix<short,1,int>,Eigen::SparseMatrix<short,1,int>,2> >,Eigen::internal::IteratorBased,short>::unary_evaluator<Eigen::SparseView<Eigen::Product<Eigen::SparseMatrix<short,1,int>,Eigen::SparseMatrix<short,1,int>,2> >,Eigen::internal::IteratorBased,short>() dans [..]\eigen\src\sparsecore\sparseproduct.h:ligne 158
    à Eigen::SparseMatrix<short,1,int>::operator=<Eigen::SparseView<Eigen::Product<Eigen::SparseMatrix<short,1,int>,Eigen::SparseMatrix<short,1,int>,2> > >() dans [..]\eigen\src\sparsecore\sparsematrix.h:ligne 1080
    à IncMethod_Tester::EigenMatrix_Tests::NullMultiplication() dans [...]\integrity_tests.cpp:ligne 856

在integrity_test中测试的代码如下：

        SparseMatrix<short, RowMajor> A(0, 1);
        SparseMatrix<short, RowMajor> B(1, 5);
        SparseMatrix<short, RowMajor> C(5, 4);

        std::vector<Triplet<short>> values;
        values.push_back(Triplet<short>(0, 4, 1));
        values.push_back(Triplet<short>(0, 3, 1));
        values.push_back(Triplet<short>(0, 2, 1));
        B.setFromTriplets(values.begin(), values.end());

        std::vector<Triplet<short>> values2;
        values2.push_back(Triplet<short>(0, 0, 1));
        values2.push_back(Triplet<short>(1, 1, 1));
        values2.push_back(Triplet<short>(2, 2, -1));
        values2.push_back(Triplet<short>(2, 3, -1));
        values2.push_back(Triplet<short>(3, 2, 1));
        values2.push_back(Triplet<short>(4, 3, 1));
        C.setFromTriplets(values2.begin(), values2.end());


        std::fstream f("NULLMULTIPLICATION", std::fstream::out);
        f << "A MATRIX (" << A.rows() << "," << A.cols() << ")" <<std::endl << A << std::endl;
        f << "B MATRIX (" << B.rows() << "," << B.cols() << ")" <<std::endl<< B << std::endl;
        f << "C MATRIX (" << C.rows() << "," << C.cols() << ")" << std::endl << C << std::endl;

        B = (A*B).pruned();
        f << "AB MATRIX (" << B.rows() << "," << B.cols() << ")" << std::endl << A << std::endl;
        B = (B*C).pruned();
        f << "B MATRIX (" << B.rows() << "," << B.cols() << ")" << std::endl << B << std::endl;

NULLMULTIPLICATION文件中的输出是：

A MATRIX (0,1)

B MATRIX (1,5)
0 0 1 1 1 

C MATRIX (5,4)
1 0 0 0 
0 1 0 0 
0 0 -1 -1 
0 0 1 0 
0 0 0 1 

AB MATRIX (0,5)

这意味着该程序正在运行：

B = (A*B).pruned();
B = (B*C).pruned();

我在这里做错了吗？我尝试创建临时变量，但它仍然崩溃：

SparseMatrix<short, RowMajor> D = (A*B).pruned();
SparseMatrix<short, RowMajor> E = (D*C).pruned();

对此产品有何建议？我对多个稀疏产品的语法可能有误？

提前致谢。

Answer 1

好的，我找到了导致错误的原因。这是由于修剪产品中的0维度，如果任何涉及的矩阵具有零维度，我只是通过在产品之前进行测试来解决它，如果是这样的话，只需做一个没有修剪的普通产品。

连续修剪过的产品导致EIGEN崩溃

1 个答案: