initialproblem.h

Go to the documentation of this file.

#ifndef INITIAL_PROBLEM_H_
#define INITIAL_PROBLEM_H_

#include "spacetimehandler.h"

using namespace dealii;

namespace DOpE
{
  template<typename PDE, typename VECTOR, int dim>
    class InitialProblem
    {
      public:
        InitialProblem(PDE& pde) :
            pde_(pde)
        {
        }

        std::string
        GetName() const
        {
          return "InitialProblem";
        }
        std::string
        GetType() const
        {
          return "initial_state";
        }

        /******************************************************/
        template<typename EDC>
          inline void
          ElementEquation(const EDC& edc,
              dealii::Vector<double> &local_vector, double scale,
              double scale_ico);

        template<typename EDC>
          inline void
          ElementRhs(const EDC& edc,
              dealii::Vector<double> &local_vector, double scale = 1.);

        inline void
        PointRhs(
            const std::map<std::string, const dealii::Vector<double>*> &param_values,
            const std::map<std::string, const VECTOR*> &domain_values,
            VECTOR& rhs_vector, double scale = 1);

        template<typename EDC>
          inline void
          ElementMatrix(const EDC& edc,
              dealii::FullMatrix<double> &local_matrix, double scale = 1.,
              double scale_ico = 1.);

        template<typename FDC>
          inline void
          FaceEquation(const FDC& fdc,
                       dealii::Vector<double> &local_vector, double scale = 1., double scale_ico = 1.);

        template<typename FDC>
          inline void
          InterfaceEquation(const FDC& fdc,
              dealii::Vector<double> &local_vector, double scale = 1., double scale_ico = 1.);

        template<typename FDC>
          inline void
          FaceRhs(const FDC& fdc,
              dealii::Vector<double> &local_vector, double scale = 1.);

       template<typename FDC>
          inline void
          FaceMatrix(const FDC& fdc,
              dealii::FullMatrix<double> &local_matrix, double scale = 1., double scale_ico = 1.);

        template<typename FDC>
          inline void
          InterfaceMatrix(const FDC& fdc,
              dealii::FullMatrix<double> &local_matrix, double scale = 1., double scale_ico = 1.);

        template<typename FDC>
          inline void
          BoundaryEquation(const FDC& fdc,
              dealii::Vector<double> &local_vector, double scale = 1., double scale_ico = 1.);

        template<typename FDC>
          inline void
          BoundaryRhs(const FDC& fdc,
              dealii::Vector<double> &local_vector, double scale = 1.);

        template<typename FDC>
          inline void
          BoundaryMatrix(const FDC& fdc,
              dealii::FullMatrix<double> &local_matrix, double scale = 1., double scale_ico = 1.);

        inline const dealii::SmartPointer<const dealii::FESystem<dim> >
        GetFESystem() const;

        inline const dealii::SmartPointer<const dealii::hp::FECollection<dim> >
        GetFECollection() const;

        inline std::string
        GetDoFType() const;

        inline bool
        HasFaces() const;
        inline bool
        HasPoints() const;
        inline bool
        HasInterfaces() const;

        inline dealii::UpdateFlags
        GetUpdateFlags() const;

        inline dealii::UpdateFlags
        GetFaceUpdateFlags() const;

        /******************************************************/
        inline void
        SetTime(double time, const TimeIterator& interval);

        template<typename SPARSITYPATTERN>
          inline void
          ComputeSparsityPattern(SPARSITYPATTERN & sparsity) const;

        inline const std::vector<unsigned int>&
        GetDirichletColors() const;
        inline const std::vector<bool>&
        GetDirichletCompMask(unsigned int color) const;
        inline const Function<dim>&
        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;
        inline const std::vector<unsigned int>&
        GetBoundaryEquationColors() const;
        inline const dealii::ConstraintMatrix&
        GetDoFConstraints() const;
       const dealii::Function<dim>&
        GetInitialValues() const;
        /******************************************************/
        DOpEOutputHandler<VECTOR>*
        GetOutputHandler()
        {
          return pde_.GetOutputHandler();
        }
        PDE&
        GetBaseProblem()
        {
          return pde_;
        } 
        template<typename ELEMENTITERATOR>
          bool AtInterface(ELEMENTITERATOR& element, unsigned int face) const;

      protected:

      private:
        PDE& pde_;
    };

