mirror of
https://github.com/paboyle/Grid.git
synced 2024-11-14 09:45:36 +00:00
91 lines
2.4 KiB
Plaintext
91 lines
2.4 KiB
Plaintext
|
Branch: develop
|
||
|
|
|
||
|
|
Files:
|
||
|
|
|
||
|
|
Grid/lattice/PaddedCell.h -- Halo exchange
|
||
|
|
tests/Test_general_stencil.cc -- test local off axis stencil addressing
|
||
|
|
tests/debug/Test_padded_cell.cc -- test PaddedCell halo exchange and the General local stencil by computing ALL plaquettes on lattice
|
||
|
|
|
||
|
|
Functionality:
|
||
|
|
|
||
|
|
-- extend a lattice field:
|
||
|
|
Grid/lattice/PaddedCell.h
|
||
|
|
|
||
|
|
// Constructor
|
||
|
|
PaddedCell(int _depth,GridCartesian *_grid)
|
||
|
|
|
||
|
|
// Expand a field "in" to depth "d"
|
||
|
|
template<class vobj>
|
||
|
|
inline Lattice<vobj> Exchange(Lattice<vobj> &in)
|
||
|
|
|
||
|
|
// Take the "apple core" of in to a smaller local volume
|
||
|
|
template<class vobj>
|
||
|
|
inline Lattice<vobj> Extract(Lattice<vobj> &in)
|
||
|
|
|
||
|
|
-- Plaquette test:
|
||
|
|
tests/debug/Test_padded_cell.cc
|
||
|
|
/////////////////////////////////////////////////
|
||
|
|
// Create a padded cell of extra padding depth=1
|
||
|
|
/////////////////////////////////////////////////
|
||
|
|
int depth = 1;
|
||
|
|
PaddedCell Ghost(depth,&GRID);
|
||
|
|
LatticeGaugeField Ughost = Ghost.Exchange(Umu);
|
||
|
|
|
||
|
|
///// Array for the site plaquette
|
||
|
|
GridBase *GhostGrid = Ughost.Grid();
|
||
|
|
LatticeComplex gplaq(GhostGrid);
|
||
|
|
|
||
|
|
std::vector<Coordinate> shifts;
|
||
|
|
for(int mu=0;mu<Nd;mu++){
|
||
|
|
for(int nu=mu+1;nu<Nd;nu++){
|
||
|
|
|
||
|
|
// Umu(x) Unu(x+mu) Umu^dag(x+nu) Unu^dag(x)
|
||
|
|
Coordinate shift_0(Nd,0);
|
||
|
|
Coordinate shift_mu(Nd,0); shift_mu[mu]=1;
|
||
|
|
Coordinate shift_nu(Nd,0); shift_nu[nu]=1;
|
||
|
|
shifts.push_back(shift_0);
|
||
|
|
shifts.push_back(shift_mu);
|
||
|
|
shifts.push_back(shift_nu);
|
||
|
|
shifts.push_back(shift_0);
|
||
|
|
}
|
||
|
|
}
|
||
|
|
GeneralLocalStencil gStencil(GhostGrid,shifts);
|
||
|
|
|
||
|
|
gplaq=Zero();
|
||
|
|
{
|
||
|
|
autoView( gp_v , gplaq, CpuWrite);
|
||
|
|
autoView( t_v , trplaq, CpuRead);
|
||
|
|
autoView( U_v , Ughost, CpuRead);
|
||
|
|
for(int ss=0;ss<gp_v.size();ss++){
|
||
|
|
int s=0;
|
||
|
|
for(int mu=0;mu<Nd;mu++){
|
||
|
|
for(int nu=mu+1;nu<Nd;nu++){
|
||
|
|
|
||
|
|
auto SE0 = gStencil.GetEntry(s+0,ss);
|
||
|
|
auto SE1 = gStencil.GetEntry(s+1,ss);
|
||
|
|
auto SE2 = gStencil.GetEntry(s+2,ss);
|
||
|
|
auto SE3 = gStencil.GetEntry(s+3,ss);
|
||
|
|
|
||
|
|
int o0 = SE0->_offset;
|
||
|
|
int o1 = SE1->_offset;
|
||
|
|
int o2 = SE2->_offset;
|
||
|
|
int o3 = SE3->_offset;
|
||
|
|
|
||
|
|
auto U0 = U_v[o0](mu);
|
||
|
|
auto U1 = U_v[o1](nu);
|
||
|
|
auto U2 = adj(U_v[o2](mu));
|
||
|
|
auto U3 = adj(U_v[o3](nu));
|
||
|
|
|
||
|
|
gpermute(U0,SE0->_permute);
|
||
|
|
gpermute(U1,SE1->_permute);
|
||
|
|
gpermute(U2,SE2->_permute);
|
||
|
|
gpermute(U3,SE3->_permute);
|
||
|
|
|
||
|
|
gp_v[ss]() =gp_v[ss]() + trace( U0*U1*U2*U3 );
|
||
|
|
s=s+4;
|
||
|
|
}
|
||
|
|
}
|
||
|
|
}
|
||
|
|
}
|
||
|
|
cplaq = Ghost.Extract(gplaq);
|