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