diff --git a/dune/solvers/iterationsteps/projectedgradientstep.hh b/dune/solvers/iterationsteps/projectedgradientstep.hh
index 1d80800b544cc994fd21dc016bbb5f55431d1b4b..cdea8d7ec3d716f9d1c72fc439c0940a20466610 100644
--- a/dune/solvers/iterationsteps/projectedgradientstep.hh
+++ b/dune/solvers/iterationsteps/projectedgradientstep.hh
@@ -52,6 +52,9 @@ public:
         computeGeneralizedCP(negativeGradient);
     }
 
+    //! Return solution object
+    virtual VectorType getSol() {return *this->x_;}
+
     //! The obstacles
     const std::vector<BoxConstraint<field_type,blocksize> >* obstacles_;
 private: