core/Matrix.hpp

Go to the documentation of this file.

/*
  This file is part of OpenFLUID software
  Copyright (c) 2007-2010 INRA-Montpellier SupAgro


 == GNU General Public License Usage ==

  OpenFLUID 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 3 of the License, or
  (at your option) any later version.

  OpenFLUID 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 OpenFLUID.  If not, see <http://www.gnu.org/licenses/>.

  In addition, as a special exception, INRA gives You the additional right
  to dynamically link the code of OpenFLUID with code not covered
  under the GNU General Public License ("Non-GPL Code") and to distribute
  linked combinations including the two, subject to the limitations in this
  paragraph. Non-GPL Code permitted under this exception must only link to
  the code of OpenFLUID dynamically through the OpenFLUID libraries
  interfaces, and only for building OpenFLUID plugins. The files of
  Non-GPL Code may be link to the OpenFLUID libraries without causing the
  resulting work to be covered by the GNU General Public License. You must
  obey the GNU General Public License in all respects for all of the
  OpenFLUID code and other code used in conjunction with OpenFLUID
  except the Non-GPL Code covered by this exception. If you modify
  this OpenFLUID, you may extend this exception to your version of the file,
  but you are not obligated to do so. If you do not wish to provide this
  exception without modification, you must delete this exception statement
  from your version and license this OpenFLUID solely under the GPL without
  exception.


 == Other Usage ==

  Other Usage means a use of OpenFLUID that is inconsistent with the GPL
  license, and requires a written agreement between You and INRA.
  Licensees for Other Usage of OpenFLUID may use this file in accordance
  with the terms contained in the written agreement between You and INRA.
*/


#ifndef __MATRIX_HPP___
#define __MATRIX_HPP___



#include <iostream>
#include <openfluid/dllexport.hpp>
#include <openfluid/base/OFException.hpp>
#include <boost/multi_array.hpp>

namespace openfluid { namespace core {


template <class T>
class DLLEXPORT Matrix
{
  protected :

    boost::multi_array<T, 2> m_Data;

    unsigned long m_ColsNbr;

    unsigned long m_RowsNbr;


  public :

    Matrix();

    Matrix(const Matrix &Matrix);

    Matrix(unsigned long ColsNbr,unsigned long RowsNbr);

    Matrix(unsigned long ColsNbr,unsigned long RowsNbr, T InitValue);

    virtual ~Matrix() {};

    inline unsigned long getColsNbr() const { return m_ColsNbr; };

    inline unsigned long getRowsNbr() const { return m_RowsNbr; };

    inline unsigned long getSize() const { return (m_ColsNbr * m_RowsNbr); };

    inline unsigned long size() const { return (m_ColsNbr * m_RowsNbr); };


    T* getData() const { return (T*)(m_Data.data()); };

    void setData(T* Data);

    T getElement(unsigned long ColIndex, unsigned long RowIndex) const;

    inline T at(unsigned long ColIndex, unsigned long RowIndex) const { return getElement(ColIndex,RowIndex); };

    inline T get(unsigned long ColIndex, unsigned long RowIndex) const { return getElement(ColIndex,RowIndex); };


    void setElement(unsigned long ColIndex, unsigned long RowIndex, T Element);

    inline void set(unsigned long ColIndex, unsigned long RowIndex, T Element) { setElement(ColIndex,RowIndex,Element); };

    Matrix<T>& operator = (const Matrix &A);

    void fill(const T& Val);

    void clear();

};

// =====================================================================
// =====================================================================

template <class T>
Matrix<T>::Matrix():
  m_Data(boost::extents[0][0],boost::fortran_storage_order()),
  m_ColsNbr(0),m_RowsNbr(0)
{
}


// =====================================================================
// =====================================================================

template <class T>
Matrix<T>::Matrix(const Matrix &A):
  m_Data(boost::extents[A.m_ColsNbr][A.m_RowsNbr],boost::fortran_storage_order()),
  m_ColsNbr(A.m_ColsNbr),m_RowsNbr(A.m_RowsNbr)
{
  m_Data = A.m_Data;
}


// =====================================================================
// =====================================================================

template <class T>
Matrix<T>::Matrix(unsigned long ColsNbr, unsigned long RowsNbr) :
  m_Data(boost::extents[ColsNbr][RowsNbr],boost::fortran_storage_order()),
  m_ColsNbr(ColsNbr),m_RowsNbr(RowsNbr)
{
  fill(0);
}


// =====================================================================
// =====================================================================


template <class T>
Matrix<T>::Matrix(unsigned long ColsNbr, unsigned long RowsNbr, T InitValue) :
  m_Data(boost::extents[ColsNbr][RowsNbr],boost::fortran_storage_order()),
  m_ColsNbr(ColsNbr),m_RowsNbr(RowsNbr)
{
  fill(InitValue);
}


// =====================================================================
// =====================================================================


template <class T>
void Matrix<T>::setData(T* Data)
{
  for (unsigned long j=0; j < m_RowsNbr;j++)
  {
    for (unsigned long i=0; i < m_ColsNbr;i++)
    {
      m_Data[i][j] = Data[i+(j*m_ColsNbr)];
    }
  }
}


// =====================================================================
// =====================================================================


template <class T>
T Matrix<T>::getElement(unsigned long ColIndex, unsigned long RowIndex) const
{
  if (ColIndex >= m_ColsNbr || RowIndex >= m_RowsNbr) throw openfluid::base::OFException("OpenFLUID framework","Matrix::getElement","element access range error");
  return m_Data[ColIndex][RowIndex];
}


// =====================================================================
// =====================================================================


template <class T>
void Matrix<T>::setElement(unsigned long ColIndex, unsigned long RowIndex, T Element)
{
  if (ColIndex >= m_ColsNbr || RowIndex >= m_RowsNbr) throw openfluid::base::OFException("OpenFLUID framework","Matrix::setElement","element access range error");
  m_Data[ColIndex][RowIndex] = Element;
}



// =====================================================================
// =====================================================================

template <class T>
Matrix<T>& Matrix<T>::operator=(const Matrix &A)
{

  if (this == &A) return *this; // in case somebody tries assign array to itself

  m_Data.resize(boost::extents[A.m_ColsNbr][A.m_RowsNbr]);
  m_Data = A.m_Data;
  m_ColsNbr = A.m_ColsNbr;
  m_RowsNbr = A.m_RowsNbr;

  return *this;
}


// =====================================================================
// =====================================================================


template <class T>
void Matrix<T>::fill(const T& Val)
{
  for (unsigned long i=0;i<m_ColsNbr;i++)
    for (unsigned long j=0;j<m_RowsNbr;j++)
      m_Data[i][j] = Val;
}


// =====================================================================
// =====================================================================

template <class T>
void Matrix<T>::clear()
{
  m_Data.resize(boost::extents[0][0]);
  m_ColsNbr=0;
  m_RowsNbr=0;
}



} }

#endif /* __MATRIX_HPP___ */