CdaggerC_Iterator.h

#ifndef QNANO_NEW_CDAGGERC_ITERATOR_DEFINED_H
#define QNANO_NEW_CDAGGERC_ITERATOR_DEFINED_H

#include "manybody/CdaggerC_Iterator_SingleI.h"

class CdaggerC_Iterator{
protected:
  Fermion_Index fi_ref;
  int n_max;
  int cur_col;
 
  CdaggerC_Iterator_SingleI it;

public: 

  bool done()const{
    return cur_col>=(int)fi_ref.size();
  }
  inline int get_n()const{return it.get_n();}
  inline int get_m()const{return it.get_m();}
  inline double get_sign()const{return it.get_sign();}
  inline const Fermion_Index& get_fi()const{return it.get_fi();}

  CdaggerC_Iterator& operator++(){
    ++it; 
    if(it.done()){
      cur_col++;
      if(cur_col<(int)fi_ref.size()){
        it.setup(fi_ref, cur_col, n_max);
      }
    }
    return *this;
  }
  void setup(const Fermion_Index &fi, int n_max_){
    fi_ref=fi;
    n_max=n_max_;
    cur_col=0;
    it.setup(fi, cur_col, n_max);
  }

  CdaggerC_Iterator(){
    setup(Fermion_Index(), 0);
  }
  CdaggerC_Iterator(const Fermion_Index &fi, int n_max_){
    setup(fi, n_max_);
  }
};

#endif