pdeprobleminterface.h

Go to the documentation of this file.

#ifndef PDEPROBLEM_INTERFACE_H_
#define PDEPROBLEM_INTERFACE_H_

#include "dopeexceptionhandler.h"
#include "outputhandler.h"
#include "controlvector.h"
#include "constraintvector.h"
#include "reducedprobleminterface.h"
#include "dwrdatacontainer.h"

#include <assert.h>

#include <deal.II/lac/vector.h>

namespace DOpE
{
  //Predeclaration necessary
  template<typename VECTOR>
    class DOpEOutputHandler;
  template<typename VECTOR>
    class DOpEExceptionHandler;

  template<typename PROBLEM, typename VECTOR, int dealdim>
    class PDEProblemInterface : public ReducedProblemInterface_Base<VECTOR>
    {
      public:
        PDEProblemInterface(PROBLEM *OP, int base_priority = 0)
            : ReducedProblemInterface_Base<VECTOR>()
        {
          OP_ = OP;
          base_priority_ = base_priority;
          post_index_ = "_" + this->GetProblem()->GetName();
        }
        virtual
        ~PDEProblemInterface()
        {
        }

        /******************************************************/

        virtual void
        ReInit()
        {
          this->GetProblem()->ReInit("reduced");
        }

        /******************************************************/

        virtual void
        ComputeReducedFunctionals()=0;

        /******************************************************/

        void
        SetProblemType(std::string type, unsigned int num = 0)
        {
          this->GetProblem()->SetType(type, num);
        }

        template<class DWRC>
          void
          InitializeDWRC(DWRC& dwrc)
          {
            dwrc.Initialize(GetProblem()->GetSpaceTimeHandler(),
                GetProblem()->GetStateNBlocks(),
                GetProblem()->GetStateBlockComponent());
          }

      protected:
        virtual const std::map<std::string, unsigned int>&
        GetFunctionalPosition() const
        {
          return GetProblem()->GetFunctionalPosition();
        }

        PROBLEM*
        GetProblem()
        {
          return OP_;
        }
        const PROBLEM*
        GetProblem() const
        {
          return OP_;
        }
        std::string
        GetPostIndex()
        {
          return post_index_;
        }
        int
        GetBasePriority()
        {
          return base_priority_;
        }
      private:
        PROBLEM* OP_;
        int base_priority_;
        std::string post_index_;
    };

}
#endif