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