ts_base.h

Go to the documentation of this file.

#ifndef TSBase_H_
#define TSBase_H_

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

namespace DOpE
{
  template<typename OPTPROBLEM, typename SPARSITYPATTERN, typename VECTOR,
      int dopedim, int dealdim,
      template <int, int> class FE = dealii::FESystem,
      template <int, int> class DH = dealii::DoFHandler>
    class TSBase
    {
      public:
        TSBase(OPTPROBLEM& OP) :
            OP_(OP)
        {
        }
        ;
        ~TSBase()
        {
        }
        ;

        /******************************************************/
        void
        SetStepPart(std::string s)
        {
          part_ = s;
        }

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

        void
        SetTime(double time,             
                unsigned int time_dof_number,
                const TimeIterator& interval, bool initial = false)
        {
          OP_.SetTime(time, time_dof_number, interval,initial);
        }

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

        template<typename DATACONTAINER>
          double
          ElementFunctional(const DATACONTAINER& dc)
          {
            return OP_.ElementFunctional(dc);
          }

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

        double
        PointFunctional(
            const std::map<std::string, const dealii::Vector<double>*> &param_values,
            const std::map<std::string, const VECTOR*> &domain_values)
        {
          return OP_.PointFunctional(param_values, domain_values);
        }

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

        template<typename FACEDATACONTAINER>
          double
          BoundaryFunctional(const FACEDATACONTAINER& fdc)
          {
            return OP_.BoundaryFunctional(fdc);
          }

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

        template<typename FACEDATACONTAINER>
          double
          FaceFunctional(const FACEDATACONTAINER& fdc)
          {
            return OP_.FaceFunctional(fdc);
          }

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

      const dealii::SmartPointer<const dealii::FESystem<dealdim> >
        GetFESystem() const
        {
          return OP_.GetFESystem();
        }

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

        bool
        HasFaces() const
        {
          return OP_.HasFaces();
        }

        /******************************************************/
        bool
        HasPoints() const
        {
          return OP_.HasPoints();
        }

        /******************************************************/
        bool
        HasInterfaces() const
        {
          return OP_.HasInterfaces();
        }

        /******************************************************/
         template<typename ELEMENTITERATOR>
         bool AtInterface(ELEMENTITERATOR& element, unsigned int face) const
         {
           return OP_.AtInterface(element,face);
         }

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

        dealii::UpdateFlags
        GetUpdateFlags() const
        {
          return OP_.GetUpdateFlags();
        }

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

        dealii::UpdateFlags
        GetFaceUpdateFlags() const
        {
          return OP_.GetFaceUpdateFlags();
        }

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

        const std::vector<unsigned int>&
        GetDirichletColors() const
        {
          return OP_.GetDirichletColors();
        }

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

        const std::vector<bool>&
        GetDirichletCompMask(unsigned int color) const
        {
          return OP_.GetDirichletCompMask(color);
        }

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

        const dealii::Function<dealdim>&
        GetDirichletValues(unsigned int color,
            const std::map<std::string, const dealii::Vector<double>*> &param_values,
            const std::map<std::string, const VECTOR*> &domain_values) const
        {
          return OP_.GetDirichletValues(color, param_values, domain_values);
        }

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

        const dealii::Function<dealdim>&
        GetInitialValues() const
        {
          return OP_.GetInitialValues();
        }

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

        const std::vector<unsigned int>&
        GetBoundaryEquationColors() const
        {
          return OP_.GetBoundaryEquationColors();
        }

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

        const std::vector<unsigned int>&
        GetBoundaryFunctionalColors() const
        {
          return OP_.GetBoundaryFunctionalColors();
        }

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

        unsigned int
        GetNFunctionals() const
        {
          return OP_.GetNFunctionals();
        }

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

        unsigned int
        GetNBlocks() const
        {
          return OP_.GetNBlocks();
        }

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

        unsigned int
        GetDoFsPerBlock(unsigned int b) const
        {
          return OP_.GetDoFsPerBlock(b);
        }

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

        const std::vector<unsigned int>&
        GetDoFsPerBlock() const
        {
          return OP_.GetDoFsPerBlock();
        }

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

        const dealii::ConstraintMatrix&
        GetDoFConstraints() const
        {
          return OP_.GetDoFConstraints();
        }

        std::string
        GetType() const
        {
          return OP_.GetType();
        }
        std::string
        GetDoFType() const
        {
          return OP_.GetDoFType();
        }

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

        std::string
        GetFunctionalType() const
        {
          return OP_.GetFunctionalType();
        }

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

        std::string
        GetFunctionalName() const
        {
          return OP_.GetFunctionalName();
        }

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

        DOpEOutputHandler<VECTOR>*
        GetOutputHandler()
        {
          return OP_.GetOutputHandler();
        }

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

        const SpaceTimeHandler<FE, DH, SPARSITYPATTERN, VECTOR, dopedim,
            dealdim>*
        GetSpaceTimeHandler() const
        {
          return OP_.GetBaseProblem().GetSpaceTimeHandler();
        }
        SpaceTimeHandler<FE, DH, SPARSITYPATTERN, VECTOR, dopedim,
            dealdim>*
        GetSpaceTimeHandler()
        {
          return OP_.GetBaseProblem().GetSpaceTimeHandler();
        }

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

        void
        ComputeSparsityPattern(SPARSITYPATTERN & sparsity) const
        {
          OP_.ComputeSparsityPattern(sparsity);
        }
      protected:
        /******************************************************/
        OPTPROBLEM&
        GetProblem()
        {
          return OP_;
        }
        /******************************************************/

        std::string
        GetPart() const
        {
          return part_;
        }

      private:
        OPTPROBLEM& OP_;
        std::string part_;
    };
}
#endif