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

uninitialised bug fix

Browse Source
This commit is contained in:
Azusa Yamaguchi 2015-06-16 14:07:05 +01:00
parent 73494a4768
commit 4e7300b68d
1 changed files with 6 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -30,6 +30,7 @@ public:
static void sitePlaquette(LatticeComplex &Plaq,const std::vector<GaugeMat> &U) static void sitePlaquette(LatticeComplex &Plaq,const std::vector<GaugeMat> &U)
{ {
LatticeComplex sitePlaq(U[0]._grid); LatticeComplex sitePlaq(U[0]._grid);
Plaq=zero;
for(int mu=1;mu<Nd;mu++){ for(int mu=1;mu<Nd;mu++){
for(int nu=0;nu<mu;nu++){ for(int nu=0;nu<mu;nu++){
traceDirPlaquette(sitePlaq,U,mu,nu); traceDirPlaquette(sitePlaq,U,mu,nu);
@ -42,6 +43,7 @@ public:
////////////////////////////////////////////////// //////////////////////////////////////////////////
static RealD sumPlaquette(const GaugeLorentz &Umu){ static RealD sumPlaquette(const GaugeLorentz &Umu){
std::vector<GaugeMat> U(4,Umu._grid); std::vector<GaugeMat> U(4,Umu._grid);
for(int mu=0;mu<Nd;mu++){ for(int mu=0;mu<Nd;mu++){
U[mu] = peekIndex<LorentzIndex>(Umu,mu); U[mu] = peekIndex<LorentzIndex>(Umu,mu);
} }
@ -72,13 +74,15 @@ public:
////////////////////////////////////////////////// //////////////////////////////////////////////////
static void Staple(GaugeMat &staple,const GaugeLorentz &Umu,int mu){ static void Staple(GaugeMat &staple,const GaugeLorentz &Umu,int mu){
std::vector<GaugeMat> U(4,Umu._grid); GridBase *grid = Umu._grid;
std::vector<GaugeMat> U(4,grid);
for(int d=0;d<Nd;d++){ for(int d=0;d<Nd;d++){
U[d] = peekIndex<LorentzIndex>(Umu,d); U[d] = peekIndex<LorentzIndex>(Umu,d);
} }
staple = zero; staple = zero;
GaugeMat tmp(U[0]._grid); GaugeMat tmp(grid);
for(int nu=0;nu<Nd;nu++){ for(int nu=0;nu<Nd;nu++){