// Defines Op class -- transforms vectors to vectors
#ifndef OPERATOR_H
#define OPERATOR_H
template <class Vector> struct Op { // returns y = A*x
  virtual Vector &operator()(const Vector &x, Vector &y) = 0;
};
#endif // OPERATOR_H