  /*****************************************************************************************/
  /********************************IMPLEMENTATION*******************************************/
  /*****************************************************************************************/

  template<typename PDE, typename VECTOR, int dim>
    template<typename EDC>
      void
      InitialProblem<PDE, VECTOR, dim>::ElementEquation(const EDC& edc,
          dealii::Vector<double> &local_vector, double scale,
          double scale_ico)
      {
        pde_.Init_ElementEquation(edc, local_vector, scale, scale_ico);
      }

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

  template<typename PDE, typename VECTOR, int dim>
    template<typename FDC>
      void
      InitialProblem<PDE, VECTOR, dim>::FaceEquation(
          const FDC& fdc,
          dealii::Vector<double> &local_vector, double scale, double scale_ico)
      {
        pde_.Init_FaceEquation(fdc, local_vector, scale, scale_ico);
      }

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

  template<typename PDE, typename VECTOR, int dim>
    template<typename FDC>
      void
      InitialProblem<PDE, VECTOR, dim>::InterfaceEquation(
          const FDC& fdc,
          dealii::Vector<double> &local_vector, double scale, double scale_ico)
      {
        pde_.Init_InterfaceEquation(fdc, local_vector, scale, scale_ico);
      }
  /******************************************************/

  template<typename PDE, typename VECTOR, int dim>
    template<typename FDC>
      void
      InitialProblem<PDE, VECTOR, dim>::BoundaryEquation(
          const FDC& fdc,
          dealii::Vector<double> &local_vector, double scale, double scale_ico)
      {
        pde_.Init_BoundaryEquation(fdc, local_vector, scale, scale_ico);
      }

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

  template<typename PDE, typename VECTOR, int dim>
    template<typename EDC>
      void
      InitialProblem<PDE, VECTOR, dim>::ElementRhs(const EDC& edc,
          dealii::Vector<double> &local_vector, double scale)
      {
        pde_.Init_ElementRhs(edc, local_vector, scale);
      }

  /******************************************************/
  template<typename PDE, typename VECTOR, int dim>
    void
    InitialProblem<PDE, VECTOR, dim>::PointRhs(
        const std::map<std::string, const dealii::Vector<double>*> &param_values,
        const std::map<std::string, const VECTOR*> &domain_values,
        VECTOR& rhs_vector, double scale)
    {
      pde_.Init_PointRhs(param_values, domain_values, rhs_vector, scale);
    }

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

  template<typename PDE, typename VECTOR, int dim>
    template<typename FDC>
      void
      InitialProblem<PDE, VECTOR, dim>::FaceRhs(
          const FDC& /*fdc*/,
          dealii::Vector<double> &/*local_vector*/, double /*scale*/)
      {
      }

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

  template<typename PDE, typename VECTOR, int dim>
    template<typename FDC>
      void
      InitialProblem<PDE, VECTOR, dim>::BoundaryRhs(
          const FDC& /*fdc*/,
          dealii::Vector<double> &/*local_vector*/, double /*scale*/)
      {
      }

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

  template<typename PDE, typename VECTOR, int dim>
    template<typename EDC>
      void
      InitialProblem<PDE, VECTOR, dim>::ElementMatrix(const EDC& edc,
          dealii::FullMatrix<double> &local_matrix, double scale, double scale_ico)
      {
        pde_.Init_ElementMatrix(edc, local_matrix, scale, scale_ico);
      }

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

  template<typename PDE, typename VECTOR, int dim>
    template<typename FDC>
      void
      InitialProblem<PDE, VECTOR, dim>::FaceMatrix(const FDC& fdc,
          FullMatrix<double> &local_matrix, double scale, double scale_ico)
      {
        pde_.Init_FaceMatrix(fdc, local_matrix, scale, scale_ico);
      }

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

  template<typename PDE, typename VECTOR, int dim>
    template<typename FDC>
      void
      InitialProblem<PDE, VECTOR, dim>::InterfaceMatrix(
          const FDC& fdc, FullMatrix<double> &local_matrix, double scale, double scale_ico)
      {
        pde_.Init_InterfaceMatrix(fdc, local_matrix, scale, scale_ico);
      }

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

  template<typename PDE, typename VECTOR, int dim>
    template<typename FDC>
      void
      InitialProblem<PDE, VECTOR, dim>::BoundaryMatrix(
          const FDC& fdc, FullMatrix<double> &local_matrix, double scale, double scale_ico)
      {
        pde_.Init_BoundaryMatrix(fdc, local_matrix, scale, scale_ico);
      }

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

