Finite_Manybody_Hilbert_Space.h

#ifndef QNANO_NEW_FINITE_MANYBODY_HILBERT_SPACE_DEFINED_H
#define QNANO_NEW_FINITE_MANYBODY_HILBERT_SPACE_DEFINED_H

#include "tools/Parameter_Map.h"
#include "tools/Organize_Blocks.h"
#include "manybody/Fermion_Double_Index.h"
#include "manybody/Manybody_Hilbert_Space.h"
#include "tools/Binomial.h"

class Finite_Manybody_Hilbert_Space: public Manybody_Hilbert_Space {
protected: 
//  int rank_e, rank_h;  ///nr of electrons/holes
//->Manybody_Hilbert_Space  int nr_e, nr_h;      ///nr of electron/hole _states_

  int dim_e, dim_h;    

public:
//  inline int get_rank_e()const{return rank_e;}
//  inline int get_rank_h()const{return rank_h;}
//->Manybody_Hilbert_Space  inline int get_nr_e()const{return nr_e;}
//->Manybody_Hilbert_Space  inline int get_nr_h()const{return nr_h;}
  inline int get_dim_e()const{return dim_e;}
  inline int get_dim_h()const{return dim_h;}
  virtual size_t size()const{return get_dim_e()*get_dim_h();}

  virtual void resize(int rk_e, int rk_h, int n_e, int n_h){
    Manybody_Hilbert_Space::resize(rk_e,rk_h,n_e,n_h);
//    rank_e=rk_e; rank_h=rk_h; nr_e=n_e; nr_h=n_h;
    recalc_dims();
  }
  virtual void setup(Parameter_Map &param_map){
    Manybody_Hilbert_Space::setup(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);
*/
    recalc_dims();
  }
   

  // recall: get_DIM inherited from Organiz_Blocks

  
  virtual void recalc_dims();

  virtual int get_Index_from_Fermion_Index(const Fermion_Index & fi, bool hole)const;
  virtual int get_index(const Fermion_Double_Index & fi)const;

  virtual Fermion_Index get_Fermion_Index_from_Index(int index, bool hole)const;
  virtual Fermion_Double_Index get_Fermion_Double_Index(int index)const;

  //test the Fermion_Double_Index <-> index relation
  virtual void dump_indices(const std::string &filename)const;

  virtual Finite_Manybody_Hilbert_Space* clone() const;
  virtual void copy(const Finite_Manybody_Hilbert_Space &other){
    rank_e=other.rank_e;
    rank_h=other.rank_h;
    nr_e=other.nr_e;
    nr_h=other.nr_h;
    dim_e=other.dim_e;
    dim_h=other.dim_h;
  }
  Finite_Manybody_Hilbert_Space &operator=(const Finite_Manybody_Hilbert_Space &other){
    copy(other);
    return *this;
  }
  Finite_Manybody_Hilbert_Space(const Finite_Manybody_Hilbert_Space &other){
    copy(other);
  }
  Finite_Manybody_Hilbert_Space(){
    rank_e=rank_h=nr_e=nr_h=dim_e=dim_h=0;
  }
  Finite_Manybody_Hilbert_Space(Parameter_Map &param_map){
    setup(param_map);
  }
  Finite_Manybody_Hilbert_Space(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;
    recalc_dims();
  }

  virtual bool is_Full_Hilbert_Space() const;
};

#endif