functionalinterface.h

Go to the documentation of this file.

#ifndef FUNCTIONAL_INTERFACE_H_
#define FUNCTIONAL_INTERFACE_H_

#include <map>
#include <string>

#include <deal.II/fe/fe_system.h>
#include <deal.II/fe/fe_values.h>
#include <deal.II/fe/mapping.h>

#include "fevalues_wrapper.h"
#include "dofhandler_wrapper.h"
#include "elementdatacontainer.h"
#include "facedatacontainer.h"
#include "multimesh_elementdatacontainer.h"
#include "multimesh_facedatacontainer.h"

namespace DOpE
{
  template<
      template<template<int, int> class DH, typename VECTOR, int dealdim> class EDC,
      template<template<int, int> class DH, typename VECTOR, int dealdim> class FDC,
      template<int, int> class DH, typename VECTOR, int dopedim, int dealdim =
          dopedim>
    class FunctionalInterface
    {
      public:
        FunctionalInterface();
        virtual
        ~FunctionalInterface();

        virtual double
        ElementValue(const EDC<DH, VECTOR, dealdim>& edc);

        virtual void
        ElementValue_U(const EDC<DH, VECTOR, dealdim>& edc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        ElementValue_Q(const EDC<DH, VECTOR, dealdim>& edc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        ElementValue_UU(const EDC<DH, VECTOR, dealdim>& edc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        ElementValue_QU(const EDC<DH, VECTOR, dealdim>& edc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        ElementValue_UQ(const EDC<DH, VECTOR, dealdim>& edc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        ElementValue_QQ(const EDC<DH, VECTOR, dealdim>& edc,
            dealii::Vector<double> &local_vector, double scale);

        virtual double
        PointValue(
            const DOpEWrapper::DoFHandler<dopedim, DH> &control_dof_handler,
            const DOpEWrapper::DoFHandler<dealdim, DH> &state_dof_handler,
            const std::map<std::string, const dealii::Vector<double>*> &param_values,
            const std::map<std::string, const VECTOR*> &domain_values);

        virtual void
        PointValue_U(
            const DOpEWrapper::DoFHandler<dopedim, DH> &control_dof_handler,
            const DOpEWrapper::DoFHandler<dealdim, DH> &state_dof_handler,
            const std::map<std::string, const dealii::Vector<double>*> &param_values,
            const std::map<std::string, const VECTOR*> &domain_values,
            VECTOR& rhs, double scale);

        virtual void
        PointValue_Q(
            const DOpEWrapper::DoFHandler<dopedim, DH> &control_dof_handler,
            const DOpEWrapper::DoFHandler<dealdim, DH> &state_dof_handler,
            const std::map<std::string, const dealii::Vector<double>*> &param_values,
            const std::map<std::string, const VECTOR*> &domain_values,
            VECTOR& rhs, double scale);

        virtual void
        PointValue_UU(
            const DOpEWrapper::DoFHandler<dopedim, DH> &control_dof_handler,
            const DOpEWrapper::DoFHandler<dealdim, DH> &state_dof_handler,
            const std::map<std::string, const dealii::Vector<double>*> &param_values,
            const std::map<std::string, const VECTOR*> &domain_values,
            VECTOR& rhs, double scale);

        virtual void
        PointValue_QU(
            const DOpEWrapper::DoFHandler<dopedim, DH> &control_dof_handler,
            const DOpEWrapper::DoFHandler<dealdim, DH> &state_dof_handler,
            const std::map<std::string, const dealii::Vector<double>*> &param_values,
            const std::map<std::string, const VECTOR*> &domain_values,
            VECTOR& rhs, double scale);

        virtual void
        PointValue_UQ(
            const DOpEWrapper::DoFHandler<dopedim, DH> &control_dof_handler,
            const DOpEWrapper::DoFHandler<dealdim, DH> &state_dof_handler,
            const std::map<std::string, const dealii::Vector<double>*> &param_values,
            const std::map<std::string, const VECTOR*> &domain_values,
            VECTOR& rhs, double scale);

        virtual void
        PointValue_QQ(
            const DOpEWrapper::DoFHandler<dopedim, DH> &control_dof_handler,
            const DOpEWrapper::DoFHandler<dealdim, DH> &state_dof_handler,
            const std::map<std::string, const dealii::Vector<double>*> &param_values,
            const std::map<std::string, const VECTOR*> &domain_values,
            VECTOR& rhs, double scale);

        virtual double
        BoundaryValue(const FDC<DH, VECTOR, dealdim>& fdc);

        virtual void
        BoundaryValue_U(const FDC<DH, VECTOR, dealdim>& fdc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        BoundaryValue_Q(const FDC<DH, VECTOR, dealdim>& fdc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        BoundaryValue_UU(const FDC<DH, VECTOR, dealdim>& fdc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        BoundaryValue_QU(const FDC<DH, VECTOR, dealdim>& fdc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        BoundaryValue_UQ(const FDC<DH, VECTOR, dealdim>& fdc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        BoundaryValue_QQ(const FDC<DH, VECTOR, dealdim>& fdc,
            dealii::Vector<double> &local_vector, double scale);

        virtual double
        FaceValue(const FDC<DH, VECTOR, dealdim>& fdc);

        virtual void
        FaceValue_U(const FDC<DH, VECTOR, dealdim>& fdc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        FaceValue_Q(const FDC<DH, VECTOR, dealdim>& fdc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        FaceValue_UU(const FDC<DH, VECTOR, dealdim>& fdc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        FaceValue_QU(const FDC<DH, VECTOR, dealdim>& fdc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        FaceValue_UQ(const FDC<DH, VECTOR, dealdim>& fdc,
            dealii::Vector<double> &local_vector, double scale);

        virtual void
        FaceValue_QQ(const FDC<DH, VECTOR, dealdim>& fdc,
            dealii::Vector<double> &local_vector, double scale);

        virtual double
        AlgebraicValue(
            const std::map<std::string, const dealii::Vector<double>*> &param_values,
            const std::map<std::string, const VECTOR*> &domain_values);
        virtual void
        AlgebraicGradient_Q(VECTOR& gradient,
            const std::map<std::string, const dealii::Vector<double>*> &param_values,
            const std::map<std::string, const VECTOR*> &domain_values);

        virtual std::string
        GetType() const;
        virtual std::string
        GetName() const;

        virtual bool
        NeedTime() const
        {
          return true;
        }

      virtual void
      SetTime(double /*t*/) const
        {
        }

        virtual dealii::UpdateFlags
        GetUpdateFlags() const;

        virtual dealii::UpdateFlags
        GetFaceUpdateFlags() const;

        virtual bool
        HasFaces() const;

        virtual bool
        HasPoints() const;

        void
        SetProblemType(std::string type);
      protected:
        std::string
        GetProblemType() const;
      private:
        std::string problem_type_;

    };
}

#endif