portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

Fix the Charge conjugate BC's

Browse Source
This commit is contained in:
paboyle 2018-03-05 13:59:02 +00:00
parent 4d53703c67
commit e199ba7e88
2 changed files with 15 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
lib/qcd/action/gauge/WilsonGaugeAction.h
View File

@ -73,17 +73,19 @@ public:
//GaugeLinkField Umu(U.Grid()); //GaugeLinkField Umu(U.Grid());
GaugeLinkField dSdU_mu(U.Grid()); GaugeLinkField dSdU_mu(U.Grid());
for (int mu = 0; mu < Nd; mu++) { for (int mu = 0; mu < Nd; mu++) {
//Umu = PeekIndex<LorentzIndex>(U, mu);
// Staple in direction mu GaugeLinkField Umu(U.Grid());
//WilsonLoops<Gimpl>::Staple(dSdU_mu, U, mu); GaugeLinkField dSdU_mu(U.Grid());
//dSdU_mu = Ta(Umu * dSdU_mu) * factor; for (int mu = 0; mu < Nd; mu++) {
Umu = PeekIndex<LorentzIndex>(U, mu);
// Staple in direction mu
WilsonLoops<Gimpl>::Staple(dSdU_mu, U, mu);
dSdU_mu = Ta(Umu * dSdU_mu) * factor;
PokeIndex<LorentzIndex>(dSdU, dSdU_mu, mu);
}
WilsonLoops<Gimpl>::StapleMult(dSdU_mu, U, mu);
dSdU_mu = Ta(dSdU_mu) * factor;
PokeIndex<LorentzIndex>(dSdU, dSdU_mu, mu);
} }
} }
private: private:

8
tests/forces/Test_gp_rect_force.cc
View File

@ -57,9 +57,9 @@ int main (int argc, char ** argv)
SU3::HotConfiguration(pRNG,U); SU3::HotConfiguration(pRNG,U);
double beta = 1.0; double beta = 1.0;
double c1 = -1.0/12.0;
//GparityPlaqPlusRectangleActionR Action(beta,c1); ConjugatePlaqPlusRectangleActionR Action(beta,c1);
ConjugateWilsonGaugeActionR Action(beta); //ConjugateWilsonGaugeActionR Action(beta);
//WilsonGaugeActionR Action(beta); //WilsonGaugeActionR Action(beta);
ComplexD S = Action.S(U); ComplexD S = Action.S(U);
@ -72,7 +72,7 @@ int main (int argc, char ** argv)
//////////////////////////////////// ////////////////////////////////////
// Modify the gauge field a little // Modify the gauge field a little
//////////////////////////////////// ////////////////////////////////////
RealD dt = 0.0001; RealD dt = 0.001;
LatticeColourMatrix mommu(&Grid); LatticeColourMatrix mommu(&Grid);
LatticeColourMatrix forcemu(&Grid); LatticeColourMatrix forcemu(&Grid);