GenericHouseholder.h

Go to the documentation of this file.

#ifndef GenericHouseholder_h
#define GenericHouseholder_h
 
#include <vector>
#include <iostream>
 
class GenericHouseholder {
public:
  GenericHouseholder(bool normalise = false);
 
  ~GenericHouseholder();
 
  std::vector<float> iterate(const std::vector<std::vector<float> >& eventMatrix,
                             const std::vector<float>& energyVector,
                             const int nIter);
 
  std::vector<float> iterate(const std::vector<std::vector<float> >& eventMatrix,
                             const std::vector<float>& energyVector);
 
private:
  bool decompose(const int m,
                 const int n,
                 std::vector<std::vector<float> >& qr,
                 std::vector<float>& alpha,
                 std::vector<int>& pivot);
 
  void solve(int m,
             int n,
             const std::vector<std::vector<float> >& qr,
             const std::vector<float>& alpha,
             const std::vector<int>& pivot,
             std::vector<float>& r,
             std::vector<float>& y);
 
  bool normaliseFlag;
};
 
#endif