TB_Hopping_Parameter.h

#ifndef QNANO_NEW_TB_HOPPING_PARAMETER_DEFINED_H
#define QNANO_NEW_TB_HOPPING_PARAMETER_DEFINED_H

#include "tools/Reader.h"


class TB_Hopping_Parameter{
public:
  static const size_t PARVEC_SIZE=(21+21);
  union {
    double parvec[PARVEC_SIZE];
    struct{
    double T_ss,T_sp,T_ps;
    double T_pp_sigma,T_pp_pi;
    double T_ssstar, T_sstars,T_psstar,T_sstarp,T_sstarsstar;
    double T_sd,T_ds, T_pd_sigma,T_dp_sigma,T_pd_pi,T_dp_pi;
    double T_sstard,T_dsstar;
    double T_dd_sigma,T_dd_pi,T_dd_delta;

    double N_ss,N_sp,N_ps;
    double N_pp_sigma,N_pp_pi;
    double N_ssstar, N_sstars,N_psstar,N_sstarp,N_sstarsstar;
    double N_sd,N_ds, N_pd_sigma,N_dp_sigma,N_pd_pi,N_dp_pi;
    double N_sstard,N_dsstar;
    double N_dd_sigma,N_dd_pi,N_dd_delta;
    };
  };

  //Get pointer to variable by name:
  double * name_to_ptr(const std::string &str);
  double * name_to_ptr_check(const std::string &str, const std::string &complain="");


  void print(int e1, int e2, std::ostream &os=std::cout)const;
  void print(int e1, int e2, const std::string &filename)const;

  void clear(){
    for(size_t i=0; i<PARVEC_SIZE; i++)parvec[i]=0;
  }
  TB_Hopping_Parameter(){
    clear();
  }

};

#endif