Eigen3:写出稀疏矩阵

时间:2017-11-13 15:47:05

标签: c++ c++11 eigen eigen3

我使用以下代码在Eigen3中创建了一个稀疏的martrix:

#include <eigen3/Eigen/Sparse>
#include <eigen3/unsupported/Eigen/src/SparseExtra/MarketIO.h>
#include <vector>

using SpMatrix = Eigen::SparseMatrix<double>;
using Eigen::SparseMatrix;
using S = Eigen::Triplet<double>;
using namespace std;
using AdjacenyMatrix = SpMatrix;

int main() {

    vector <S> nonzero_compenents;
    AdjacenyMatrix am(10, 10);
    nonzero_compenents.push_back(move(S(1, 1, 1.0)));

    am.setFromTriplets(nonzero_compenents.begin(), nonzero_compenents.end());

    bool s = saveMarket(am, "/home/morris/Schreibtisch/sparse_matrices");

    return 0;

}

不幸的是,代码无法编译。 MarketIO.h

似乎有问题
In file included from /home/morris/sparse/main.cpp:2:0:
/usr/include/eigen3/unsupported/Eigen/src/SparseExtra/MarketIO.h: In function ‘void Eigen::internal::PutMatrixElt(Scalar, int, int, std::ofstream&)’:
/usr/include/eigen3/unsupported/Eigen/src/SparseExtra/MarketIO.h:87:9: error: no match for ‘operator<<’ (operand types are ‘std::ofstream {aka std::basic_ofstream<char>}’ and ‘int’)
     out << row << " "<< col << " " << value << "\n";
         ^

... 有什么建议吗?

2 个答案:

答案 0 :(得分:0)

正如@ggael所说:不要包含单个文件,而是包含它们的模块/目录。请注意,在这种情况下,您还需要从路径中删除src

此前:

#include <eigen3/unsupported/Eigen/src/SparseExtra/MarketIO.h>

现在:

#include <eigen3/unsupported/Eigen/SparseExtra>

完整代码:

#include <eigen3/Eigen/Sparse>
#include <eigen3/unsupported/Eigen/SparseExtra>
#include <vector>

using SpMatrix = Eigen::SparseMatrix<double>;
using Eigen::SparseMatrix;
using S = Eigen::Triplet<double>;
using namespace std;
using AdjacenyMatrix = SpMatrix;

int main() {

    vector <S> nonzero_compenents;
    AdjacenyMatrix am(10, 10);
    nonzero_compenents.push_back(move(S(1, 1, 1.0)));

    am.setFromTriplets(nonzero_compenents.begin(), nonzero_compenents.end());

    bool s = saveMarket(am, "/home/morris/Schreibtisch/sparse_matrices");

    return 0;

}

现在它应该使用例如

进行编译
g++ -std=c++11 mycode.cc

答案 1 :(得分:0)

使用

#include <fstream>
#include <Eigen/Sparse>
#include <unsupported/Eigen/SparseExtra>

代替

#include <eigen3/Eigen/Sparse>
#include <eigen3/unsupported/Eigen/src/SparseExtra/MarketIO.h>

代码编译

注意#include <eigen3/unsupported/Eigen/SparseExtra>无法编译