Vector.h

Go to the documentation of this file.

/***************************************************************************
                          Vector.h  -  description
                             -------------------
    begin                : Thu Jun 7 2001
    copyright            : (C) 2001 by Olivier Roussel & Alexei Tsigulin
    email                : roussel@ict.uni-karlsruhe.de, lpsoft@mail.ru
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   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.                                   *
 *                                                                         *
 ***************************************************************************/


class Vector
{

/*
______________________________________________________________________________________________

    Constructor and distructor
______________________________________________________________________________________________

*/


public:
  Vector();

Vector(const int n);

Vector(const real x, const real y);

Vector(const real x, const real y, const real z);

Vector(const Vector& V);

~Vector();   
/*
______________________________________________________________________________________________

    Set and get
______________________________________________________________________________________________

*/


inline void setValue(const int n, const real a);

void setZero();

void setDimension(const int n);

inline real value(const int n) const;

inline int dimension() const;


/*
______________________________________________________________________________________________

    Operators
______________________________________________________________________________________________

*/

bool operator== (const Vector& V ) const;

void operator= (const Vector& V );

void operator+=(const Vector& V );

Vector operator+ (const Vector& V ) const;

void operator-=(const Vector& V );

Vector operator- (const Vector& V ) const;

Vector operator- () const;

void operator*=(const real a );

Vector operator* (const real a ) const;

void operator/=(const real a );

Vector operator/ (const real a ) const;

real operator* (const Vector& V ) const;
    
Vector operator| (const Vector& V) const;

Vector operator^ (const Vector& V) const;

bool isNaN() const;

/*
______________________________________________________________________________________________

PRIVATE VARIABLES
______________________________________________________________________________________________

*/
//parallel modification
//private:
public:
    int  Columns;   
    /* real *U; // Components */
    real U[9];  
};
/*
______________________________________________________________________________________________

EXTERNAL FUNCTIONS
______________________________________________________________________________________________

*/

Vector operator* (const real a, const Vector& V);

Vector abs (const Vector& V);

int dim (const Vector& V);

real N1(const Vector& V);

real N2(const Vector& V);

real NMax(const Vector& V);

ostream& operator<<(ostream& out, const Vector& V);

/*
______________________________________________________________________________________________

INLINE FUNCTIONS
______________________________________________________________________________________________

*/
inline int Vector::dimension() const
{
    return Columns;
}

/*
______________________________________________________________________________________________

*/

inline void Vector::setValue(const int n, const real a)
{

#ifdef DEBUG
    if ( n <= 0 || n > Columns)
    {
        cout << "Vector.cpp: In method `void Vector::setValue(int, real)':\n";
        cout << "Vector.cpp: first argument out of range\n";
        cout << "carmen: *** [Vector.o] Execution error\n";
        cout << "carmen: abort execution.\n";
        exit(1);
    }
#endif

        *(U+n-1) = a;
}
/*
______________________________________________________________________________________________

*/
inline real Vector::value(const int n) const
{

#ifdef DEBUG

    if ( n <= 0 || n > Columns)
    {
        cout << "Vector.cpp: In method `void Vector::value(int)':\n";
        cout << "Vector.cpp: argument out of range\n";
        cout << "carmen: *** [Vector.o] Execution error\n";
        cout << "carmen: abort execution.\n";
        exit(1);
    }
#endif

    return *(U+n-1);
}