diff --git a/lib/algorithms/LinearOperator.h b/lib/algorithms/LinearOperator.h
new file mode 100644
index 00000000..5929570f
--- /dev/null
+++ b/lib/algorithms/LinearOperator.h
@@ -0,0 +1,18 @@
+#ifndef  GRID_ALGORITHM_LINEAR_OP_H
+#define  GRID_ALGORITHM_LINEAR_OP_H
+
+#include <Grid.h>
+
+namespace Grid {
+
+  // Red black cases?
+    template<class vtype> class LinearOperatorBase {
+    public:
+      void M(const Lattice<vtype> &in, Lattice<vtype> &out){ assert(0);}
+      void Mdag(const Lattice<vtype> &in, Lattice<vtype> &out){ assert(0);}
+      void MdagM(const Lattice<vtype> &in, Lattice<vtype> &out){ assert(0);}
+    };
+
+}
+
+#endif