userdefineddofconstraints.h

Go to the documentation of this file.

#ifndef CONSTRAINTMAKER_H_
#define CONSTRAINTMAKER_H_

#include <deal.II/dofs/dof_tools.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/constraint_matrix.h>

#include "dofhandler_wrapper.h"
#include "mapping_wrapper.h"
#include "parameterreader.h"
#include "dopeexception.h"

namespace DOpE
{
  template<template<int, int> class DH, int dopedim, int dealdim = dopedim>
    class UserDefinedDoFConstraints
    {
      public:
        UserDefinedDoFConstraints()
        {
        }
        virtual
        ~UserDefinedDoFConstraints()
        {
        }
        virtual void
        MakeStateDoFConstraints(
            const DOpEWrapper::DoFHandler<dealdim, DH> & dof_handler,
            dealii::ConstraintMatrix& dof_constraints) const;

        virtual void
        MakeControlDoFConstraints(
            const DOpEWrapper::DoFHandler<dopedim, DH> & dof_handler,
            dealii::ConstraintMatrix& dof_constraints) const;

        void
        RegisterMapping(const typename DOpEWrapper::Mapping<dealdim, DH>& mapping)
        {
          mapping_ = &mapping;
        }

    protected:
      const DOpEWrapper::Mapping<dealdim, DH> &
        GetMapping() const
        {
          return *mapping_;
        }
    private:
      const DOpEWrapper::Mapping<dealdim, DH> * mapping_;
    };

  template<template<int, int> class DH, int dopedim, int dealdim>
    void
    UserDefinedDoFConstraints<DH, dopedim, dealdim>::MakeStateDoFConstraints(
        const DOpEWrapper::DoFHandler<dealdim, DH> & /*dof_handler*/,
        dealii::ConstraintMatrix& /*dof_constraints*/) const
    {
      throw DOpEException("Not Implemented.",
          "UserDefinedDoFConstraints::MakeStateDoFConstraints");
    }

  template<template<int, int> class DH, int dopedim, int dealdim>
    void
    UserDefinedDoFConstraints<DH, dopedim, dealdim>::MakeControlDoFConstraints(
        const DOpEWrapper::DoFHandler<dopedim, DH> & /*dof_handler*/,
        dealii::ConstraintMatrix& /*dof_constraints*/) const
    {
      throw DOpEException("Not Implemented.",
          "UserDefinedDoFConstraints::MakeControlDoFConstraints");
    }

} //end of namespace
#endif /* CONSTRAINTMAKER_H_ */