쓰기 BOOST 희소 행렬

Question

친애 유래 커뮤니티 :

내가 밀어주는 아주 새로운 해요, 나는 그 uBLAS 라이브러리의 스파 스 매트릭스를 사용하는 것을 시도하고있다. 필자가 마주 치게되는 한 가지 문제는 스파 스 매트릭스를 계산 한 후이를 바이너리 파일에 저장하고 다른 프로그램에서 파일을 읽고 스파 스 매트릭스를 복구하는 것입니다. 일반적으로 나는 다음과 같은 방법으로 파일을 작성 (예를 들어 내가 100 개 수레와 배열 A를 경우) :

std::ofstream ofsDat("file.dat", std::ofstream::out); 
ofsDat.write((char*)A, sizeof(float)*100); 
ofsDat.close(); 

나는 BOOST 희소 행렬 비슷한 쓰기 작업을 할 수 있을지 궁금하네요을? 그렇다면 ofstream :: write의 두 번째 인수는 무엇이되어야 하는가? (데이터 청크의 크기 여야 함)

Answer 1

부스트 스파 스 매트릭스는 0이 아닌 값만 저장합니다. 행렬의 밀집한 표현에 대해서는 &A이 없습니다. 바이너리 출력이 필요한 경우, 직접 작성해야합니다. ios::out이 필요하지 않습니다. 당신은

std::ofstream test("./file.dat", std::ios::binary); 

또는 ofstream는/문자 스트림으로 처리하여 출력을 난도질 할 수 있어야합니다. 그 다음 :

#include <iostream> 
#include <fstream> 
#include <boost/numeric/ublas/matrix_sparse.hpp> 
#include <boost/numeric/ublas/storage.hpp> 

namespace ublas = boost::numeric::ublas; 

int main() 
{ 
    size_t width= 10; 
    size_t depth= 10; 
    ublas::compressed_matrix<double> m(width, depth); 
    m(0, 0) = 9; 
    m(1, 0) = 2; 
    m(0, 1) = 3; 
    m(5, 5) = 7; 

    std::ofstream test("./file.dat", std::ios::binary); 
    double zero= 0.0; 
    for(int i=0; i < width; ++i) 
     for (int j=0; j< depth; ++j) 
     { 
      double* temp= m.find_element(i, j); 
      if(temp) 
       test.write((char*)temp, sizeof(double)); 
      else 
       test.write((char*)&zero, sizeof(double)); 
     } 
    test.close(); 
} 

그러나 더 많은 정보를 파일에 쓰지 않으면 매트릭스는 알려진 너비와 깊이 여야합니다. 나는 실제로 상대방의 입력을 조사하여 텍스트 파일로 처리 할 수 ​​없는지 확인합니다. 그러면 다음을 수행하면됩니다.

test << m; 

Answer 2

마지막으로 몇 가지 검색 및 시도를 한 후에는이 작업을 수행하고 스파 스 매트릭스를 읽는 방법을 발견했습니다. 내 작업은 실제로 비교적 간단하므로 더 복잡하고 좀 더 일반적인 목적을 위해이 조잡한 방법이 효과적인지 아닌지를 알지 못합니다.

기본적인 아이디어는 부스트의 희소 행렬에있는 모든 0이 아닌 요소를 const_iterator (자세한 내용은 link 참조)를 통해 반복하여 ofstream에 쓰는 것입니다. 그리고 ifstream에서 읽으면, 나는 가난한 사람의 방법을 사용하고 있습니다 : 반복적으로 쓰기 형식에 따라 읽고 스파 스 매트릭스에 삽입하십시오. 여기에 내 코드 내 테스트하기위한 목적 :

#include <iostream> 
#include <fstream> 
#include <boost/numeric/ublas/matrix_sparse.hpp> 
#include <boost/numeric/ublas/io.hpp> 


    int main(int argc, char** argv) 
    { 
     using std::cerr; 
     using std::cout; 
     using std::endl; 
     using namespace boost::numeric::ublas; 
     typedef compressed_matrix<float, row_major> cMatrix; 

     const size_t size = 5; 
     const size_t rowInd[5] = { 0, 0, 1, 2, 4 }; 
     const size_t colInd[5] = { 0, 2, 0, 4, 4 }; 

     cMatrix sparseMat(size,size); 
     for (size_t i=0; i<size; ++i) 
      sparseMat.insert_element(rowInd[i], colInd[i], 1.0); 

     cout << sparseMat << endl; 

     // Try writing to file 
     std::ofstream ofsDat("temp.dat", std::ios::out | std::ios::binary); 
     for(cMatrix::const_iterator1 rowIter = sparseMat.begin1(); rowIter != sparseMat.end1(); ++rowIter) { 
      for(cMatrix::const_iterator2 colIter = rowIter.begin(); colIter != rowIter.end(); ++colIter) { 
       ofsDat << " " << colIter.index1() << " " << colIter.index2() << " " << *colIter; 
      }  // end for colIter 
     }  // end for rowIter 
     ofsDat.close(); 

     cout << "Writing ended, starting to read" << endl; 

     // Try reading the file 
     cMatrix sparseMat_2(size, size); 
     std::ifstream ifsDat("temp.dat", std::ios::in | std::ios::binary); 
     size_t rowTemp, colTemp; 
     float valTemp; 
     while(!ifsDat.eof()) { 
      ifsDat >> rowTemp >> colTemp >> valTemp; 
      cout << "row " << rowTemp << " column " << colTemp << " value " << valTemp << endl; 
      sparseMat_2.insert_element(rowTemp, colTemp, valTemp); 
     } 

     cout << sparseMat_2 << endl; 

     return 0; 
    } 

내가 구분 기호로 데이터 사이에 공백을 추가, 난이 할 수있는 더 좋은, 또는 표준 방법이 있는지 몰라? 모든 의견을 부탁드립니다!