QRDecomposition.h

#ifndef MATH_QR_DECOMPOSITION_H
#define MATH_QR_DECOMPOSITION_H

#include "MatrixTemplate.h"

namespace Math {

template <class T>
class QRDecomposition
{
public:
  typedef MatrixTemplate<T> MatrixT;
  typedef VectorTemplate<T> VectorT;
  bool set(const MatrixT& A);
  void backSub(const VectorT& x,VectorT& b) const;
  void leastSquares(const VectorT& x,VectorT& b,VectorT& residual) const;
  void QMul(const VectorT& b,VectorT& x) const;
  void QtMul(const VectorT& b,VectorT& x) const;
  void RBackSub(const VectorT& b,VectorT& x) const;
  void getQ(MatrixT& Q) const;
  void getR(MatrixT& R) const;
  
  MatrixT QR;
  VectorT tau;
};

template <class T>
class NRQRDecomposition
{
public:
  typedef MatrixTemplate<T> MatrixT;
  typedef VectorTemplate<T> VectorT;
  bool set(const MatrixT& A);
  void backSub(const VectorT& x,VectorT& b) const;
  void QBackSub(const VectorT& x,VectorT& b) const;
  void getQ(MatrixT& Q) const;
  
  MatrixT QR;
  VectorT c;
  bool singular;
};

} //namespace Math

#endif