MatrixMultiplication_main.cpp

/**
\file     MatrixMultipliation_main.cpp
\brief    An program using the Zenautics Matrix Class to demonstrate matrix multiplication.
\author   Glenn D. MacGougan (GDM)
\date     2008-04-25

\b License \b Information \n
Copyright (c) 2008, Glenn D. MacGougan, Zenautics Technologies Inc. \n
This file may be used, copied, modified and redistributed without restriction.
*/

#include <iostream>
#include <string>
#include <complex>
#include "Matrix.h"
#include <time.h>

#if defined _MSC_VER && _MSC_VER < 1400
#define _CRT_SECURE_NO_DEPRECATE
#endif

#ifndef _MSC_VER
#define _CRT_SECURE_NO_DEPRECATE
#endif


using namespace Zenautics;
using namespace std;
 
int main( int argc, char* argv[] )                            
{  
  try
  {
    // e.g.
    Matrix A;
    Matrix B;
    Matrix C;
    char c;
    int i;
    int nr;
    char linestr[256];
    char MatrixBuffer[65536];    
    int scount = 0;
    int nrTimes;
     time_t  t0, t1; // time_t is defined on <time.h> and <sys/types.h> as long
    clock_t c0, c1; // clock_t is defined on <time.h> and <sys/types.h> as int 


    cout << "A = B*C" << endl;
    cout << "(Press U): User input matrix B and matrix C." << endl;
    cout << "(Press S): B = 1:1:N         Column Matrix" << endl;
    cout << "           C = 0.1:0.1:N/10  Column Matrix" << endl;
    cout << "           B = B*B.Transpose()" << endl;
    cout << "           C = C*C.Transpose()" << endl;
    cout << "Enter Choice:";
    cin.get( c );
    if( cin.fail() )
    {
      cout << "Bad input. Aborting." << endl;
      return 1;
    }
    c = toupper( c );

    cout << "Enter the number of times to perform the multiplication:";
    cin >> nrTimes;
    if( nrTimes <= 0 )
    {
      cout << "Bad input. Aborting." << endl;
      return 1;
    }

    if( c == 'U' )
    {
      cout << "You will enter Matrix B and Matrix C." << endl;
      cout << "Ensure that the matrices are conformal for multiplication." << endl;

      cout << "Enter the number of rows B:";
      cin >> nr;
      if( nr < 0 )
      {
        cout << "Bad input. Aborting." << endl;
        return 1;
      }
      cout << "\nThe number of columns is determined automatically." << endl;
      cout << "Enter each row of the matrix then press Enter." << endl;
      cin.clear();
      for( i = 0; i < nr; i++ )
      {
        cin.getline( linestr, 256 );
        if( cin.fail() )
        {
          cout << "Bad input. Aborting." << endl;
          return 1;
        }
        if( strlen( linestr ) == 0 )
        {
          // cin is weird
          i--;
          continue;
        }
#ifndef _CRT_SECURE_NO_DEPRECATE
        scount += sprintf_s( MatrixBuffer+scount, 65536-scount, "%s\r\n", linestr );
#else
        scount += sprintf( MatrixBuffer+scount, "%s\r\n", linestr );
#endif
      }
      if( !B.SetFromMatrixString( MatrixBuffer ) )
      {
        cout << "Bad input. Aborting." << endl;
        return 1;
      }

      cout << "Enter the number of rows C:";
      cin >> nr;
      if( nr < 0 )
      {
        cout << "Bad input. Aborting." << endl;
        return 1;
      }
      cout << "\nThe number of columns is determined automatically." << endl;
      cout << "Enter each row of the matrix then press Enter." << endl;
      scount = 0;
      cin.clear();
      for( i = 0; i < nr; i++ )
      {
        cin.getline( linestr, 256 );
        if( cin.fail() )
        {
          cout << "Bad input. Aborting." << endl;
          return 1;
        }
        if( strlen( linestr ) == 0 )
        {
          // cin is weird
          i--;
          continue;
        }
#ifndef _CRT_SECURE_NO_DEPRECATE
        scount += sprintf_s( MatrixBuffer+scount, 65536-scount, "%s\n", linestr );
#else
        scount += sprintf( MatrixBuffer+scount, "%s\n", linestr );
#endif
      }
      if( !C.SetFromMatrixString( MatrixBuffer ) )
      {
        cout << "Bad input. Aborting." << endl;
        return 1;
      }

      cout << "B = " << endl;
      B.PrintStdout();
      cout << "C = " << endl;
      C.PrintStdout();
    }
    else
    {
      cout << "Enter the square dimensions of Matrix B and C: " << endl;
      cin >> nr;
      if( nr < 0 )
      {
        cout << "Bad input. Aborting." << endl;
        return 1;
      }
      if( !B.Inplace_colon( 1.0, 1.0, nr ) )
      {
        cout << "B.Inplace_colon returned false. Bad input. Aborting." << endl;
        return 1;
      }
      if( !C.Inplace_colon( 0.1, 0.1, nr/10.0+1e-09 ) )
      {
        cout << "B.Inplace_colon returned false. Bad input. Aborting." << endl;
        return 1;
      }
      B = B * B.Transpose();
      C = C * C.Transpose();

      cout << "A = B*C;" << endl;
      if( B.nrows() < 20 )
      {
        cout << "B = " << endl;
        B.PrintStdout();
      }
      if( C.nrows() < 20 )
      {
        cout << "C = " << endl;
        C.PrintStdout();
      }

     
    }
    
    t0 = time(NULL); 
    c0 = clock();

    for( i = 0; i < nrTimes; i++ )
    {
      A = B*C;
      if( i == 0 )
      {
        cout << "A = B*C;" << endl;
        cout << "A = " << endl;
        A.PrintStdout();
      }
    }

    t1 = time(NULL); 
    c1 = clock(); 
    printf( "Multipliciation testing took %10.1f ms\n", (float) (c1 - c0)/CLOCKS_PER_SEC*1000.0 );     
  }                                                           
  catch( MatrixException& matrix_exception )
  {
    cout << matrix_exception << endl;
  }
  catch ( ... )
  {
    cout << "Caught unknown exception" << endl;
  }

  char anykey;
  cout << "\nProgram complete.\n" << endl;
  cout << "Type something and hit enter: " << endl;
  cin.get( anykey );
  cin.get( anykey );
  return 0;
}                                   

---