LUDecomposition.h

#ifndef MATH_LU_DECOMPOSITION_H
#define MATH_LU_DECOMPOSITION_H

#include "MatrixTemplate.h"
#include <vector>

namespace Math {

template <class T>
class LUDecomposition
{
public:
  typedef MatrixTemplate<T> MatrixT;
  typedef VectorTemplate<T> VectorT;

  LUDecomposition();
  LUDecomposition(const MatrixT& A);
  
  bool set(const MatrixT& A);
  void backSub(const VectorT& b, VectorT& x) const;
  void PBackSub(const VectorT& b, VectorT& x) const;
  void LBackSub(const VectorT& b, VectorT& x) const;
  void UBackSub(const VectorT& b, VectorT& x) const;
  void getInverse(MatrixT& Ainv) const;
  
  void getL(MatrixT& L) const;
  void getU(MatrixT& U) const;

  MatrixT LU;
  std::vector<int> P;
  T zeroTolerance;
};

}
#endif