diagonalMatrix.h

Go to the documentation of this file.

#ifndef _oblas_diagonalMatrix_h_
#define _oblas_diagonalMatrix_h_

#ifndef NDEBUG
#  define NDEBUG
#endif
#include <boost/numeric/ublas/banded.hpp>


namespace oBLAS {

  namespace ublas = boost::numeric::ublas;


  class diagonalMatrix : public ublas::banded_matrix<double> {

  public:
    diagonalMatrix(int i);
    diagonalMatrix(const ublas::matrix<double>& r);
    ~diagonalMatrix();
    void operator=(const ublas::matrix<double>& r);
    const diagonalMatrix* GetInverse();

  private:
    void Init();
    diagonalMatrix();             // no default constuctor ...
    diagonalMatrix* fInverse;

  };

}


#endif