From 9394450c1a5da139ccd55295893d6f333cb99db9 Mon Sep 17 00:00:00 2001
From: Quadro <Peter Boyle>
Date: Wed, 9 Jun 2021 13:30:42 -0400
Subject: [PATCH] Verbose changes

---
 Grid/qcd/utils/MixedPrecisionOperatorFunction.h | 11 ++++++-----
 1 file changed, 6 insertions(+), 5 deletions(-)

diff --git a/Grid/qcd/utils/MixedPrecisionOperatorFunction.h b/Grid/qcd/utils/MixedPrecisionOperatorFunction.h
index 7ef26c85..a8eefaba 100644
--- a/Grid/qcd/utils/MixedPrecisionOperatorFunction.h
+++ b/Grid/qcd/utils/MixedPrecisionOperatorFunction.h
@@ -56,7 +56,7 @@ class MixedPrecisionConjugateGradientOperatorFunction : public OperatorFunction<
     Integer TotalOuterIterations; //Number of restarts
     Integer TotalFinalStepIterations; //Number of CG iterations in final patch-up step
 
-    MixedPrecisionConjugateGradientOperatorFunction(RealD tol, 
+     MixedPrecisionConjugateGradientOperatorFunction(RealD tol, RealD tolInner,
 						    Integer maxinnerit, 
 						    Integer maxouterit,
 						    GridBase *_SinglePrecGrid,
@@ -69,12 +69,12 @@ class MixedPrecisionConjugateGradientOperatorFunction : public OperatorFunction<
       FermOpF(_FermOpF),
       FermOpD(_FermOpD),
       Tolerance(tol), 
-      InnerTolerance(tol), 
+      InnerTolerance(tolInner), 
       MaxInnerIterations(maxinnerit), 
       MaxOuterIterations(maxouterit), 
       SinglePrecGrid(_SinglePrecGrid),
       OuterLoopNormMult(100.) 
-    {     };
+  { assert(tolInner<0.01);    };
 
     void operator()(LinearOperatorBase<FieldD> &LinOpU, const FieldD &src, FieldD &psi)
     {
@@ -101,8 +101,9 @@ class MixedPrecisionConjugateGradientOperatorFunction : public OperatorFunction<
       // Make a mixed precision conjugate gradient
       //////////////////////////////////////////////////////////////////////////////////////////
       // Could assume red black solver here and remove the SinglePrecGrid parameter???
-      MixedPrecisionConjugateGradient<FieldD,FieldF> MPCG(Tolerance,MaxInnerIterations,MaxOuterIterations,SinglePrecGrid,LinOpF,LinOpD);
-      std::cout << GridLogMessage << "Calling mixed precision Conjugate Gradient src "<<norm2(src)<< " U "<<norm2(Umu_d) <<std::endl;
+      MixedPrecisionConjugateGradient<FieldD,FieldF> MPCG(Tolerance, InnerTolerance,MaxInnerIterations,MaxOuterIterations,SinglePrecGrid,LinOpF,LinOpD);
+
+      std::cout << GridLogMessage << "Calling mixed precision Conjugate Gradient src "<<norm2(src) <<std::endl;
       psi=Zero();
       MPCG(src,psi);
     }