diff --git a/tests/solver/Test_dwf_mr_unprec.cc b/tests/solver/Test_dwf_mr_unprec.cc deleted file mode 100644 index 3ca35e32..00000000 --- a/tests/solver/Test_dwf_mr_unprec.cc +++ /dev/null @@ -1,68 +0,0 @@ -/************************************************************************************* - -Grid physics library, www.github.com/paboyle/Grid - -Source file: ./tests/solver/Test_dwf_mr_unprec.cc - -Copyright (C) 2015 - -Author: Daniel Richtmann - -This program is free software; you can redistribute it and/or modify -it under the terms of the GNU General Public License as published by -the Free Software Foundation; either version 2 of the License, or -(at your option) any later version. - -This program is distributed in the hope that it will be useful, -but WITHOUT ANY WARRANTY; without even the implied warranty of -MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the -GNU General Public License for more details. - -You should have received a copy of the GNU General Public License along -with this program; if not, write to the Free Software Foundation, Inc., -51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. - -See the full license in the file "LICENSE" in the top level distribution -directory -*************************************************************************************/ -/* END LEGAL */ -#include - -using namespace Grid; -using namespace Grid::QCD; - -int main (int argc, char ** argv) -{ - Grid_init(&argc,&argv); - - const int Ls=8; - - GridCartesian * UGrid = SpaceTimeGrid::makeFourDimGrid(GridDefaultLatt(), GridDefaultSimd(Nd,vComplex::Nsimd()),GridDefaultMpi()); - GridRedBlackCartesian * UrbGrid = SpaceTimeGrid::makeFourDimRedBlackGrid(UGrid); - GridCartesian * FGrid = SpaceTimeGrid::makeFiveDimGrid(Ls,UGrid); - GridRedBlackCartesian * FrbGrid = SpaceTimeGrid::makeFiveDimRedBlackGrid(Ls,UGrid); - - std::vector seeds4({1,2,3,4}); - std::vector seeds5({5,6,7,8}); - GridParallelRNG RNG5(FGrid); RNG5.SeedFixedIntegers(seeds5); - GridParallelRNG RNG4(UGrid); RNG4.SeedFixedIntegers(seeds4); - - LatticeFermion src(FGrid); random(RNG5,src); - LatticeFermion result(FGrid); result=zero; - LatticeGaugeField Umu(UGrid); SU3::HotConfiguration(RNG4,Umu); - - std::vector U(4,UGrid); - for(int mu=0;mu(Umu,mu); - } - - RealD mass=0.1; - RealD M5=1.8; - DomainWallFermionR Ddwf(Umu,*FGrid,*FrbGrid,*UGrid,*UrbGrid,mass,M5); - - MdagMLinearOperator HermOp(Ddwf); - MinimalResidual MR(1.0e-6,10000,0.8); - MR(HermOp,src,result); - - Grid_finalize(); -}