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

Merge commit '6ad73145bc9754a5f26093eee5a34473ba0cff82' into feature/qed-fvol

Browse Source
This commit is contained in:
Antonin Portelli 2016-12-15 19:48:58 +00:00
commit b10ae00c8a
1 changed files with 31 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

47
programs/qed-fvol/qed-fvol.cc
View File

@ -31,6 +31,9 @@ typedef Photon<QedGimplR> PhotonR;
typedef PhotonR::GaugeField EmField; typedef PhotonR::GaugeField EmField;
typedef PhotonR::GaugeLinkField EmComp; typedef PhotonR::GaugeLinkField EmComp;
const int NCONFIGS = 10;
const int NWILSON = 10;
int main(int argc, char *argv[]) int main(int argc, char *argv[])
{ {
// parse command line // parse command line
@ -64,27 +67,39 @@ int main(int argc, char *argv[])
EmField a(&grid); EmField a(&grid);
EmField expA(&grid); EmField expA(&grid);
Real wlA, logWlA; Complex imag_unit(0, 1);
Real wlA;
std::vector<Real> logWlAvg(NWILSON, 0.0), logWlTime(NWILSON, 0.0), logWlSpace(NWILSON, 0.0);
pRNG.SeedRandomDevice(); pRNG.SeedRandomDevice();
photon.StochasticField(a, pRNG);
// Exponentiate photon field LOG(Message) << "Wilson loop calculation beginning" << std::endl;
Complex imag_unit(0, 1); for(int ic = 0; ic < NCONFIGS; ic++){
expA = exp(imag_unit*a); LOG(Message) << "Configuration " << ic <<std::endl;
photon.StochasticField(a, pRNG);
// Exponentiate photon field
expA = exp(imag_unit*a);
// Calculate Wilson loops
for(int iw=1; iw<=NWILSON; iw++){
wlA = NewWilsonLoops<QedPeriodicGimplR>::avgWilsonLoop(expA, iw, iw) * 3;
logWlAvg[iw-1] -= 2*log(wlA);
wlA = NewWilsonLoops<QedPeriodicGimplR>::avgTimelikeWilsonLoop(expA, iw, iw) * 3;
logWlTime[iw-1] -= 2*log(wlA);
wlA = NewWilsonLoops<QedPeriodicGimplR>::avgSpatialWilsonLoop(expA, iw, iw) * 3;
logWlSpace[iw-1] -= 2*log(wlA);
}
}
LOG(Message) << "Wilson loop calculation completed" << std::endl;
// Calculate Wilson loops // Calculate Wilson loops
for(int i=1; i<=10; i++){ for(int iw=1; iw<=10; iw++){
LOG(Message) << i << 'x' << i << " Wilson loop" << std::endl; LOG(Message) << iw << 'x' << iw << " Wilson loop" << std::endl;
wlA = NewWilsonLoops<QedPeriodicGimplR>::avgWilsonLoop(expA, i, i) * 3; LOG(Message) << "-2log(W) average: " << logWlAvg[iw-1]/NCONFIGS << std::endl;
logWlA = -2*log(wlA); LOG(Message) << "-2log(W) timelike: " << logWlTime[iw-1]/NCONFIGS << std::endl;
LOG(Message) << "-2log(W) average: " << logWlA << std::endl; LOG(Message) << "-2log(W) spatial: " << logWlSpace[iw-1]/NCONFIGS << std::endl;
wlA = NewWilsonLoops<QedPeriodicGimplR>::avgTimelikeWilsonLoop(expA, i, i) * 3;
logWlA = -2*log(wlA);
LOG(Message) << "-2log(W) timelike: " << logWlA << std::endl;
wlA = NewWilsonLoops<QedPeriodicGimplR>::avgSpatialWilsonLoop(expA, i, i) * 3;
logWlA = -2*log(wlA);
LOG(Message) << "-2log(W) spatial: " << logWlA << std::endl;
} }
// epilogue // epilogue