00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013 class Instanton4D {
00014 private:
00015 static int epsilon123(int i, int j, int k) {
00016 if(i==j || j==k || i==k) return 0;
00017 if(i==1 && j==2 && k==3) return 1;
00018 if(i==1 && j==3 && k==2) return -1;
00019 if(i==2 && j==1 && k==3) return -1;
00020 if(i==2 && j==3 && k==1) return 1;
00021 if(i==3 && j==1 && k==2) return 1;
00022 if(i==3 && j==2 && k==1) return -1;
00023 }
00024 public:
00025 vector<mdp_real> p;
00026 int nc;
00027 int sub_i, sub_j;
00028 mdp_real charge;
00029 mdp_real lambda;
00030 mdp_matrix eta[4][4];
00031 Instanton4D(int nc, int sub_i, int sub_j, mdp_real charge, mdp_real lambda, vector<mdp_real> &p) {
00032 this->nc=nc;
00033 this->sub_i=sub_i;
00034 this->sub_j=sub_j;
00035 this->lambda=lambda;
00036 this->charge=charge;
00037 this->p=p;
00038 for(int mu=0; mu<4; mu++)
00039 for(int nu=0; nu<4; nu++) {
00040 this->eta[mu][nu].dimension(nc,nc);
00041 for(int i=0; i<nc; i++)
00042 for(int j=0; j<nc; j++)
00043 if(i==j && (i!=sub_i && i!=sub_j)) {
00044 this->eta[mu][nu](i,j) = 1;
00045 } else if((i!=sub_i && i!=sub_j) || (j!=sub_i && j!=sub_j)) {
00046 this->eta[mu][nu](i,j) = 0;
00047 } else if(mu==0 && nu==0) {
00048 this->eta[mu][nu](i,j) = 0;
00049 } else {
00050 int i0=(i==sub_i?0:1);
00051 int j0=(j==sub_i?0:1);
00052 for(int a=1; a<4; a++) {
00053 if(mu==0) {
00054
00055 this->eta[mu][nu](i,j) += sigma[a](i0,j0)*((a==nu)?-1:0);
00056 } else if(nu==0) {
00057
00058 this->eta[mu][nu](i,j) += sigma[a](i0,j0)*((a==mu)?+1:0);
00059 } else {
00060
00061 this->eta[mu][nu](i,j) += sigma[a](i0,j0)*epsilon123(a,mu,nu);
00062 }
00063 }
00064 }
00065 }
00066 }
00067 mdp_matrix operator()(mdp_site &x, int mu) {
00068 int v[4];
00069 mdp_lattice &lattice=x.lattice();
00070 for(int nu=0; nu<4; nu++) v[nu] = x(nu)-this->p[nu];
00071 float d2 = v[0]*v[0]+v[1]*v[1]+v[2]*v[2]+v[3]*v[3];
00072 mdp_matrix A(this->nc,this->nc);
00073 A=0;
00074 for(int nu=0; nu<4; nu++)
00075 A += this->eta[mu][nu]*v[nu];
00076 return (2.0*this->charge/(d2+this->lambda*this->lambda))*A;
00077 }
00078 };
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097
00098
00099
00100
00101
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114