Selected_Manybody_Hilbert_Space.h

#ifndef QNANO_NEW_SELECTED_MANYBODY_HILBERT_SPACE_DEFINED_H
#define QNANO_NEW_SELECTED_MANYBODY_HILBERT_SPACE_DEFINED_H

#include "manybody/Manybody_Hilbert_Space.h"
#include "manybody/Selected_Manybody_States_with_Hash_Table.h"

class Selected_Manybody_Hilbert_Space: public Manybody_Hilbert_Space{
public:
  Selected_Manybody_States_with_Hash_Table states;      

  virtual void setup(Parameter_Map &param_map, bool do_print=false);
  virtual size_t size()const{return states.size();}
  virtual int get_index(const Fermion_Double_Index & fi)const;
  virtual Fermion_Double_Index get_Fermion_Double_Index(int index)const;

  void add_singles(){
    states.add_singles(nr_e,nr_h);
  }

  virtual Selected_Manybody_Hilbert_Space* clone() const;
  //Copy from  Selected_Manybody_Hilbert_Space:
  virtual void copy(const Selected_Manybody_Hilbert_Space &other){
    rank_e=other.rank_e;
    rank_h=other.rank_h;
    nr_e=other.nr_e;
    nr_h=other.nr_h;
    states=other.states;
  }
  Selected_Manybody_Hilbert_Space &operator=(const Selected_Manybody_Hilbert_Space &other){
    copy(other);
    return *this;
  }
  Selected_Manybody_Hilbert_Space(const Selected_Manybody_Hilbert_Space &other){
    copy(other);
  }

  virtual void copy(const Manybody_Hilbert_Space &other){
    rank_e=other.rank_e;
    rank_h=other.rank_h;
    nr_e=other.nr_e;
    nr_h=other.nr_h;
    states.clear();
    for(size_t i=0; i<other.size(); i++){
      Fermion_Double_Index fi=other.get_Fermion_Double_Index(i);
      states.add_if_not_contained(fi);
    }
  }
  Selected_Manybody_Hilbert_Space &operator=(const Manybody_Hilbert_Space &other){
    copy(other);
    return *this;
  }
  Selected_Manybody_Hilbert_Space(const Manybody_Hilbert_Space &other){
    copy(other);
  }
 

  //Constructor from parameters
  Selected_Manybody_Hilbert_Space(Parameter_Map &param_map, bool do_print=false){
    setup(param_map, do_print);
  }
};

#endif