ts_base.h

Go to the documentation of this file.

#ifndef _TSBase_H_
#define _TSBase_H_

#include "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, const TimeIterator& interval)
        {
          _OP.SetTime(time, interval);
        }

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

        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();
        }

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

        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