lowerTriangularMatrix.h

Go to the documentation of this file.

#ifndef _oblas_lowerTriangularMatrix_h_
#define _oblas_lowerTriangularMatrix_h_

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


namespace oBLAS {

  namespace ublas = boost::numeric::ublas;


  class lowerTriangularMatrix : public ublas::triangular_matrix<double, ublas::lower> {

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

  private:
    mutable lowerTriangularMatrix* fInverse = nullptr;

  };
}


#endif