RowEchelon.h

#ifndef MATH_ROW_ECHELON_H
#define MATH_ROW_ECHELON_H

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

namespace Math {

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

  RowEchelon();
  RowEchelon(const MatrixT& A);
  RowEchelon(const MatrixT& A,const VectorT& b);
  RowEchelon(const MatrixT& A,const MatrixT& B);
  
  void set(const MatrixT& A);
  void set(const MatrixT& A,const VectorT& b);
  void set(const MatrixT& A,const MatrixT& B);
  void backSub(VectorT& x) const;
  int getRank() const;
  int getNull() const;
  void getNullspace(MatrixT& N) const;
  void getAllSolutions(VectorT& x0,MatrixT& N) const;

  //helpers
  void calcFirstEntries();
  
  MatrixT R;
  MatrixT EB;
  std::vector<int> firstEntry;  
};

template <class T>
bool IsRowEchelon(const MatrixTemplate<T>& A);
template <class T>
bool IsReducedRowEchelon(const MatrixTemplate<T>& A);

} //namespace Math


#endif