matrixFunctions.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>
#include <math.h>
using namespace std;

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

double myFunction1( double x ){ return pow( x, 0.5 ) ; }
double myFunction2( double x ){ return sin( x ) ; }

int main(int argc, char *argv[]){
        
        cout << "************************************************" << endl ;
        cout << "*                                              *" << endl ;
        cout << "*            MATRIX FUNCTIONS TEST             *" << endl ;
        cout << "*             Nestor Aguirre  2006             *" << endl ;
        cout << "*            nfaguirrec@unal.edu.co            *" << endl ;
        cout << "*                                              *" << endl ;
        cout << "************************************************" << endl ;
        cout << endl ;
        
        /**********************************************************
        * Se inicializa la matriz cargandola desde un stream
        */
        SawMatrix matrix( 7, 7 ) ;
        ifstream file("overlapH2O.data") ;
        file >> matrix ;
        file.close() ;
        
        cout << endl ;
        cout << "-----------------------------------------" << endl ;
        cout << "   Usando funciones predefinidas en Saw  " << endl ;
        cout << "-----------------------------------------" << endl ;
        cout << "A = " << matrix << endl ;
        cout << "A^(0.5) = " << SawMatrix::pow( matrix, -0.5 ) << endl ;
        cout << "sin(A) = " << SawMatrix::sin( matrix ) << endl ;
        
        cout << "-----------------------------------------" << endl ;
        cout << "   Usando SawMatrix::genericFunction()   " << endl ;
        cout << "-----------------------------------------" << endl ;
        cout << endl ;
        cout << "A^(-0.5) = " << SawMatrix::genericFunction( matrix, *myFunction1 ) << endl ;
        cout << "sin(A) = " << SawMatrix::genericFunction( matrix, *myFunction2 ) << endl ;
        
        return EXIT_SUCCESS;
}

---