I have tried to construct a sparse matrix using set_from_triplets method but when I tried with large data I get this error :

a.out: /usr/local/include/eigen3/Eigen/src/SparseCore/SparseMatrix.h:952: void Eigen::internal::set_from_triplets(const InputIterator&, const InputIterator&, SparseMatrixType&, int) [with InputIterator = __gnu_cxx::__normal_iterator<Eigen::Triplet<long int>*, std::vector<Eigen::Triplet<long int> > >; SparseMatrixType = Eigen::SparseMatrix<long int, 1>]: Assertion `it->row()>=0 && it->row()<mat.rows() && it->col()>=0 && it->col()<mat.cols()' failed.
Aborted


actually I read a data form file and push them to triplet
Could you help me ?

thanks,
Mohammad

you have to resize the sparse matrix (rows/cols) before inserting the elements

ggael wrote:you have to resize the sparse matrix (rows/cols) before inserting the elements

I still get this error
It is my code :

Code: Select all

#include <Eigen/Sparse>
#include <Eigen/SparseCore>
#include <iostream>
#include <sstream>
#include <fstream>
#include <string>
#include <iostream>
#include <cstring>
#include <vector>

typedef Eigen::Triplet<long int> T;

using namespace Eigen;

using namespace std;
int main() {
    int ni = 0,no = 0,Nx = 6955272, Ny = 17354796 ;
    SparseMatrix<long int,RowMajor> S;
         long int value,*B,txid,userid;
         std::vector<T> coefficients;
         coefficients.reserve(60000000);
    std::ifstream infile2("txout__modified.txt");
    std::string line2;
    while (std::getline(infile2,line2))
    {
    std::istringstream iss(line2);
    iss >> txid >> userid >> value;
             coefficients.push_back(T(txid,userid,value));   
        std::cout << "The Line is :"<<no<<"\n";
        no++;}
        Eigen::VectorXd b(Nx);
        S.resize(Nx,Ny);
        S.setFromTriplets(coefficients.begin(), coefficients.end());
       return 0;
}


maybe your file format is not 0-based indexed...

ggael wrote:maybe your file format is not 0-based indexed...

what is 0-based indexed file ?

many file formats assume that the index of the first row/column is 1 while in C++/Eigen it's 0.

But when i run this code for small lines it worked

This assert is telling you that at least one of the coefficient in the triplet list is outside the range you specified. You can check that in your while loop:

txid < Nx && userid < Ny

I have got it.
Thank you.