constraintsparameters.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
 
#ifndef DUNE_PDELAB_CONSTRAINTS_COMMON_CONSTRAINTSPARAMETERS_HH
#define DUNE_PDELAB_CONSTRAINTS_COMMON_CONSTRAINTSPARAMETERS_HH
 
#include <dune/common/fvector.hh>
#include <dune/typetree/typetree.hh>
 
namespace Dune {
  namespace PDELab {
 
 
    struct DirichletConstraintsParameters :
      public TypeTree::LeafNode
    {
 
      template<typename I>
      bool isDirichlet(const I & intersection, const FieldVector<typename I::ctype, I::mydimension> & coord) const
      {
        return true;
      }
 
      template<typename I>
      bool isNeumann(const I & intersection, const FieldVector<typename I::ctype, I::mydimension> & coord) const
      {
        return false;
      }
 
 
      template<typename T>
      void setTime(const T& time)
      {}
 
    };
 
    struct NoDirichletConstraintsParameters :
      public TypeTree::LeafNode
    {
 
      template<typename I>
      bool isDirichlet(const I & intersection, const FieldVector<typename I::ctype, I::mydimension> & coord) const
      {
        return false;
      }
 
      template<typename I>
      bool isNeumann(const I & intersection, const FieldVector<typename I::ctype, I::mydimension> & coord) const
      {
        return true;
      }
 
      template<typename T>
      void setTime(const T& time)
      {}
 
    };
 
 
    struct FluxConstraintsParameters :
      public TypeTree::LeafNode
    {
 
      template<typename I>
      bool isNeumann(const I & intersection, const FieldVector<typename I::ctype, I::mydimension> & coord) const
      {
        return true;
      }
 
      template<typename T>
      void setTime(const T& time)
      {}
 
    };
 
    struct NoFluxConstraintsParameters :
      public TypeTree::LeafNode
    {
 
      template<typename I>
      bool isNeumann(const I & intersection, const FieldVector<typename I::ctype, I::mydimension> & coord) const
      {
        return false;
      }
 
      template<typename T>
      void setTime(const T& time)
      {}
 
    };
 
    template<typename DirichletConstraintsParameters>
    struct FluxFromDirichletConstraintsAdapter :
      public TypeTree::LeafNode
    {
 
      template<typename I>
      bool isDirichlet(const I & intersection, const FieldVector<typename I::ctype, I::mydimension> & coord) const
      {
        return _dirichletConstraintsParameters.isDirichlet(intersection,coord);
      }
 
      template<typename I>
      bool isNeumann(const I & intersection, const FieldVector<typename I::ctype, I::mydimension> & coord) const
      {
        return !_dirichletConstraintsParameters.isDirichlet(intersection,coord);
      }
 
      template<typename T>
      void setTime(const T& time)
      {
        _dirichletConstraintsParameters.setTime(time);
      }
 
      FluxFromDirichletConstraintsAdapter(DirichletConstraintsParameters& dirichletConstraintsParameters)
        : _dirichletConstraintsParameters(dirichletConstraintsParameters)
      {}
 
    private:
 
      DirichletConstraintsParameters& _dirichletConstraintsParameters;
 
    };
 
    template<typename FluxConstraintsParameters>
    struct DirichletFromFluxConstraintsAdapter :
      public TypeTree::LeafNode
    {
 
      template<typename I>
      bool isDirichlet(const I & intersection, const FieldVector<typename I::ctype, I::mydimension> & coord) const
      {
        return !_fluxConstraintsParameters.isNeumann(intersection,coord);
      }
 
      template<typename I>
      bool isNeumann(const I & intersection, const FieldVector<typename I::ctype, I::mydimension> & coord) const
      {
        return _fluxConstraintsParameters.isNeumann(intersection,coord);
      }
 
      template<typename T>
      void setTime(const T& time)
      {
        _fluxConstraintsParameters.setTime(time);
      }
 
      DirichletFromFluxConstraintsAdapter(FluxConstraintsParameters& fluxConstraintsParameters)
        : _fluxConstraintsParameters(fluxConstraintsParameters)
      {}
 
    private:
 
      FluxConstraintsParameters& _fluxConstraintsParameters;
 
    };
 
    class NoConstraintsParameters : public TypeTree::LeafNode {};
 
 
  }
}
 
#endif // DUNE_PDELAB_CONSTRAINTS_COMMON_CONSTRAINTSPARAMETERS_HH