portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-09-20 09:15:38 +01:00
Code Releases Activity

Gparity now accepting twists through constructor

Browse Source
This commit is contained in:
Peter Boyle 2015-08-19 11:26:01 +01:00
parent 4e085dd0ed
commit 8b070ae54c
1 changed files with 7 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
tests/Test_gparity.cc
View File

@ -13,6 +13,9 @@ int main (int argc, char ** argv)
Grid_init(&argc,&argv);
std::vector<int> twists(Nd,0);
twists[nu] = 1;
const int Ls=4;
const int L =4;
std::vector<int> latt_2f(Nd,L);
@ -95,8 +98,10 @@ int main (int argc, char ** argv)
SchurDiagMooeeOperator<DomainWallFermionR,LatticeFermion> HermOpEO(Ddwf);
ConjugateGradient<LatticeFermion> CG(1.0e-8,10000);
CG(HermOpEO,src_o_1f,result_o_1f);
GparityDomainWallFermionR GPDdwf(Umu_2f,*FGrid_2f,*FrbGrid_2f,*UGrid_2f,*UrbGrid_2f,mass,M5);
GparityDomainWallFermionR::ImplParams params;
params.twists = twists;
GparityDomainWallFermionR GPDdwf(Umu_2f,*FGrid_2f,*FrbGrid_2f,*UGrid_2f,*UrbGrid_2f,mass,M5,params);
for(int disp=-1;disp<=1;disp+=2)
for(int mu=0;mu<5;mu++)