---

more/num/matrixio.h

Go to the documentation of this file.

//  Copyright (C) 1998--2002  Petter Urkedal
//
//  This file is free software; you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation; either version 2 of the License, or
//  (at your option) any later version.
//
//  This file is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.
//
//  You should have received a copy of the GNU General Public License
//  along with this program; if not, write to the Free Software
//  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
//  As a special exception, you may use this file as part of a free
//  software library without restriction.  Specifically, if other files
//  instantiate templates or use macros or inline functions from this
//  file, or you compile this file and link it with other files to
//  produce an executable, this file does not by itself cause the
//  resulting executable to be covered by the GNU General Public
//  License.  This exception does not however invalidate any other
//  reasons why the executable file might be covered by the GNU General
//  Public License.

//  $Id: matrixio.h,v 1.2 2002/09/01 12:02:10 petter_urkedal Exp $


#ifndef MORE_NUM_MATRIXIO_H
#define MORE_NUM_MATRIXIO_H

///\if unmaintained
#include <more/num/matrix.h>
#include <iostream>
#include <iomanip>

namespace more {
namespace num {

  template <class T, class Sg>
    std::ostream& operator<<(std::ostream& out, const matrix<T, Sg>& x) {
        std::ostream::fmtflags flags = out.flags();
        std::streamsize prec = out.precision();
        std::streamsize width = out.width();
        out.width(0);
        out << x.size(0) << ' ' << x.size(1) << '\n';
        for(int i = 0; i < x.size(0); i++) {
            for(int j = 0; j < x.size(1); j++) {
                out.flags(flags);
                out.precision(prec);
                out.width(width);
                out << x(i, j) << ' ';
                if(out.fail()) return out;
            }
            out << '\n';
        }
        return out;
    }

  template <class T, class Sg>
    std::istream& operator>>(std::istream& in, matrix<T, Sg> x) {
        int n, m;
        in >> n >> m;
        matrix<T, Sg> y(n, m);
        for(int i = 0; i < x.size(0); i++)
            for(int j = 0; j < x.size(1); j++) {
                in >> x(i, j);
                if(in.fail()) return in;
            }
        x = y;
        return in;
    }

}}
///\endif

#endif

---