portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-06-19 16:27:05 +01:00
Code Releases Activity

Updates to ldop tests

Browse Source
This commit is contained in:
Peter Boyle
2015-06-10 12:26:25 +01:00
parent 622261b1ea
commit 5cce44edb4
5 changed files with 137 additions and 22 deletions
Download Patch File Download Diff File Expand all files Collapse all files

17
tests/Test_cayley_ldop_cr.cc
View File

@ -64,21 +64,10 @@ int main (int argc, char ** argv)
//random(RNG4,Umu);
NerscField header;
std::string file("./ckpoint_lat.400");
std::string file("./ckpoint_lat.4000");
readNerscConfiguration(Umu,header,file);
#if 0
LatticeColourMatrix U(UGrid);
Complex cone(1.0,0.0);
for(int nn=0;nn<Nd;nn++){
if (nn==3) {
U=zero; std::cout << "zeroing gauge field in dir "<<nn<<std::endl;
// else { U[nn]= cone;std::cout << "unit gauge field in dir "<<nn<<std::endl; }
pokeIndex<LorentzIndex>(Umu,U,nn);
}
}
#endif
RealD mass=0.5;
RealD mass=0.01;
RealD M5=1.8;
DomainWallFermion Ddwf(Umu,*FGrid,*FrbGrid,*UGrid,*UrbGrid,mass,M5);
@ -100,7 +89,7 @@ int main (int argc, char ** argv)
// filter(HermIndefOp,noise,subspace[b]);
// inverse iteration
MdagMLinearOperator<DomainWallFermion,LatticeFermion> HermDefOp(Ddwf);
ConjugateGradient<LatticeFermion> CG(1.0e-5,10000);
ConjugateGradient<LatticeFermion> CG(1.0e-4,10000);
for(int i=0;i<4;i++){