diff --git a/tests/debug/Test_general_coarse_hdcg.cc b/tests/debug/Test_general_coarse_hdcg.cc index b6508ecb..f8587b78 100644 --- a/tests/debug/Test_general_coarse_hdcg.cc +++ b/tests/debug/Test_general_coarse_hdcg.cc @@ -77,13 +77,13 @@ int main (int argc, char ** argv) { Grid_init(&argc,&argv); - const int Ls=24; - const int nbasis = 60; + const int Ls=16; + const int nbasis = 40; const int cb = 0 ; - RealD mass=0.00078; + RealD mass=0.01; RealD M5=1.8; - RealD b=1.5; - RealD c=0.5; + RealD b=1.0; + RealD c=0.0; GridCartesian * UGrid = SpaceTimeGrid::makeFourDimGrid(GridDefaultLatt(), GridDefaultSimd(Nd,vComplex::Nsimd()), @@ -167,7 +167,7 @@ int main (int argc, char ** argv) std::cout << "**************************************"< coarseCG(4.0e-2,20000,true); - const int nrhs=12; + const int nrhs=8; Coordinate mpi=GridDefaultMpi(); Coordinate rhMpi ({1,1,mpi[0],mpi[1],mpi[2],mpi[3]}); @@ -185,7 +185,7 @@ int main (int argc, char ** argv) std::cout << "**************************************"< MrhsHermMatrix; - Chebyshev IRLCheby(0.01,42.0,301); // 1 iter + Chebyshev IRLCheby(0.5,40.0,101); // 1 iter MrhsHermMatrix MrhsCoarseOp (mrhs); CoarseVector pm_src(CoarseMrhs); @@ -193,8 +193,8 @@ int main (int argc, char ** argv) PowerMethod cPM; cPM(MrhsCoarseOp,pm_src); - int Nk=192; - int Nm=384; + int Nk=nrhs; + int Nm=Nk*3; int Nstop=Nk; int Nconv_test_interval=1; @@ -250,9 +250,9 @@ int main (int argc, char ** argv) // Extra HDCG parameters ////////////////////////// int maxit=3000; - ConjugateGradient CG(5.0e-2,maxit,false); + ConjugateGradient CG(2.0e-1,maxit,false); RealD lo=2.0; - int ord = 7; + int ord = 9; DoNothingGuesser DoNothing; HPDSolver HPDSolveMrhs(MrhsCoarseOp,CG,DoNothing);