iterationstep.hh

#ifndef ITERATIONSTEP_HH
#define ITERATIONSTEP_HH

#include <dune-solvers/numproc.hh>
#include <vector>
#include <string>

#include <dune/common/bitsetvector.hh>
#include <dune-solvers/common/canignore.hh>

    //! Base class for iteration steps being called by an iterative solver
template<class VectorType, class BitVectorType = Dune::BitSetVector<VectorType::block_type::dimension> >
class IterationStep : public NumProc, public CanIgnore<BitVectorType>
    {
    public:

        //! Default constructor
        IterationStep()
        {}

        /** \brief Destructor */
        virtual ~IterationStep() {}

        //! Constructor being given the current iterate
        IterationStep(VectorType& x) :
            x_(&x)
        {}

        //! Set the current iterate
        virtual void setProblem(VectorType& x) {
            x_ = &x;
        }

        //! to be called before iteration
        virtual void preprocess() {}

        //! Do the actual iteration
        virtual void iterate() = 0;

        //! Return solution object
        virtual VectorType getSol() = 0;

        /** \brief Checks whether all relevant member variables are set
         * \exception SolverError if the iteration step is not set up properly
         */
        virtual void check() const {
#if 0
            if (!x_)
                DUNE_THROW(SolverError, "Iteration step has no solution vector");
#endif
        }

        //! get output of last iteration step
        virtual std::string getOutput() const
        {
            return "";
        }

        //! The solution container
        VectorType* x_;


        using CanIgnore<BitVectorType>::ignoreNodes_;

    };

#endif