diff --git a/tests/lanczos/Test_wilson_DWFKernel.cc b/tests/lanczos/Test_wilson_DWFKernel.cc index ffd01b52..ab67e251 100644 --- a/tests/lanczos/Test_wilson_DWFKernel.cc +++ b/tests/lanczos/Test_wilson_DWFKernel.cc @@ -112,6 +112,7 @@ public: struct LanczosParameters: Serializable { GRID_SERIALIZABLE_CLASS_MEMBERS(LanczosParameters, RealD, mass , + RealD, resid, RealD, ChebyLow, RealD, ChebyHigh, Integer, ChebyOrder) @@ -200,9 +201,9 @@ int main(int argc, char** argv) { } */ - int Nstop = 10; - int Nk = 20; - int Np = 80; + int Nstop = 5; + int Nk = 10; + int Np = 90; int Nm = Nk + Np; int MaxIt = 10000; RealD resid = 1.0e-5; @@ -227,6 +228,7 @@ int main(int argc, char** argv) { } mass=LanParams.mass; + resid=LanParams.resid; while ( mass > - 5.0){ @@ -235,16 +237,16 @@ while ( mass > - 5.0){ //SchurDiagTwoOperator HermOp(WilsonOperator); // Gamma5HermitianLinearOperator HermOp2(WilsonOperator); /// <----- - std::vector Coeffs{0, 1.}; + std::vector Coeffs{0, 0, 1.}; Polynomial PolyX(Coeffs); -// Chebyshev Cheby(LanParams.ChebyLow,LanParams.ChebyHigh,LanParams.ChebyOrder); + Chebyshev Cheby(LanParams.ChebyLow,LanParams.ChebyHigh,LanParams.ChebyOrder); -// FunctionHermOp OpCheby(Cheby,HermOp); + FunctionHermOp OpCheby(Cheby,HermOp); // InvHermOp Op(WilsonOperator,HermOp); PlainHermOp Op (HermOp); // PlainHermOp Op2 (HermOp2); - ImplicitlyRestartedLanczos IRL(Op, Op, Nstop, Nk, Nm, resid, MaxIt); + ImplicitlyRestartedLanczos IRL(OpCheby, Op, Nstop, Nk, Nm, resid, MaxIt); std::vector eval(Nm); FermionField src(FGrid);