eigensolver_IRLanczos.h

Go to the documentation of this file.

#ifndef EIGENSOLVER_IRLANCZOS_INCLUDED
#define EIGENSOLVER_IRLANCZOS_INCLUDED

#include "eigensolver.h"
#include "Fopr/fopr.h"
#include "Tools/sorter.h"
#include "bridge_complex.h"

#include "IO/bridgeIO.h"
using Bridge::vout;



class Eigensolver_IRLanczos : public Eigensolver
{
 public:
  static const std::string class_name;

 private:

  int    m_Nk;
  int    m_Np;
  int    m_Niter_eigen;
  double m_Enorm_eigen;
  double m_Vthreshold;

  Fopr   *m_fopr;
  Sorter *m_sorter;

 public:
  Eigensolver_IRLanczos(Fopr *fopr)
    : Eigensolver(), m_fopr(fopr), m_sorter(0) {}

  Eigensolver_IRLanczos(unique_ptr<Fopr>& fopr)
    : Eigensolver(), m_fopr(fopr.get()), m_sorter(0) {}

  ~Eigensolver_IRLanczos();

  void set_parameters(const Parameters& params);
  void set_parameters(int Nk, int Np, int Niter_eigen, double Enorm_eigen,
                      double Vthreshold);
  void set_parameters(const std::string& sort_type,
                      int Nk, int Np, int Niter_eigen, double Enorm_eigen,
                      double Vthreshold);

  void solve(std::vector<double>& TDa, std::vector<Field>& vk,
             int& Nsbt, int& Nconv, const Field& b);

 private:

  void step(int Nm, int k,
            std::vector<double>& TDa,
            std::vector<double>& TDb,
            std::vector<Field>& vk, Field& f);

  void qrtrf(std::vector<double>& TDa, std::vector<double>& TDb,
             int Nk, int Nm, std::vector<double>& Qt,
             double Dsh, int kmin, int kmax);

  void tqri(std::vector<double>& TDa, std::vector<double>& TDb,
            int Nk, int Nm, std::vector<double>& Qt);

  void setUnit_Qt(int Nm, std::vector<double>& Qt);

  void schmidt_orthogonalization(Field& w, std::vector<Field>& vk, int k);
};
#endif