upperTriangularMatrix.h

Go to the documentation of this file.

#ifndef _oblas_upperTriangularMatrix_h_
#define _oblas_upperTriangularMatrix_h_

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


namespace oBLAS {

  namespace ublas = boost::numeric::ublas;


  class upperTriangularMatrix : public ublas::triangular_matrix<double, ublas::upper> {
    
  public:
    upperTriangularMatrix(int i);
    upperTriangularMatrix(const ublas::matrix<double>& r);
    ~upperTriangularMatrix();
    void operator=(const ublas::matrix<double>& r);
    const upperTriangularMatrix* GetInverse();

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

  };

}


#endif