source: git/factory/templates/ftmpl_matrix.h @ 6881f9b

// emacs edit mode for this file is -*- C++ -*-
// $Id: ftmpl_matrix.h,v 1.5 1997-04-30 13:08:28 schmidt Exp $

#ifndef INCL_MATRIX_H
#define INCL_MATRIX_H

/*
 * $Log: not supported by cvs2svn $
 * Revision 1.4  1997/04/18 16:42:01  schmidt
 * o class Matrix: typedef T* T_ptr; added
 *
 * Revision 1.3  1997/04/15 10:13:42  schmidt
 * #include <config.h> added
 * the header config.h will be included be makeheader
 *
 * Revision 1.2  1997/03/27 10:34:15  schmidt
 * stream-io wrapped by NOSTREAMIO
 *
 * Revision 1.1  1996/12/18 15:04:42  schmidt
 * Initial revision
 *
 */

#include <factoryconf.h>

#ifndef NOSTREAMIO
#include <iostream.h>
#endif /* NOSTREAMIO */

template <class T>
class SubMatrix;

template <class T>
class Matrix
{
private:
    int NR, NC;
    T ** elems;
#ifndef NOSTREAMIO
    void printrow ( ostream & s, int i ) const;
#endif /* NOSTREAMIO */
    typedef T* T_ptr;
public:
    Matrix() : NR(0), NC(0), elems(0) {}
    Matrix( int nr, int nc );
    Matrix( const Matrix<T>& M );
    ~Matrix();
    Matrix<T>& operator= ( const Matrix<T>& M );
    int rows() const { return NR; }
    int columns() const { return NC; }
    SubMatrix<T> operator[] ( int i );
    const SubMatrix<T> operator[] ( int i ) const;
    T& operator() ( int row, int col );
    T operator() ( int row, int col ) const;
    SubMatrix<T> operator() ( int rmin, int rmax, int cmin, int cmax );
    const SubMatrix<T> operator() ( int rmin, int rmax, int cmin, int cmax ) const;
    void swapRow( int i, int j );
    void swapColumn( int i, int j );
    friend Matrix<T> operator+ ( const Matrix<T>& lhs, const Matrix<T>& rhs );
    friend Matrix<T> operator- ( const Matrix<T>& lhs, const Matrix<T>& rhs );
    friend Matrix<T> operator* ( const Matrix<T>& lhs, const Matrix<T>& rhs );
    friend Matrix<T> operator* ( const Matrix<T>& lhs, const T& rhs );
    friend Matrix<T> operator* ( const T& lhs, const Matrix<T>& rhs );
#ifndef NOSTREAMIO
    void print( ostream& s ) const;
    friend ostream & operator<< ( ostream & s, const Matrix<T>& M )
    {
        M.print( s );
        return s;
    }
#endif /* NOSTREAMIO */
    friend class SubMatrix<T>;
};

template <class T>
class SubMatrix
{
private:
    int r_min, r_max, c_min, c_max;
    Matrix<T>& M;
    // we do not provide a default ctor, so nobody can declare an empty SubMatrix
    SubMatrix( int rmin, int rmax, int cmin, int cmax, const Matrix<T> & m );
public:
    SubMatrix( const SubMatrix<T> & S );
    SubMatrix<T>& operator= ( const SubMatrix<T>& S );
    SubMatrix<T>& operator= ( const Matrix<T>& S );
    operator Matrix<T>() const;
    T operator[] ( int i ) const;
    T& operator[] ( int i );
    friend class Matrix<T>;
};




#endif /* INCL_MATRIX_H */