Matrix.h

Go to the documentation of this file.

/***************************************************************************
                          Matrix.h  -  description
                             -------------------
    begin                : mer fév 11 2004
    copyright            : (C) 2004 by Olivier Roussel and 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 Matrix
{

/*
______________________________________________________________________________________________

    Constructor and distructor
______________________________________________________________________________________________

*/

public:
  Matrix();

Matrix(const int i, const int j);

Matrix(const int i);

Matrix(const Matrix& M);

Matrix(const Vector& V);

~Matrix();
/*
______________________________________________________________________________________________

    Set and get
______________________________________________________________________________________________

*/

inline void setValue(const int i, const int j, const real a);

void setZero();

inline real value(const int i, const int j) const;

inline int lines() const;

inline int columns() const;
/*
______________________________________________________________________________________________

    Operators
______________________________________________________________________________________________

*/

bool operator== (const Matrix& M ) const;

void operator= (const Matrix& M );

void operator+=(const Matrix& M );

Matrix operator+ (const Matrix& M ) const;

void operator-=(const Matrix& M );

Matrix operator- (const Matrix& M ) const;

Matrix operator- () const;

void operator*=(const real a );

Matrix operator* (const real a ) const;

void operator/=(const real a );

Matrix operator/ (const real a ) const;

Matrix operator* (const Matrix& M ) const;

Vector operator* (const Vector& V) const;

void setEigenMatrix(const bool isLeft, const int AxisNo, const Vector V, const real c, const real h=0.);
/*
______________________________________________________________________________________________

PRIVATE VARIABLES
______________________________________________________________________________________________

*/
//parallel modification
//private:
public:
  int  Lines, Columns;  
    real *U;    
};

/*
______________________________________________________________________________________________

EXTERNAL FUNCTIONS
______________________________________________________________________________________________

*/

Matrix operator* (const real a, const Matrix& M);

ostream& operator<<(ostream& out, const  Matrix& M);

/*
______________________________________________________________________________________________

INLINE FUNCTIONS
______________________________________________________________________________________________
*/

inline int Matrix::lines() const
{
    return Lines;
}
/*
______________________________________________________________________________________________

*/

inline int Matrix::columns() const
{
    return Columns;
}

/*
______________________________________________________________________________________________

*/

inline void Matrix::setValue(const int i, const int j, const real a)
{
        *( U + (i-1)*Columns + (j-1) ) = a;
}

/*
______________________________________________________________________________________________

*/

inline real Matrix::value(const int i, const int j) const
{
        return *(U+(i-1)*Columns+(j-1));
}