diff --git a/lib/algorithms/LinearOperator.h b/lib/algorithms/LinearOperator.h
index 0d32cc15..2a68a7b9 100644
--- a/lib/algorithms/LinearOperator.h
+++ b/lib/algorithms/LinearOperator.h
@@ -346,6 +346,13 @@ namespace Grid {
       virtual void operator() (const Field &in, Field &out) = 0;
     };
 
+    template<class Field> class IdentityLinearFunction : public LinearFunction<Field> {
+    public:
+      void operator() (const Field &in, Field &out){
+	out = in;
+      };
+    };
+
 
     /////////////////////////////////////////////////////////////
     // Base classes for Multishift solvers for operators