SawVector.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 GPPVECTOR_H
#define GPPVECTOR_H

/***********************************
* gsl Library
*/
#include <gsl/gsl_vector.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_blas.h>

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

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

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

class SawVector : public SawObject
{
        public:
                SawVector( uint size = 1, double a = 0.0 );
                SawVector( double* array, uint size );
                SawVector( const SawVector& vector );
                void operator = ( const SawVector& vector );
                ~SawVector();
                
                void swapElements( uint i, uint j ) ;
                void reverse() ;
                
                gsl_vector* getGSLVector() const ;
                void setGSLVector( gsl_vector* vector ) ;
                uint getSize() const ;
                double get( uint i ) const;
                void set( uint i, double value );
                
                bool operator == ( SawVector vector );
                bool operator != ( SawVector vector );
                
                double getMax( uint* i=NULL ) const ;
                double getMin( uint* i=NULL ) const ;
                
                bool isNull() const ;
                
                SawVector operator + ( SawVector vector );
                SawVector operator - ( SawVector vector );
                double operator * ( SawVector vector );
                double getEuclideanNorm() ;
                
                friend SawVector operator * ( SawVector vector, double a );
                friend SawVector operator * ( double a, SawVector vector );
                friend SawVector operator + ( SawVector vector, double a );
                friend SawVector operator + ( double a, SawVector vector );
                friend SawVector operator ^ ( double a, SawVector vector );
                
                friend ostream& operator << ( ostream& os, const SawVector& e ) ;
                friend istream& operator >> ( istream& is, const SawVector& e ) ;
                
        private:
                gsl_vector* vector ;
                uint size ;
};

#endif

---