matrix.cpp

/***************************************************************************
 *   Copyright (C) 2005 by Nestor Aguirre                                  *
 *   nfaguirrec@unal.edu.co                                                *
 *                                                                         *
 *   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.                                   *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU General Public License     *
 *   along with this program; if not, write to the                         *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/

// C++
#include <fstream>
using namespace std;

// Project
#include <SawMatrix.h>
#include <SawVector.h>
#include <SawVectorList.h>
#include <SawUtils.h>

int main(int argc, char *argv[]){
        
        cout << "****************************************************" << endl ;
        cout << "*                                                  *" << endl ;
        cout << "*               MATRIX MANIPULATION                *" << endl ;
        cout << "*                                                  *" << endl ;
        cout << "*               Nestor Aguirre  2006               *" << endl ;
        cout << "*              nfaguirrec@unal.edu.co              *" << endl ;
        cout << "*                                                  *" << endl ;
        cout << "****************************************************" << endl ;
        cout << endl ;
        
        cout << "---------------------------------------" << endl ;
        cout << "   Initializing to constant value      " << endl ;
        cout << "---------------------------------------" << endl ;
        SawMatrix matrix1(4, 4, 2.36) ;
        cout << "A = " << matrix1 << endl ;
        cout << endl ;
        
        cout << "-------------------------------" << endl ;
        cout << "   Getting data of C array     " << endl ;
        cout << "-------------------------------" << endl ;
        double mat[] = { 2.350, 0.156, 4.521, 2.534,
                         0.789, 5.423, 6.156, 5.123,
                         5.489, 6.464, 1.453, 6.156,
                         5.156, 2.456, 6.156, 9.156 } ;
        SawMatrix matrix2( mat, 4, 4 ) ;
        cout << "B = " << matrix2 << endl ;
        cout << endl ;
                
        cout << "---------------------------------------" << endl ;
        cout << "   Probando el constructor de copia    " << endl ;
        cout << "---------------------------------------" << endl ;
        SawMatrix matrix3( matrix2 ) ;
        cout << "C = "<< matrix3 << endl ;
        cout << endl ;
        
        cout << "----------------------------" << endl ;
        cout << "   Probando el operador ==  " << endl ;
        cout << "----------------------------" << endl ;
        cout << "C == B ? " << ( (matrix3 == matrix2) ? "TRUE" : "FALSE" ) << endl ;
        cout << "A == B ? " << ( (matrix1 == matrix2) ? "TRUE" : "FALSE" ) << endl ;
        cout << endl ;
        
        cout << "-------------------------------------------" << endl ;
        cout << "   Obteniendo el elemento maximo y minimo  " << endl ;
        cout << "-------------------------------------------" << endl ;
        uint i, j ;
        double max = matrix3.getMax( &i, &j ) ;
        cout << "MAX(C):   vec[" << i << "][" << j << "] = " << max << endl ;
        double min = matrix3.getMin( &i, &j ) ;
        cout << "MIN(C):   vec[" << i << "][" << j << "] = " << min << endl ;
        cout << endl ;
        
        cout << "------------------------------------------" << endl ;
        cout << "   Probando el operador de asignacion =   " << endl ;
        cout << "------------------------------------------" << endl ;
        SawMatrix matrix4( 4, 4 ) ;
        matrix4 = matrix3 ;
        cout << "D = " << matrix4 << endl ;
        cout << endl ;
        
        cout << "--------------------------------------------" << endl ;
        cout << "   Probando intercambiar columnas 1 <--> 3  " << endl ;
        cout << "--------------------------------------------" << endl ;
        matrix4.swapColumns( 1, 3 ) ;
        cout << "D_{1c,3c} = " << matrix4 << endl ;
        cout << endl ;
        
        cout << "-----------------------------------------" << endl ;
        cout << "   Probando intercambiar filas 1 <--> 3  " << endl ;
        cout << "-----------------------------------------" << endl ;
        matrix4.swapRows( 1, 3 ) ;
        cout << "D_{1r,3r} = " << matrix4 << endl ;
        cout << endl ;
        
        cout << "-------------------------------" << endl ;
        cout << "  Probando sumar dos matrices  " << endl ;
        cout << "-------------------------------" << endl ;
        cout << "A   = " << matrix3 << endl ;
        cout << "B   = " << matrix4 << endl ;
        cout << "A+B = " << matrix3 + matrix4 << endl ;
        cout << endl ;
        
        cout << "--------------------------------" << endl ;
        cout << "  Probando restar dos matrices  " << endl ;
        cout << "--------------------------------" << endl ;
        cout << "A   = " << matrix3 << endl ;
        cout << "B   = " << matrix4 << endl ;
        cout << "A-B = " << matrix3 - matrix4 << endl ;
        cout << endl ;
        
        cout << "--------------------------------------------------" << endl ;
        cout << "  Probando multiplicar un escalar con una matriz  " << endl ;
        cout << "--------------------------------------------------" << endl ;
        cout << "A   = " << matrix4 << endl ;
        cout << "B   = " << 0.456 << endl << endl ;
        cout << "A*B = " << matrix4*0.456 << endl ;
        cout << "B*A = " << 0.456*matrix4 << endl ;
        cout << endl ;
        
        cout << "------------------------------------" << endl ;
        cout << "  Obteniendo la matriz transpuesta  " << endl ;
        cout << "------------------------------------" << endl ;
        cout << "A^t   = " << matrix4.getTransPose() << endl ;
        cout << endl ;
        
        cout << "-----------------------------------------" << endl ;
        cout << "   Obteniendo los datos desde un stream  " << endl ;
        cout << "-----------------------------------------" << endl ;
        SawMatrix matrix5( 4, 4 ) ;
        ifstream file("matrix.data") ;
        file >> matrix5 ;
        file.close() ;
        cout << matrix5 << endl ;
        
        cout << "-------------------------------" << endl ;
        cout << "   Multiplicando dos matrices  " << endl ;
        cout << "-------------------------------" << endl ;
        
        double dataA[] = { 0.116, 0.126, 0.138,
                           0.213, 0.223, 0.235,
                           0.566, 6.456, 8.321,
                           6.546, 9.156, 0.456 } ;
                       
        double dataB[] = { 1.011, 7.012,
                           6.216, 8.022,
                           2.031, 5.032 } ;

        SawMatrix A( dataA, 4, 3 ) ;
        SawMatrix B( dataB, 3, 2 ) ;

        cout << "A = " << A << endl ;
        cout << "B = " << B << endl ;
        cout << "A*B = " << A*B << endl ;
        
        cout << "------------------------------------------------------------" << endl ;
        cout << "   Obteniendo los valores y vectores propios de una matriz  " << endl ;
        cout << "------------------------------------------------------------" << endl ;

        SawMatrix matrix6( mat, 4, 4 ) ;
        ifstream file2("matrix.data") ;
        file2 >> matrix6 ;
        file2.close() ;
        
        SawVector eigenValues(4) ;
        SawVectorList eigenVectors(4, 4) ;
        eigenValues = matrix6.getEigenValues( &eigenVectors ) ;
        
        cout << "A = " << matrix6 << endl ;
        
        cout << "  Valores propios  " << endl ;
        cout << "-------------------" << endl ;
        SawUtils::writeWithFormatMapStyle( cout, eigenValues, eigenVectors, 10 ) ;
        cout << endl ;
        
        cout << "---------------------------------------------------" << endl ;
        cout << "   Calculando la matriz inversa y el determinante  " << endl ;
        cout << "---------------------------------------------------" << endl ;
        cout << "A = " << matrix6 << endl ;
        SawMatrix invers(4,4) ;
        invers = matrix6.getInverse() ;
        cout << "A^(-1) = " << invers << endl ;
        cout << "A*A^(-1) = " << matrix6*invers << endl ;
        cout << "Det(A) = " << matrix6.getDeterminant() << endl ;
        
        return EXIT_SUCCESS;
}

---