portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-26 05:35:55 +01:00
Code Releases Activity

Reverting whitespace changes

Browse Source
This commit is contained in:
Mashy Green 2025-02-17 17:16:28 +00:00
parent e99d42404e
commit 7726ee4b16
2 changed files with 20 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
Grid/qcd/action/gauge/PlaqPlusRectangleAction.h
View File

@ -72,29 +72,29 @@ public:
}; };
virtual void deriv(const GaugeField &U, GaugeField &dSdU) { virtual void deriv(const GaugeField &U, GaugeField &dSdU) {
// extend Ta to include Lorentz indexes //extend Ta to include Lorentz indexes
RealD factor_p = c_plaq / RealD(Nc) * 0.5; RealD factor_p = c_plaq/RealD(Nc)*0.5;
RealD factor_r = c_rect / RealD(Nc) * 0.5; RealD factor_r = c_rect/RealD(Nc)*0.5;
GridBase *grid = U.Grid(); GridBase *grid = Umu.Grid();
std::vector<GaugeLinkField> Umu(Nd, grid); std::vector<GaugeLinkField> U (Nd,grid);
for (int mu = 0; mu < Nd; mu++) { for(int mu=0;mu<Nd;mu++){
Umu[mu] = PeekIndex<LorentzIndex>(U, mu); U[mu] = PeekIndex<LorentzIndex>(Umu,mu);
} }
std::vector<GaugeLinkField> RectStaple(Nd, grid), Staple(Nd, grid); std::vector<GaugeLinkField> RectStaple(Nd,grid), Staple(Nd,grid);
WilsonLoops<Gimpl>::StapleAndRectStapleAll(Staple, RectStaple, Umu, WilsonLoops<Gimpl>::StapleAndRectStapleAll(Staple, RectStaple, U, workspace);
workspace);
GaugeLinkField dSdU_mu(grid); GaugeLinkField dSdU_mu(grid);
GaugeLinkField staple(grid); GaugeLinkField staple(grid);
for (int mu = 0; mu < Nd; mu++) { for (int mu=0; mu < Nd; mu++){
dSdU_mu = Ta(Umu[mu] * Staple[mu]) * factor_p; dSdU_mu = Ta(U[mu]*Staple[mu])*factor_p;
dSdU_mu = dSdU_mu + Ta(Umu[mu] * RectStaple[mu]) * factor_r; dSdU_mu = dSdU_mu + Ta(U[mu]*RectStaple[mu])*factor_r;
PokeIndex<LorentzIndex>(dSdU, dSdU_mu, mu); PokeIndex<LorentzIndex>(dSdU, dSdU_mu, mu);
} }
}; };
}; };

11
Grid/qcd/utils/WilsonLoops.h
View File

@ -318,8 +318,11 @@ public:
// //
staple += Gimpl::ShiftStaple( staple += Gimpl::ShiftStaple(
Gimpl::CovShiftForward(Umu[nu], nu, Gimpl::CovShiftForward(
Gimpl::CovShiftBackward(Umu[mu], mu, Gimpl::CovShiftIdentityBackward(Umu[nu], nu))), mu); U[nu], nu,
Gimpl::CovShiftBackward(
U[mu], mu, Gimpl::CovShiftIdentityBackward(U[nu], nu))),
mu);
// __ // __
// | // |
@ -328,8 +331,8 @@ public:
// //
staple += Gimpl::ShiftStaple( staple += Gimpl::ShiftStaple(
Gimpl::CovShiftBackward(Umu[nu], nu, Gimpl::CovShiftBackward(U[nu], nu,
Gimpl::CovShiftBackward(Umu[mu], mu, Umu[nu])), mu); Gimpl::CovShiftBackward(U[mu], mu, U[nu])), mu);
} }
} }
} }