#ifndef MATIO_H
#define MATIO_H

#include<iostream>

namespace Eigen{
    template<class Matrix>
    void write_stream(std::ostream &ostr, const Matrix& matrix){
        typename Matrix::Index rows=matrix.rows(), cols=matrix.cols();
        ostr<<rows<<" "<<cols<<std::endl;
        for (int r=0;r<rows;++r)
        {
            for (int c=0;c<cols;++c)
                ostr<<matrix(r,c)<<" ";
            ostr<<std::endl;
        }
    }
    template<class Matrix>
    void read_stream(std::istream &istr, Matrix& matrix){
        typename Matrix::Index rows=0, cols=0;
        istr>>rows>>cols;
        matrix.resize(rows, cols);
        for (int r=0;r<rows;++r)
            for (int c=0;c<cols;++c)
                istr>>matrix(r,c);
    }
} // Eigen::


#endif