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

Debug of the field strength

Browse Source
This commit is contained in:
Guido Cossu 2017-10-24 10:20:13 +01:00
parent 6579dd30ff
commit cbda4f66e0
2 changed files with 10 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
lib/qcd/utils/WilsonLoops.h
View File

@ -337,7 +337,9 @@ static void StapleMult(GaugeMat &staple, const GaugeLorentz &Umu, int mu) {
GaugeMat v = Vup - Vdn; GaugeMat v = Vup - Vdn;
GaugeMat u = PeekIndex<LorentzIndex>(Umu, mu); // some redundant copies GaugeMat u = PeekIndex<LorentzIndex>(Umu, mu); // some redundant copies
GaugeMat vu = v*u; GaugeMat vu = v*u;
FS = 0.25*Ta(u*v + Cshift(vu, mu, -1)); //FS = 0.25*Ta(u*v + Cshift(vu, mu, -1));
FS = (u*v + Cshift(vu, mu, -1));
FS = 0.125*(FS - adj(FS));
} }
static Real TopologicalCharge(GaugeLorentz &U){ static Real TopologicalCharge(GaugeLorentz &U){

9
tests/qdpxx/Test_qdpxx_wilson.cc
View File

@ -351,6 +351,7 @@ int main(int argc, char **argv)
std::cout << "Norm of difference " << Grid::norm2(difference) << std::endl; std::cout << "Norm of difference " << Grid::norm2(difference) << std::endl;
// Isolate Clover term // Isolate Clover term
/*
calc_grid(Wilson, Ug, src, only_wilson, dag); // Wilson term calc_grid(Wilson, Ug, src, only_wilson, dag); // Wilson term
res_grid -= only_wilson; res_grid -= only_wilson;
res_chroma -= only_wilson; res_chroma -= only_wilson;
@ -359,6 +360,7 @@ int main(int argc, char **argv)
std::cout << "Grid :" << res_grid << std::endl; std::cout << "Grid :" << res_grid << std::endl;
difference = (res_grid-res_chroma); difference = (res_grid-res_chroma);
std::cout << "Difference :" << difference << std::endl; std::cout << "Difference :" << difference << std::endl;
*/
} }
} }
@ -428,7 +430,8 @@ void make_gauge(GaugeField &Umu, FermionField &src)
Grid::QCD::SU3::HotConfiguration(RNG4, Umu); Grid::QCD::SU3::HotConfiguration(RNG4, Umu);
// Fermion field // Fermion field
//Grid::gaussian(RNG4, src); Grid::gaussian(RNG4, src);
/*
Grid::QCD::SpinColourVector F; Grid::QCD::SpinColourVector F;
Grid::Complex c; Grid::Complex c;
@ -450,7 +453,7 @@ void make_gauge(GaugeField &Umu, FermionField &src)
for (int j = 0; j < 3; j++) // colours for (int j = 0; j < 3; j++) // colours
{ {
F()(sp)(j) = Grid::Complex(0.0,0.0); F()(sp)(j) = Grid::Complex(0.0,0.0);
if (((sp == 0)|| (sp==3)) && (j==0)) if (((sp == 0)|| (sp==3)) && (j==2))
{ {
c = Grid::Complex(1.0, 0.0); c = Grid::Complex(1.0, 0.0);
F()(sp)(j) = c; F()(sp)(j) = c;
@ -463,6 +466,8 @@ void make_gauge(GaugeField &Umu, FermionField &src)
} }
} }
} }
*/
} }
void calc_grid(ChromaAction action, Grid::QCD::LatticeGaugeField &Umu, Grid::QCD::LatticeFermion &src, Grid::QCD::LatticeFermion &res, int dag) void calc_grid(ChromaAction action, Grid::QCD::LatticeGaugeField &Umu, Grid::QCD::LatticeFermion &src, Grid::QCD::LatticeFermion &res, int dag)