diff --git a/src/TNL/Matrices/MatrixReader_impl.h b/src/TNL/Matrices/MatrixReader_impl.h index 418e6f5b3eda29a659c4487dfaf34d88c12fea1d..6d6b3eb555150ce516f11781fc12ba41f2ea4e3a 100644 --- a/src/TNL/Matrices/MatrixReader_impl.h +++ b/src/TNL/Matrices/MatrixReader_impl.h @@ -69,7 +69,13 @@ bool MatrixReader< Matrix >::readMtxFileHostMatrix( std::istream& file, if( ! computeCompressedRowLengthsFromMtxFile( file, rowLengths, columns, rows, symmetricMatrix, verbose ) ) return false; + std::cout << " rowLengths sizeof: " << sizeof( rowLengths ) << std::endl; + std::cout << " rowLengths element sizeof: " << sizeof( rowLengths[0] ) << std::endl; + std::cout << " rowLengths getSize(): " << rowLengths.getSize() << std::endl; + matrix.setCompressedRowLengths( rowLengths ); + + std::cout << "->CompressedRowLengths SET" << std::endl; if( ! readMatrixElementsFromMtxFile( file, matrix, symmetricMatrix, verbose, symReader ) ) return false; @@ -340,6 +346,9 @@ bool MatrixReader< Matrix >::readMatrixElementsFromMtxFile( std::istream& file, IndexType processedElements( 0 ); Timer timer; timer.start(); + + std::cout << "\nBefore while..." << std::endl; + while( std::getline( file, line ) ) { if( line[ 0 ] == '%' ) continue; @@ -370,6 +379,9 @@ bool MatrixReader< Matrix >::readMatrixElementsFromMtxFile( std::istream& file, processedElements++; } } + + std::cout << "\nAfter while..." << std::endl; + file.clear(); long int fileSize = file.tellg(); timer.stop(); @@ -377,6 +389,9 @@ bool MatrixReader< Matrix >::readMatrixElementsFromMtxFile( std::istream& file, std::cout << " Reading the matrix elements ... " << processedElements << " / " << matrix.getNumberOfMatrixElements() << " -> " << timer.getRealTime() << " sec. i.e. " << fileSize / ( timer.getRealTime() * ( 1 << 20 )) << "MB/s." << std::endl; + + std::cout << "->END of reading matrix elements from file" << std::endl; + return true; }