diff --git a/lib/algorithms/iterative/ConjugateGradientMultiShift.h b/lib/algorithms/iterative/ConjugateGradientMultiShift.h
index a9ccfd2c..fca75b06 100644
--- a/lib/algorithms/iterative/ConjugateGradientMultiShift.h
+++ b/lib/algorithms/iterative/ConjugateGradientMultiShift.h
@@ -243,6 +243,8 @@ void operator() (LinearOperatorBase<Field> &Linop, const Field &src, std::vector
       if ( (!converged[s]) ){
 	
 	RealD css  = c * z[s][iz]* z[s][iz];
+    if((k%100)==0 && (s==0) )
+	    std::cout<<GridLogMessage<<"ConjugateGradientMultiShift k="<<k<<" css " <<css <<std::endl;
 	
 	if(css<rsq[s]){
 	  if ( ! converged[s] )