source: git/factory/templates/ftmpl_matrix.h @ c5f2101

// emacs edit mode for this file is -*- C++ -*-
// $Id: ftmpl_matrix.h,v 1.1 1996-12-18 15:04:42 schmidt Exp $

#ifndef INCL_MATRIX_H
#define INCL_MATRIX_H

/*
 * $Log: not supported by cvs2svn $
 */

#include <iostream.h>

template <class T>
class SubMatrix;

template <class T>
class Matrix
{
private:
    int NR, NC;
    T ** elems;
    void printrow ( ostream & s, int i ) const;
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 );
    void print( ostream& s ) const;
    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 );
    friend ostream & operator<< ( ostream & s, const Matrix<T>& M )
    {
        M.print( s );
        return s;
    }
    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 */