/Users/barrand/private/dev/softinex/old/inexlib-1.2/inlib/inlib/a4z

Go to the documentation of this file.

// Copyright (C) 2010, Guy Barrand. All rights reserved.
// See the file inlib.license for terms.

#ifndef inlib_a4z
#define inlib_a4z

#include "a4"

#include <complex>

namespace inlib {
namespace a4 {

inline sqm< std::complex<double> > to_complex(const sqm<double>& a_from) {
  typedef std::complex<double> z;
  sqm< std::complex<double> > to;
  const std::vector<double>& fvec = a_from.vector();
  std::vector<z>& tvec = to.vector();
  typedef std::vector<double>::const_iterator const_vec_it_t;
  const_vec_it_t fit = fvec.begin();
  typedef std::vector<z>::iterator vec_it_t;
  vec_it_t tit = tvec.begin();
  for(;fit!=fvec.end();++fit,++tit) *tit = *fit;
  return to;
}

}}

namespace inlib {
namespace a4z {


class id : public a4::sqm< std::complex<double> > {
public:
  id(){
    typedef std::complex<double> z;
    v00(z(1,0));
    v11(z(1,0));
    v22(z(1,0));
    v33(z(1,0));
  }
};

//  Standard Dirac matrices :
//  D0               D1               D2              D3
//    1  0  0  0       0  0  0  1       0  0  0 -i      0  0  1  0
//    0  1  0  0       0  0  1  0       0  0  i  0      0  0  0 -1
//    0  0 -1  0       0 -1  0  0       0  i  0  0     -1  0  0  0
//    0  0  0 -1      -1  0  0  0      -i  0  0  0      0  1  0  0
class D0 : public a4::sqm< std::complex<double> > {
public:
  D0(){
    v00( 1);
    v11( 1);
    v22(-1);
    v33(-1);
  }
};
class D1 : public a4::sqm< std::complex<double> > {
public:
  D1(){
    v03( 1);
    v12( 1);
    v21(-1);
    v30(-1);
  }
};
class D2 : public a4::sqm< std::complex<double> > {
public:
  D2(){
    typedef std::complex<double> z;
    z i(0,1);
    v03(-i);
    v12( i);
    v21( i);
    v30(-i);
  }
};
class D3 : public a4::sqm< std::complex<double> > {
public:
  D3(){
    v02( 1);
    v13(-1);
    v20(-1);
    v31( 1);
  }
};

// Lorentz rep (1/2,1/2)
// HH1             HH2              HH3
//  0  1  1  0       0 -i -i  0       2  0  0  0
//  1  0  0  1       i  0  0 -i       0  0  0  0
//  1  0  0  1       i  0  0 -i       0  0  0  0
//  0  1  1  0       0  i  i  0       0  0  0 -2
// the uppers must be multiplied by : -i*0.5

// HH4             HH5              HH6
//  0  0  0  0       0  i -i  0       0 -1  1  0
//  0  2  0  0      -i  0  0 -i      -1  0  0  1
//  0  0 -2  0       i  0  0  i       1  0  0 -1
//  0  0  0  0       0  i -i  0       0  1 -1  0
// the uppers must be multiplied by : 0.5

class HH1 : public a4::sqm< std::complex<double> > {
public:
  HH1(){
    typedef std::complex<double> z;
    z i(0,1);
    v01( 1);
    v02( 1);
    v10( 1);
    v13( 1);
    v20( 1);
    v23( 1);
    v31( 1);
    v32( 1);
    multiply(-i*0.5);
  }
};
class HH2 : public a4::sqm< std::complex<double> > {
public:
  HH2(){
    typedef std::complex<double> z;
    z i(0,1);
    v01(-i);
    v02(-i);
    v10( i);
    v13(-i);
    v20( i);
    v23(-i);
    v31( i);
    v32( i);
    multiply(-i*0.5);  
  }
};
class HH3 : public a4::sqm< std::complex<double> > {
public:
  HH3(){
    typedef std::complex<double> z;
    z i(0,1);
    v00( 2);
    v11( 0);
    v22( 0);
    v33(-2);
    multiply(-i*0.5);  
  }
};
class HH4 : public a4::sqm< std::complex<double> > {
public:
  HH4(){
    v00( 0);
    v11( 2);
    v22(-2);
    v33( 0);
    multiply(0.5);  
  }
};
class HH5 : public a4::sqm< std::complex<double> > {
public:
  HH5(){
    typedef std::complex<double> z;
    z i(0,1);
    v01( i);
    v02(-i);
    v10(-i);
    v13(-i);
    v20( i);
    v23( i);
    v31( i);
    v32(-i);
    multiply(0.5);  
  }
};
class HH6 : public a4::sqm< std::complex<double> > {
public:
  HH6(){
    typedef std::complex<double> z;
    z i(0,1);
    v01(-1);
    v02( 1);
    v10(-1);
    v13( 1);
    v20( 1);
    v23(-1);
    v31( 1);
    v32(-1);
    multiply(0.5);  
  }
};

}}

#endif