implement method setProblem()

[[Imported from SVN: r2952]]

dune-solvers/iterationsteps/amgstep.hh

@@ -52,8 +52,16 @@ public:
              VectorType& rhs)
     {
         setProblem(*stiffnessMatrix, x, rhs);
+    }
 
-        fop_ = std::auto_ptr<Operator>(new Operator(*stiffnessMatrix));
+    /** \brief Initialize the iteration step but don't actually build the matrix hierarchy yet */
+    virtual void setProblem(const MatrixType& stiffnessMatrix, 
+             VectorType& x, 
+             VectorType& rhs)
+    {
+        LinearIterationStep<MatrixType, VectorType>::setProblem(stiffnessMatrix, x, rhs);
+
+        fop_ = std::auto_ptr<Operator>(new Operator(stiffnessMatrix));
     }
 
     /** \brief Sets up an algebraic hierarchy