const.hh

// -*- tab-width: 8; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=8 sw=2 sts=2:
#ifndef DUNE_PDELAB_FUNCTION_CONST_HH
#define DUNE_PDELAB_FUNCTION_CONST_HH
 
#include <dune/pdelab/common/function.hh>
 
namespace Dune {
  namespace PDELab {
 
 
    template<typename GV, typename RF, unsigned dimR = 1>
    class ConstGridFunction
      : public AnalyticGridFunctionBase<
          AnalyticGridFunctionTraits<GV,RF,dimR>,
          ConstGridFunction<GV,RF,dimR>
        >,
        public InstationaryFunctionDefaults
    {
    public:
      typedef AnalyticGridFunctionTraits<GV,RF,dimR> Traits;
      typedef AnalyticGridFunctionBase<
        Traits,
        ConstGridFunction<GV,RF, dimR> > BaseT;
 
 
      ConstGridFunction(const GV& gv,
                        const typename Traits::RangeType& val_ = 1)
        : BaseT(gv)
        , val(val_)
      {}
 
 
      inline void
      evaluateGlobal(const typename Traits::DomainType& x,
                     typename Traits::RangeType& y) const
      {
        y = val;
      }
 
    private:
      typename Traits::RangeType val;
    };
 
 
    template<typename GV, typename RF, unsigned dimR = 1>
    class ConstBoundaryGridFunction
      : public BoundaryGridFunctionBase<
          BoundaryGridFunctionTraits<
            GV,
            RF,dimR,FieldVector<RF,dimR> >,
          ConstBoundaryGridFunction<GV, RF, dimR> >
    {
    public:
      typedef BoundaryGridFunctionTraits<
        GV,
        RF,dimR,FieldVector<RF,dimR> > Traits;
 
    private:
      typedef BoundaryGridFunctionBase<
        Traits,
        ConstBoundaryGridFunction<GV,RF,dimR> > BaseT;
 
    public:
 
      ConstBoundaryGridFunction(const GV& gv_,
                                const typename Traits::RangeType& val_ = 1)
        : gv(gv_)
        , val(val_)
      {}
 
 
      template<typename I>
      inline void
      evaluate(const IntersectionGeometry<I>& ig,
               const typename Traits::DomainType& x,
               typename Traits::RangeType& y) const
      {
        y = val;
      }
 
      inline const GV& getGridView () const
      {
        return gv;
      }
 
    private:
      const GV gv;
      const typename Traits::RangeType val;
    };
 
  } // namspace PDELab
} // namspace Dune
 
#endif // DUNE_PDELAB_FUNCTION_CONST_HH