Ticket #3547: compressed_mtx_ser.cpp

#include <boost/lexical_cast.hpp>
#include <boost/numeric/ublas/io.hpp>
#include <boost/archive/binary_iarchive.hpp>
#include <boost/archive/binary_oarchive.hpp>
#include <iostream>
#include <string>
#include <fstream>
#include "throw_with_info.hpp"
//#include "boost_extras/serialization/numeric/ublas/include/matrix_sparse.hpp"
#include <boost/numeric/ublas/matrix_sparse.hpp>


using namespace std;
using namespace boost;
using namespace numeric;
using namespace ublas;

#undef  __FUNCT__
#define __FUNCT__       "main"

int main(int argc,char *argv[]) {

        if (argc > 2 || argc < 1) {
                cerr << "\nSyntax:\n\t" << argv[0] << " <file name>\n" << endl;
                return  EXIT_FAILURE;
        }

        string          fname(lexical_cast<string>(argv[1]));

        compressed_matrix<double>                                       M(4,3);

        M(0,0) = 1.1;
        M(1,1) = 2.2;

        cout << "\nM =\n" << M << endl;


        ofstream        ofs(fname.c_str(),ios::out | ios::binary);

        {
                boost::archive::binary_oarchive oa(ofs);
                oa << M;
        }

        ofs.close();

        M.clear();

        compressed_matrix<double>                                       M2;

        ifstream        ifs(fname.c_str(),ios::in | ios::binary);

        {
                boost::archive::binary_iarchive ia(ifs);
                ia >> M2;
        }

        ifs.close();

        cout << "\nre-read matrix:\n" << M2 << endl;

        cout << "\nDone.\n" << endl;

        return EXIT_SUCCESS;
}