diff --git a/lib/algorithms/iterative/ConjugateGradientMultiShift.h b/lib/algorithms/iterative/ConjugateGradientMultiShift.h
index 9d2719fc..4c311bc1 100644
--- a/lib/algorithms/iterative/ConjugateGradientMultiShift.h
+++ b/lib/algorithms/iterative/ConjugateGradientMultiShift.h
@@ -208,6 +208,7 @@ void operator() (LinearOperatorBase<Field> &Linop, const Field &src, std::vector
     cp=c;
     MatrixTimer.Start();  
     //Linop.HermOpAndNorm(p,mmp,d,qq); // d is used
+    // The below is faster on KNL
     Linop.HermOp(p,mmp); 
     d=real(innerProduct(p,mmp));