Manybody_Hamiltonian.h

#ifndef QNANO_NEW_MANYBODY_HAMILTONIAN_DENSE_DEFINED_H
#define QNANO_NEW_MANYBODY_HAMILTONIAN_DENSE_DEFINED_H

#include "tools/QNANO_Matrix.h"
#include "tools/Sparse_Matrix_Generator.h"
#include "manybody/Finite_Manybody_Hilbert_Space.h"
#include "manybody/Selected_Manybody_Hilbert_Space.h"
#include "manybody/Coulomb_Matrix_Elements.h"
#include "manybody/Manybody_Eigenvector.h"
#include "tools/Map.h"
#include "manybody/CdaggerCdaggerCC_Iterator.h"


class Manybody_Hamiltonian: public virtual Sparse_Matrix_Generator, public virtual Organize_Blocks_Single_Forward {
public:
  bool no_dir, no_exc, no_ee, no_hh;
  Coulomb_Matrix_Elements Vdir, Vexc, Vee, Vhh;

  bool use_offdiag_energies_e; 
  bool use_offdiag_energies_h; 
  std::vector<double> cb_energies;
  std::vector<double> vb_energies;
  QNANO_Matrix cb_energies_mat;
  QNANO_Matrix vb_energies_mat;


  Manybody_Hilbert_Space *hilbert;

  void setup(Parameter_Map &param_map, bool do_print=false);
  void setup_energies(Parameter_Map &param_map);
  void setup_Organize_Blocks();
  void replace_hilbert_space(const Manybody_Hilbert_Space &hs);

  virtual void get_Matrix_row_block(  std::vector< col_val > & col_val_vec, const Fermion_Double_Index &fi);
  virtual void get_Matrix_row_block(  std::vector< col_val > & col_val_vec, const int atom_index);

  //helper functions for get_Matrix_row_block
  void add_eh_terms(Map<std::complex<double> > &col_val_map, const Fermion_Double_Index &fi)const;
  void add_ee_terms(Map<std::complex<double> > &col_val_map, const Fermion_Double_Index &fi)const;
  void add_hh_terms(Map<std::complex<double> > &col_val_map, const Fermion_Double_Index &fi)const;

  double get_diagonal(const Fermion_Double_Index &fi)const;

  Manybody_Eigenvector apply(const Manybody_Eigenvector &EV);

  //Constructors
  void copy(const Manybody_Hamiltonian &other){
    //Sparse_Matrix_Generator
    Organize_Blocks_Single_Forward::copy(other);
    
    no_dir=other.no_dir;
    no_exc=other.no_exc;
    no_ee=other.no_ee;
    no_hh=other.no_hh;
    Vdir.copy(other.Vdir);
    Vexc.copy(other.Vexc);
    Vee.copy(other.Vee);
    Vhh.copy(other.Vhh);
    cb_energies=other.cb_energies;
    vb_energies=other.vb_energies;

    hilbert=other.hilbert->clone();
  }
  Manybody_Hamiltonian &operator=(const Manybody_Hamiltonian &other){
    copy(other);
    return *this;
  }
  Manybody_Hamiltonian(const Manybody_Hamiltonian &other){
    copy(other);
  }
  Manybody_Hamiltonian() : hilbert(NULL){
    Parameter_Map param_map;
    setup(param_map);
  }
  Manybody_Hamiltonian(Parameter_Map &param_map) : hilbert(NULL){
    setup(param_map);
  }
  
  virtual ~Manybody_Hamiltonian(){
    if(hilbert!=NULL)delete hilbert;
  }
};

#endif