portelli/Grid
1
0 Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2026-10-07 08:18:06 +01:00
Code Releases Activity

HMC for Wilson Gauge action works

Browse Source 
Fixed bug in momenta generation
This commit is contained in:
neo committed 2015-07-06 12:58:49 +09:00
1 parent 6261770f59
commit 12e1682a87
9 files changed
+24 -40
Download Patch File Download Diff File Expand all files Collapse all files

No files matched your search

lib/qcd/action/gauge/WilsonGaugeAction.h
-3
View File
@@ -27,16 +27,13 @@ namespace Grid{
//not optimal implementation FIXME
//extend Ta to include Lorentz indexes
RealD factor = 0.5*beta/RealD(Nc);
std::cout << "Debug force Wilson "<< factor << "\n";
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);
std::cout << "Staple norm ["<<mu<<"] = "<< norm2(dSdU_mu) <<"\n";
dSdU_mu = Ta(Umu*adj(dSdU_mu))*factor;
std::cout << "Force norm ["<<mu<<"] = "<< norm2(dSdU_mu) << " : Umu "<< norm2(Umu)<<"\n";
pokeLorentz(dSdU, dSdU_mu, mu);
}
};