Selected_Manybody_States_with_Hash_Table.h

#ifndef QNANO_NEW_SELECTED_MANYBODY_STATES_WITH_HASH_TABLE_DEFINED_H
#define QNANO_NEW_SELECTED_MANYBODY_STATES_WITH_HASH_TABLE_DEFINED_H

#include "manybody/Selected_Manybody_States.h"
#include "tools/Index_Hash_Table.h"


class Selected_Manybody_States_with_Hash_Table: public Selected_Manybody_States{
public:
  Index_Hash_Table<Fermion_Double_Index> HT;


  virtual size_t size()const{return HT.size();}
  void clear(){HT.clear();}
  bool is_contained(const Fermion_Double_Index &fi)const{
    return HT.get_index(fi)>=0;
  }

  int find(const Fermion_Double_Index &fi)const{
    return HT.get_index(fi);
  }
  int find_check(const Fermion_Double_Index &fi)const;
  Fermion_Double_Index find_reverse(int index)const{
    return HT.get_key(index);
  }
  
  void check_consistency()const{}

  void add_if_not_contained(const Fermion_Double_Index &fi){
    HT.add(fi);
  }
  void add_from_file(const std::string &fname, int rank_e, int rank_h, int max_nr_e, int max_nr_h);

  void print(std::ostream &os=std::cout)const;
  void print(const std::string &fname)const{
    std::ofstream ofs(fname.c_str());
    print(ofs);
  }

  void add_all(int rank_e, int rank_h, int nr_e, int nr_h, int start_e=0, int start_h=0);
  void add_singles(int nr_e, int nr_h);


  Selected_Manybody_States_with_Hash_Table(const Selected_Manybody_States_with_Hash_Table &other) 
   : HT(other.HT){
  }
  Selected_Manybody_States_with_Hash_Table(size_t start_size=1024) : HT(Hash_Fermion_Double_Index,start_size){
  }

};
#endif