hamiltonian_decorator.hpp

Go to the documentation of this file.

#ifndef H_DECORATOR_HEADER
#define H_DECORATOR_HEADER

#include "hamiltonian.hpp"
#include "matrix.hpp"

class Hamiltonian_Decorator: public Hamiltonian{

public:

    Hamiltonian_Decorator() : my_m(NULL), my_pos(NULL) {};
    ~Hamiltonian_Decorator() throw() {};

    
protected:
    
    Hamiltonian::const_ref get_H() const { return my_H;};
    Matrix* get_m() const { return my_m;};
    const geometry_const_Iterator* get_pos() const {return my_pos;};

private:

    std::unique_ptr<pencil> get_pencil_() const;
    const std::string& get_pencil_name_() const;

    void createH0(Matrix&,const geometry_const_Iterator& pos) const;
    void createHI(Matrix&,const geometry_const_Iterator& pos) const;
    void createHIt(Matrix&,const geometry_const_Iterator& pos) const;
    size_t layers_() const {return my_H->layers();};

    int ny_(const geometry_const_Iterator& pos)        const {return my_H->ny(pos);};
    int ny_layer_(const geometry_const_Iterator& pos)  const {return my_H->ny_layer(pos);};
    double dx_(const geometry_const_Iterator& pos) const {return my_H->dx(pos);};
    double dy_(const geometry_const_Iterator& pos) const {return my_H->dy(pos);};
    double width_(const geometry_const_Iterator& pos)  const {return my_H->width(pos);};

    virtual point_Positions::ref rxy_(const geometry_const_Iterator& pos) const;

    template <class T> 
    friend class custom_HD_Matrix;

    virtual void Inserting(row r, column c, complex v) const;
    virtual void Adding(row r, column c, complex v) const;
    
    virtual complex Get(row r, column c) const;

    virtual void prae_create_H0() const= 0;
    virtual void post_create_H0() const= 0;
    virtual void prae_create_HI() const= 0;
    virtual void post_create_HI() const= 0;
    virtual void prae_create_HIt() const = 0;
    virtual void post_create_HIt() const = 0;

    void init();
    virtual void init_Decorator() = 0;

    Hamiltonian::ref my_H;
    mutable Matrix* my_m;
    mutable const geometry_const_Iterator* my_pos;
   
};

class simple_Hamiltonian_Decorator: public simple_Hamiltonian{

public:

    simple_Hamiltonian_Decorator() : my_m(NULL) {};
    ~simple_Hamiltonian_Decorator() throw() {};
    
    int lead_system_handler() const {return get_H()->lead_system_handler();};
    
protected:
    
    simple_Hamiltonian::const_ref get_H() const { return my_H;};
    Matrix* get_m() const { return my_m;};

private:

    bool do_generalized_() const { return my_H->do_generalized();}

    std::unique_ptr<pencil> get_pencil_() const;
    const std::string& get_pencil_name_() const;

    void createperiodicH0(Matrix&) const;
    void createperiodicHI(Matrix&) const;

    size_t layers_() const {return my_H->layers();};

    int ny_layer_sh_(size_t i) const { return my_H->ny_layer(i); }
    double dx_sh_() const {return my_H->dx();};
    double dy_sh_() const {return my_H->dy();};
    double width_sh_()  const {return my_H->width();};

    virtual point_Positions::ref rxy_sh_() const;

    template <class T> 
    friend class custom_HD_Matrix;

    virtual void Inserting(row r, column c, complex v) const;
    virtual void Adding(row r, column c, complex v) const;

    virtual complex Get(row r, column c) const;

    virtual void prae_create_H0() const= 0;
    virtual void post_create_H0() const= 0;
    virtual void prae_create_HI() const= 0;
    virtual void post_create_HI() const= 0;

    void init_H_params();
    void init_ny_dependent_H_params();

    virtual void init_Decorator() = 0;
    virtual void init_ny_dependent_Decorator() = 0;

    simple_Hamiltonian::ref my_H;
    mutable Matrix* my_m;
   
};

#endif