#include <iostream.h>
#include "Matrix3.h"

Matrix3::Matrix3()
{  int i, j;

   for (i = 0; i < 3; i=i+1)
      for (j = 0; j < 3; j=j+1)
         A[i][j] = i + j;
}

/************************************************
  Matrix vector multiply member FRIEND function
 *************************************************/
Vector3 Matrix3::operator*(Vector3 & v)
{
   Vector3 z;
   int i, j;

   for (i = 0; i < 3; i=i+1)
      z.x[i] = 0;

   for (i = 0; i < 3; i=i+1)
      for (j = 0; j < 3; j=j+1)
         z.x[i] = z.x[i] + A[i][j]*v.x[j];

   return z;
}


/* -------------------------------------------------------
   Friend function to print a matrix
   ------------------------------------------------------- */
ostream & operator<<(ostream & cout, Matrix3 & M)
{
   int i, j;

   for (i = 0; i < 3; i=i+1)
   {
      for (j = 0; j < 3; j=j+1)
         cout << M.A[i][j] << "\t";
      cout << endl;
   }

   return(cout);
}