Manybody_Eigenvector.h

#ifndef QNANO_NEW_MANYBODY_EIGENVECTOR_DEFINED_H
#define QNANO_NEW_MANYBODY_EIGENVECTOR_DEFINED_H

#include "manybody/Fermion_Double_Index.h"
#include "tools/Reader.h"
#include "tools/List_Class.h"

typedef std::pair<Fermion_Double_Index, std::complex<double> > Manybody_Eigenvector_Element;

class Manybody_Eigenvector: public List_Class<Manybody_Eigenvector_Element>{
public:

  int nr_e, nr_h, rank_e, rank_h;
  double eigenvalue;
  
  void sort();
  void read(const std::string &filename);
  void print(const std::string &filename, int figure=8)const;
  int find_sorted(const Fermion_Double_Index &fi)const;
  void add_value(const Fermion_Double_Index &fi, std::complex<double> val);

  //Add perturbative corrections to eigenstate
  void correctPT(const std::string &PTfile, int line, double tolerance=1e-5);

  std::complex<double> dot(const Manybody_Eigenvector &other)const;
  
  void set_defaults(){
    nr_e=0; nr_h=0; rank_e=0; rank_h=0; eigenvalue=0;
  }
  void set_vacuum(){
    list.resize(1);
    list[0]=std::make_pair(Fermion_Double_Index(), 1.);
  }
  Manybody_Eigenvector(){
    set_defaults();
  }
  Manybody_Eigenvector(int nr_e_, int nr_h_, int rank_e_, int rank_h_, double eval=0.){
    nr_e=nr_e_;
    nr_h=nr_h_;
    rank_e=rank_e_;
    rank_h=rank_h_;
    eigenvalue=eval;
  }
  Manybody_Eigenvector(const std::string &filename){
    set_defaults();
    read(filename);
  }


};

#endif