Program Listing for File lattice_matrix_element.hpp

↰ Return to documentation for file (src_qcm/lattice_matrix_element.hpp)

#ifndef lattice_matrix_element_h
#define lattice_matrix_element_h

#include "types.hpp"

using namespace std;

struct lattice_matrix_element
{
  size_t site1;
  size_t spin1;
  size_t site2;
  size_t spin2;
  size_t neighbor;
  complex<double> v;

  lattice_matrix_element() : site1(0), site2(0), spin1(0), spin2(0), neighbor(0), v(0.0) {}

  lattice_matrix_element(const int &_site1, const int &_spin1, const int &_site2, const int &_spin2, const int _neighbor, complex<double> _v)
  : site1(_site1), spin1(_spin1), site2(_site2), spin2(_spin2), neighbor(_neighbor), v(_v){}

  friend std::ostream & operator<<(std::ostream &flux, const lattice_matrix_element &s){
    if(s.v.imag() == 0) flux << '(' << s.site1+1 << (s.spin1? '-' : '+') <<  ',' << s.site2+1 << (s.spin2? '-' : '+') << ';' << s.neighbor << ") : " << s.v.real();
    else flux << '(' << s.site1+1 << (s.spin1? '-' : '+') <<  ',' << s.site2+1 << (s.spin2? '-' : '+') << ';' << s.neighbor << ") : " << s.v;
    return flux;
  }

  string str() const{
    ostringstream sout;
    sout << *this;
    return sout.str();
  }
};


struct lattice_index_pair
{
  static size_t Nc;
  size_t site1;
  size_t spin1;
  size_t site2;
  size_t spin2;
  size_t neighbor;
  size_t rank;

  lattice_index_pair() : site1(0), site2(0), spin1(0), spin2(0), neighbor(0), rank(0) {}

  lattice_index_pair(const int &_site1, const int &_spin1, const int &_site2, const int &_spin2, const int _neighbor)
  : site1(_site1), spin1(_spin1), site2(_site2), spin2(_spin2), neighbor(_neighbor){
    rank = site1+Nc*(spin1 + 2*(site2 + Nc*(spin2 + 2*neighbor)));
  }

  friend std::ostream & operator<<(std::ostream &flux, const lattice_index_pair &s){
    flux << '(' << s.site1+1 << (s.spin1? '-' : '+') <<  ',' << s.site2+1 << (s.spin2? '-' : '+') << ';' << s.neighbor << ")";
    return flux;
  }

  string str() const{
    ostringstream sout;
    sout << *this;
    return sout.str();
  }

  friend bool operator<(const lattice_index_pair &x, const lattice_index_pair &y){
    return x.rank < y.rank;
  }

};




template<typename HilbertField>
struct GF_matrix_element
{

  size_t r;
  size_t c;
  size_t n;
  HilbertField v;
};

#endif