Manybody_Hilbert_Space.h

#ifndef QNANO_NEW_MANYBODY_HILBERT_SPACE_DEFINED_H
#define QNANO_NEW_MANYBODY_HILBERT_SPACE_DEFINED_H

#include "tools/Parameter_Map.h"
#include "manybody/Fermion_Double_Index.h"
#include "tools/Size_Provider.h"

class Manybody_Hilbert_Space: public Size_Provider{
public:
  int rank_e, rank_h;  
  int nr_e, nr_h; 

public:
  inline int get_rank_e()const{return rank_e;}
  inline int get_rank_h()const{return rank_h;}
  inline int get_nr_e()const{return nr_e;}
  inline int get_nr_h()const{return nr_h;}
  //virtual size_t size()const=0;

  virtual void resize(int rk_e, int rk_h, int n_e, int n_h){
    rank_e=rk_e; rank_h=rk_h; nr_e=n_e; nr_h=n_h;
  }
  virtual void setup(Parameter_Map &param_map){
    rank_e=param_map.get_as_double("rank_e",0);
    rank_h=param_map.get_as_double("rank_h",0);
    nr_e=param_map.get_as_double("nr_e",0);
    nr_h=param_map.get_as_double("nr_h",0);
  }
 
  virtual int get_index(const Fermion_Double_Index & fi)const=0;
  virtual Fermion_Double_Index get_Fermion_Double_Index(int index)const=0;
  static std::string print_index(size_t line_counter, void* spi_ctx);


  virtual Manybody_Hilbert_Space* clone()const =0;
  Manybody_Hilbert_Space(){}
  Manybody_Hilbert_Space(Parameter_Map &param_map){
    setup(param_map);
  }
  virtual ~Manybody_Hilbert_Space(){}

  virtual bool is_Full_Hilbert_Space() const;
};

#endif