/***************************************************************************
                          matrix.h  -  description
                             -------------------
    begin                : Fri Jan 17 19:41:58 CST 2003
    copyright            : (C) 2003 by Daniel Eric Smith
    email                : 
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program 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.                                   *
 *                                                                         *
 ***************************************************************************/


#include <math.h>
#include <iomanip.h>

const int SIZE=20;

class matrix{
  int nrow,ncol;
  float elt[SIZE][SIZE];
 public:
  //matrix(); //{nrow=1; ncol=1;}
  matrix(int r, int c) {nrow=r; ncol=c;}  //parameterized construtor
  matrix trans();
/*
ARITIMETIC
*/
  matrix operator+(matrix N); //rational addition r+s
  matrix operator-(matrix N);
  matrix operator*(matrix N);
  //matrix operator/(matrix N);
/*
INPUT/OUTPUT
*/
  friend istream& operator>>(istream& s, matrix& N);
  friend ostream& operator<<(ostream& s, matrix N);
};

matrix matrix::trans()
{
  matrix t(ncol,nrow);

  for(int i=0;i<nrow;i++)
    {
      for(int j=0;j<ncol;j++)
	{
	  t.elt[j][i]=elt[i][j];
	}
    }
  return t;
}
/************
* ARITMETIC *
************/

/****************
*               *
* r+s  ADDITION * 
*               *
****************/
matrix matrix::operator+(matrix N)
{
  matrix t(N.nrow,N.ncol);

  for(int j=0;j<nrow;j++)
    {
      for(int i=0;i<ncol;i++)
	{
	  t.elt[j][i]=elt[j][i]+N.elt[j][i];
	}
    }

  return t;	
}
/******************
*                 *
* r-s SUBTRACTION *
*                 *
******************/
matrix matrix::operator-(matrix N)
{
  matrix t(N.nrow,N.ncol);

  for(int j=0;j<nrow;j++)
    {
      for(int i=0;i<ncol;i++)
	{
	  t.elt[j][i]=elt[j][i]-N.elt[j][i];
	}
    }

  return t;		
}
/*********************
*                    *
* r*s MULTIPLICATION *
*                    *
*********************/
matrix matrix::operator*(matrix N)
{
  matrix t(nrow,N.ncol);
  float s=0;
  
  for(int j=0;j<nrow;j++)
    {
      for(int i=0;i<N.ncol;i++)
	{
	  s=0;
	  for(int k=0;k<ncol;k++)
	    {
	      s+=elt[j][k]*N.elt[k][i];
	    }
	  t.elt[j][i]=s;	  
	}
    }

  return t;
}
/***************
*              *
* r/s DIVISION *
*              *
***************/
/*matrix matrix::operator/(matrix N)
{
  rational t;

  t.num=num*r.den;
  t.den=den*r.num;

  t.reduce();

  return t;
}*/

/***************
* INPUT/OUTPUT *
***************/
istream& operator>>(istream& s, matrix& N)
{
  for(int j=0;j<N.nrow;j++)
    {
      cout <<"Enter the "<<j<<" row:"<<endl;
      for(int i=0;i<N.ncol;i++)
	{
	  s>>N.elt[j][i];
	}
    }
  return s;
}

ostream& operator<<(ostream& s, matrix N)
{
  for(int j=0;j<N.nrow;j++)
    {
      s<<endl<<"[";
      for(int i=0;i<N.ncol;i++)
	{
	  s<<N.elt[j][i]<<" ";
	}
      s<<"]";
    }

  return s;
}