SawMatrix.h

/***************************************************************************
 *   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.             *
 ***************************************************************************/
#ifndef GPPMATRIX_H
#define GPPMATRIX_H

/***********************************
* AnSII C++
*/
#include <math.h>

//**********************************

/***********************************
* gsl Library
*/
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_blas.h>
#include <gsl/gsl_eigen.h>
#include <gsl/gsl_linalg.h>

//**********************************

/***********************************
* Project
*/
#include <SawObject.h>

class SawVector ;
class SawVectorList ;

//**********************************

class SawMatrix : public SawObject
{
        public:
                SawMatrix( uint rows , uint cols , double a = 0.0 );
                SawMatrix( double* array, uint rows, uint cols );
                SawMatrix( const SawMatrix& m );
                void operator = ( const SawMatrix& m );
                ~SawMatrix();
                
                void swapRows( uint i, uint j ) ;
                void swapColumns( uint i, uint j ) ;
                
                gsl_matrix* getGSLMatrix() const ;
                void setGSLMatrix( gsl_matrix* matrix ) ;
                uint getRows() const ;
                uint getColumns() const ;
                double get( uint i, uint j ) const;
                void set( uint i, uint j, double value );
                
                void setIdentity() ;
                void setNull() ;
                
                bool operator == ( SawMatrix m );
                bool operator != ( SawMatrix m );
                
                double getMax( uint* i=NULL, uint *j=NULL ) const ;
                double getMin( uint* i=NULL, uint *j=NULL ) const ;
                
                double getDeterminant() ;
                SawMatrix getInverse() ;
                SawVector getEigenValues( SawVectorList* eigenVectors = NULL ) ;
                
                bool isNull() const ;
                
                SawMatrix getTransPose() ;
                
                friend SawMatrix operator + ( SawMatrix a, SawMatrix b ) ;
                friend SawMatrix operator - ( SawMatrix a, SawMatrix b ) ;
                friend SawMatrix operator * ( SawMatrix a, SawMatrix b ) ;
                
                friend SawMatrix operator * ( SawMatrix m, double a ) ;
                friend SawMatrix operator * ( double a, SawMatrix m ) ;
                
                friend ostream& operator << ( ostream& os, const SawMatrix& m ) ;
                friend istream& operator >> ( istream& is, const SawMatrix& m ) ;
                
                static SawMatrix genericFunction( SawMatrix a, double(* func)(double) ) ;
                static SawMatrix genericFunction( SawMatrix a, double(* func)(double, double), double param = 0.0 ) ;
                static SawMatrix sin( SawMatrix a ) ;
                static SawMatrix cos( SawMatrix a ) ;
                static SawMatrix tan( SawMatrix a ) ;
                static SawMatrix asin( SawMatrix a ) ;
                static SawMatrix acos( SawMatrix a ) ;
                static SawMatrix atan( SawMatrix a ) ;
                static SawMatrix sinh( SawMatrix a ) ;
                static SawMatrix cosh( SawMatrix a ) ;
                static SawMatrix exp( SawMatrix a ) ;
                static SawMatrix ln( SawMatrix a ) ;
                static SawMatrix log( SawMatrix a ) ;
                static SawMatrix sqrt( SawMatrix a ) ;
                static SawMatrix pow( SawMatrix a, double e ) ;
                
        private:
                static SawMatrix genericFunc( SawMatrix a, double(* func)(double) ) ;
                
        private:
                gsl_matrix* matrix ;
                uint rows ;
                uint cols ;
};

#endif

---