1
0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 15:55:37 +00:00
Grid/lib/qcd/action/gauge/WilsonGaugeAction.h

46 lines
1.3 KiB
C
Raw Normal View History

2015-06-05 10:02:36 +01:00
#ifndef QCD_WILSON_GAUGE_ACTION_H
#define QCD_WILSON_GAUGE_ACTION_H
2015-07-03 18:43:14 +01:00
namespace Grid{
namespace QCD{
////////////////////////////////////////////////////////////////////////
// Wilson Gauge Action .. should I template the Nc etc..
////////////////////////////////////////////////////////////////////////
template<class GaugeField,class MatrixField>
class WilsonGaugeAction : public Action<GaugeField> {
private:
RealD beta;
public:
WilsonGaugeAction(RealD b):beta(b){};
2015-07-04 09:47:50 +01:00
virtual void init(const GaugeField &U, GridParallelRNG& pRNG) {};
2015-07-03 18:43:14 +01:00
virtual RealD S(const GaugeField &U) {
RealD plaq = WilsonLoops<MatrixField,GaugeField>::avgPlaquette(U);
std::cout << "Plaq : "<<plaq << "\n";
double vol = U._grid->gSites();
return beta*(1.0 -plaq)*(Nd*(Nd-1.0))*vol*0.5;
2015-07-03 18:43:14 +01:00
};
virtual void deriv(const GaugeField &U,GaugeField & dSdU) {
//not optimal implementation FIXME
//extend Ta to include Lorentz indexes
RealD factor = 0.5*beta/RealD(Nc);
2015-07-03 18:43:14 +01:00
MatrixField dSdU_mu(U._grid);
MatrixField Umu(U._grid);
for (int mu=0; mu < Nd; mu++){
Umu = PeekIndex<LorentzIndex>(U,mu);
// Staple in direction mu
WilsonLoops<MatrixField,GaugeField>::Staple(dSdU_mu,U,mu);
dSdU_mu = Ta(Umu*adj(dSdU_mu))*factor;
pokeLorentz(dSdU, dSdU_mu, mu);
}
2015-07-03 18:43:14 +01:00
};
};
}
}
2015-06-05 10:02:36 +01:00
#endif