potential_hamiltonian.hpp

#ifndef POT_HAMILT_DEC_HPP
#define POT_HAMILT_DEC_HPP

#include "rectangle.hpp"
#include "hamiltonian_decorator.hpp"
#include "funcs_base.hpp"

class Potential;
class Rectangle;


class ExtPot_Hamiltonian: public Hamiltonian_Decorator{
public:

    ExtPot_Hamiltonian(){};
    ~ExtPot_Hamiltonian() throw() {}

private:

    const char* name() const {return "ExtPot_Decorator";}; 

    typedef std::list<STD_TR1::shared_ptr< Potential > > T_Potentiale;
    STD_TR1::shared_ptr<T_Potentiale> Potentiale;

    void init_Decorator() ;

    void prae_create_H0() const;
    void post_create_H0() const;
    void prae_create_HI() const;
    void post_create_HI() const;
    void prae_create_HIt() const;
    void post_create_HIt() const;

};


class ChangeHop_Hamiltonian: public Hamiltonian_Decorator{
public:

    ChangeHop_Hamiltonian(){};
    ~ChangeHop_Hamiltonian() throw() {}

private:
    vector_Object::ref hopvec_v;
    vector_Object::ref hopvec_h;
    bool add_not_mult;
    mutable bool nowh0;

    const char* name() const {return "ChangeHop_Decorator";}; 

    void init_Decorator() ;

    void Inserting(row r, column c, complex v) const;
    
    void prae_create_H0() const;
    void post_create_H0() const;
    void prae_create_HI() const;
    void post_create_HI() const;
    void prae_create_HIt() const;
    void post_create_HIt() const;

};


class ChangeHop2D_Hamiltonian: public Hamiltonian_Decorator{
public:

    ChangeHop2D_Hamiltonian(){};
    ~ChangeHop2D_Hamiltonian() throw() {}

private:
    function2d_Object::ref hopfunc;
    mutable STD_TR1::shared_ptr<std::vector<double> > values;

    bool always_mult;
    bool add_not_mult;
    mutable bool nowh0;

    const char* name() const {return "ChangeHop2D_Decorator";}; 

    void init_Decorator() ;

    void update_values() const;

    void Inserting(row r, column c, complex v) const;
    
    void prae_create_H0() const;
    void post_create_H0() const;
    void prae_create_HI() const;
    void post_create_HI() const;
    void prae_create_HIt() const;
    void post_create_HIt() const;

};


class periodic_Hamiltonian: public simple_Hamiltonian{
public:
    
    periodic_Hamiltonian(){};
    ~periodic_Hamiltonian() throw() {};

private:
    int sub_ny;
    int nx;
    double big_dx;
    Rectangle::ref my_R;
    Hamiltonian::ref my_H;
    std::unique_ptr<pencil> get_pencil_() const;
    const std::string& get_pencil_name_() const;

    double dx_sh_() const;
    double dy_sh_() const;
    double width_sh_() const;
    point_Positions::ref rxy_sh_() const;
    size_t layers_() const {return my_H->layers();}
    int ny_layer_sh_(size_t) const;
    void createperiodicH0(Matrix&) const;
    void createperiodicHI(Matrix&) const;


    const char* name() const { return "periodic_Hamiltonian";}

    void init_H_params();
    void init_ny_dependent_H_params(){};
};

#endif