/* File: matrix.H         Copyright glh 10/21/92 */
#ifndef MATRIXH
#define MATRIXH
#include <math.h>
#include <stdlib.h>
#include <iostream.h>
#include <assert.h>

typedef enum {false, true} Boolean;

template<class T>
class Matrix {
  private:
    int  row_dim, col_dim;
    T**  element;
  public:
    inline Matrix();                              // default constructor
    Matrix(const Matrix<T>&);                     // copy constructor
    Matrix(int rdim, int cdim);                   // constructor
    Matrix(int rdim, int cdim, const T& initval); // constructor
    Matrix(int rdim, int cdim, T* initval);       // constructor
    ~Matrix();                                    // destructor
    inline int Row() const;
    inline int Col() const;
    inline T Retrieve(int row, int col) const;
    inline void Assign(int row, int col, const T& val);
    Matrix<T>& operator=(const Matrix<T>&);
    friend ostream& operator<<(ostream& os, const Matrix<T>&); 
    friend Matrix<T>& operator+(const Matrix<T>&, const Matrix<T>&);
    friend Matrix<T>& operator-(const Matrix<T>&, const Matrix<T>&);
    friend Matrix<T>& operator-(const Matrix<T>&);              
    friend Matrix<T>& operator*(const Matrix<T>&, const Matrix<T>&);
    friend Matrix<T>& operator*(const T& scalar, const Matrix<T>&);
    friend Matrix<T>& Transpose(const Matrix<T>&);
    // remaining member functions not shown
};

template<class T>
inline Matrix<T>::Matrix() {row_dim=col_dim=0; element=NULL;}

template<class T>
inline int Matrix<T>::Row() const {return row_dim;}

template<class T>
inline int Matrix<T>::Col() const {return col_dim;}

template<class T>
inline T Matrix<T>::Retrieve(int row, int col) const 
{return element[row][col];}

template<class T>
inline void Matrix<T>::Assign(int row, int col, const T& val) 
{element[row][col] = val;}

#endif