  template<typename PDE, typename VECTOR, int dim>
    std::string
    InitialProblem<PDE, VECTOR, dim>::GetDoFType() const
    {
      return pde_.GetDoFType();
    }

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

  template<typename PDE, typename VECTOR, int dim>
    const SmartPointer<const dealii::FESystem<dim> >
    InitialProblem<PDE, VECTOR, dim>::GetFESystem() const
    {
      return pde_.GetFESystem();
    }

  /******************************************************/
  template<typename PDE, typename VECTOR, int dim>
    const SmartPointer<const dealii::hp::FECollection<dim> >
    InitialProblem<PDE, VECTOR, dim>::GetFECollection() const
    {
      return pde_.GetFECollection();
    }

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

  template<typename PDE, typename VECTOR, int dim>
    UpdateFlags
    InitialProblem<PDE, VECTOR, dim>::GetUpdateFlags() const
    {
      UpdateFlags r;
      r = pde_.GetUpdateFlags();
      return r | update_JxW_values;
    }

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

  template<typename PDE, typename VECTOR, int dim>
    UpdateFlags
    InitialProblem<PDE, VECTOR, dim>::GetFaceUpdateFlags() const
    {
      UpdateFlags r;
      r = pde_.GetFaceUpdateFlags();
      return r | update_JxW_values;
    }

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

  template<typename PDE, typename VECTOR, int dim>
    void
    InitialProblem<PDE, VECTOR, dim>::SetTime(double time,
        const TimeIterator& interval)
    {
      pde_.SetTime(time, interval);
    }

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

  template<typename PDE, typename VECTOR, int dim>
    template<typename SPARSITYPATTERN>
      void
      InitialProblem<PDE, VECTOR, dim>::ComputeSparsityPattern(
          SPARSITYPATTERN & sparsity) const
      {
        pde_.ComputeStateSparsityPattern(sparsity);
      }

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

  template<typename PDE, typename VECTOR, int dim>
    bool
    InitialProblem<PDE, VECTOR, dim>::HasFaces() const
    {
      return pde_.HasFaces();
    }

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

  template<typename PDE, typename VECTOR, int dim>
    bool
    InitialProblem<PDE, VECTOR, dim>::HasPoints() const
    {
      return pde_.HasPoints();
    }

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

  template<typename PDE, typename VECTOR, int dim>
    bool
    InitialProblem<PDE, VECTOR, dim>::HasInterfaces() const
    {
      return pde_.HasInterfaces();
    }

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

  template<typename PDE, typename VECTOR, int dim>
    const std::vector<unsigned int>&
    InitialProblem<PDE, VECTOR, dim>::GetDirichletColors() const
    {
      return pde_.GetDirichletColors();
    }

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

  template<typename PDE, typename VECTOR, int dim>
    const std::vector<bool>&
    InitialProblem<PDE, VECTOR, dim>::GetDirichletCompMask(
        unsigned int color) const
    {
      return pde_.GetDirichletCompMask(color);
    }

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

  template<typename PDE, typename VECTOR, int dim>
    const Function<dim>&
    InitialProblem<PDE, VECTOR, dim>::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 pde_.GetDirichletValues(color, param_values, domain_values);
    }

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

  template<typename PDE, typename VECTOR, int dim>
    const std::vector<unsigned int>&
    InitialProblem<PDE, VECTOR, dim>::GetBoundaryEquationColors() const
    {
      return pde_.GetBoundaryEquationColors();
    }

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

  template<typename PDE, typename VECTOR, int dim>
    const dealii::ConstraintMatrix&
    InitialProblem<PDE, VECTOR, dim>::GetDoFConstraints() const
    {
      return pde_.GetDoFConstraints();
    }

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

  template<typename PDE, typename VECTOR, int dim>
    const dealii::Function<dim>&
    InitialProblem<PDE, VECTOR, dim>::GetInitialValues() const
    {
      return pde_.GetInitialValues();
    }

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

  template<typename PDE, typename VECTOR, int dim>
    template<typename ELEMENTITERATOR>
    bool InitialProblem<PDE, VECTOR, dim>
                      ::AtInterface(ELEMENTITERATOR& element, unsigned int face) const
  {
    return pde_.AtInterface(element,face);
  }

}
#endif