gridfunctionspace.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_PDELAB_GRIDFUNCTIONSPACE_GRIDFUNCTIONSPACE_HH
#define DUNE_PDELAB_GRIDFUNCTIONSPACE_GRIDFUNCTIONSPACE_HH
 
#include <cstddef>
#include <map>
#include <ostream>
#include <set>
#include <vector>
#include <memory>
 
#include <dune/common/deprecated.hh>
#include <dune/common/exceptions.hh>
#include <dune/common/stdstreams.hh>
#include <dune/common/typetraits.hh>
#include <dune/common/shared_ptr.hh>
 
#include <dune/geometry/referenceelements.hh>
#include <dune/geometry/type.hh>
 
#include <dune/localfunctions/common/interfaceswitch.hh>
#include <dune/localfunctions/common/localkey.hh>
 
#include <dune/common/version.hh>
#include <dune/typetree/typetree.hh>
 
// This alias should be removed after a PDELab 2.7 release.
#if DUNE_VERSION_LT_REV(DUNE_TYPETREE,2,7,1)
namespace Dune {
  namespace TypeTree {
    template<std::size_t... i>
    using StaticTreePath = TreePath<i...>;
  }
}
#endif
 
#include <dune/pdelab/common/partitionviewentityset.hh>
#include <dune/pdelab/backend/interface.hh>
 
// we just want the descriptors here, so we temporarily switch off the warning for
// directly including ISTL backend headers
#define _DUNE_PDELAB_SUPPRESS_ISTL_HH_WARNING
#include <dune/pdelab/backend/istl/descriptors.hh>
#undef _DUNE_PDELAB_SUPPRESS_ISTL_HH_WARNING
 
#include <dune/pdelab/constraints/noconstraints.hh>
#include <dune/pdelab/gridfunctionspace/compositegridfunctionspace.hh>
#include <dune/pdelab/gridfunctionspace/datahandleprovider.hh>
#include <dune/pdelab/gridfunctionspace/gridfunctionspacebase.hh>
#include <dune/pdelab/gridfunctionspace/gridfunctionspaceutilities.hh>
#include <dune/pdelab/gridfunctionspace/powergridfunctionspace.hh>
#include <dune/pdelab/gridfunctionspace/tags.hh>
#include <dune/pdelab/ordering/gridviewordering.hh>
#include <dune/pdelab/ordering/lexicographicordering.hh>
 
namespace Dune {
  namespace PDELab {
 
 
#ifndef DOXYGEN
 
    namespace impl {
 
      // Helper structs to avoid compilation failures in the
      // backwards compatibility mode where users stuff a
      // GridView into a GridFunctionSpace.
      // In that case, we cannot extract the GridView type from
      // the GridView itself, so we use a std::conditional in the
      // Traits class to pick either one of the following structs
      // and then use the correct class to do the lookup.
 
      struct _lazy_identity
      {
        template<typename T>
        struct evaluate
        {
          using type = T;
        };
      };
 
      struct _lazy_extract_gridview
      {
        template<typename T>
        struct evaluate
        {
          using type = typename T::GridView;
        };
      };
 
      // Returns a GridView, regardless of whether GV_or_ES is a GridView or an EntitySet
      template<typename GV_or_ES>
      using GridView = typename std::conditional<
        isEntitySet<GV_or_ES>::value,
        impl::_lazy_extract_gridview,
        impl::_lazy_identity
        >::type::template evaluate<GV_or_ES>::type;
 
 
      // Returns an EntitySet, regardless of whether GV_or_ES is a GridView or an EntitySet
      template<typename GV_or_ES>
      using EntitySet = typename std::conditional<
        isEntitySet<GV_or_ES>::value,
        GV_or_ES,
        AllEntitySet<GV_or_ES>
        >::type;
 
    }
 
#endif // DOXYGEN
 
    //=======================================
    // grid function space : single component case
    //=======================================
 
 
    template<typename G, typename L, typename C, typename B, typename O>
    struct GridFunctionSpaceTraits
    {
      static const bool isComposite = false;
 
      using GridView = impl::GridView<G>;
 
      using EntitySet = impl::EntitySet<G>;
 
      using GridViewType = GridView;
 
      typedef B BackendType;
 
      typedef B Backend;
 
      typedef typename B::size_type SizeType;
 
      typedef L FiniteElementMapType;
 
      typedef L FiniteElementMap;
 
      typedef typename L::Traits::FiniteElementType FiniteElementType;
 
      typedef typename L::Traits::FiniteElementType FiniteElement;
 
      typedef C ConstraintsType;
 
 
      typedef O OrderingTag;
 
    };
 
    template<typename GV, typename FEM, typename CE=NoConstraints,
             typename B=ISTL::VectorBackend<>, typename O=DefaultLeafOrderingTag>
    class GridFunctionSpace
      : public TypeTree::LeafNode
      , public GridFunctionSpaceBase<
                 GridFunctionSpace<GV,FEM,CE,B,O>,
                 GridFunctionSpaceTraits<GV,FEM,CE,B,O>
                 >
      , public GridFunctionOutputParameters
      , public DataHandleProvider<GridFunctionSpace<GV,FEM,CE,B,O> >
    {
 
      typedef TypeTree::TransformTree<GridFunctionSpace,gfs_to_ordering<GridFunctionSpace> > ordering_transformation;
 
      template<typename,typename>
      friend class GridFunctionSpaceBase;
 
    public:
      typedef GridFunctionSpaceTraits<GV,FEM,CE,B,O> Traits;
 
    private:
 
      typedef GridFunctionSpaceBase<GridFunctionSpace,Traits> BaseT;
 
    public:
 
      typedef typename GV::Traits::template Codim<0>::Entity Element;
      typedef typename GV::Traits::template Codim<0>::Iterator ElementIterator;
 
      DUNE_DEPRECATED
      typedef O SizeTag;
 
      typedef O OrderingTag;
 
      typedef LeafGridFunctionSpaceTag ImplementationTag;
 
      typedef typename ordering_transformation::Type Ordering;
 
      template<typename E>
      struct ConstraintsContainer
      {
 
        typedef typename std::conditional<
          std::is_same<
            CE,
            NoConstraints
            >::value,
          EmptyTransformation,
          ConstraintsTransformation<typename Ordering::Traits::DOFIndex,typename Ordering::Traits::ContainerIndex,E>
          >::type Type;
 
      private:
        ConstraintsContainer () {}
      };
 
      // ****************************************************************************************************
      // Construct from GridView
      // ****************************************************************************************************
 
      GridFunctionSpace (const typename Traits::GridView& gridview, const FEM& fem, const CE& ce, const B& backend = B(), const OrderingTag& ordering_tag = OrderingTag())
#if DUNE_PDELAB_WARN_ON_GRIDVIEW_BASED_GFS
        DUNE_DEPRECATED_MSG("GridFunctionSpaces now internally use an EntitySet instead of a GridView, please replace the template parameter and the first constructor parameter by an EntitySet").
#endif
        : BaseT(backend,ordering_tag)
        , _es(gridview)
        , pfem(stackobject_to_shared_ptr(fem))
        , _pce(stackobject_to_shared_ptr(ce))
      {
      }
 
      GridFunctionSpace (const typename Traits::GridView& gridview, const std::shared_ptr<const FEM>& fem, const std::shared_ptr<const CE>& ce, const B& backend = B(), const OrderingTag& ordering_tag = OrderingTag())
#if DUNE_PDELAB_WARN_ON_GRIDVIEW_BASED_GFS
        DUNE_DEPRECATED_MSG("GridFunctionSpaces now internally use an EntitySet instead of a GridView, please replace the template parameter and the first constructor parameter by an EntitySet").
#endif
        : BaseT(backend,ordering_tag)
        , _es(gridview)
        , pfem(fem)
        , _pce(ce)
      {}
 
      GridFunctionSpace (const typename Traits::GridView& gridview, const FEM& fem, const B& backend = B(), const OrderingTag& ordering_tag = OrderingTag())
#if DUNE_PDELAB_WARN_ON_GRIDVIEW_BASED_GFS
        DUNE_DEPRECATED_MSG("GridFunctionSpaces now internally use an EntitySet instead of a GridView, please replace the template parameter and the first constructor parameter by an EntitySet").
#endif
        : BaseT(backend,ordering_tag)
        , _es(gridview)
        , pfem(stackobject_to_shared_ptr(fem))
        , _pce(std::make_shared<CE>())
      {}
 
      GridFunctionSpace (const typename Traits::GridView& gridview, const std::shared_ptr<const FEM>& fem, const B& backend = B(), const OrderingTag& ordering_tag = OrderingTag())
#if DUNE_PDELAB_WARN_ON_GRIDVIEW_BASED_GFS
        DUNE_DEPRECATED_MSG("GridFunctionSpaces now internally use an EntitySet instead of a GridView, please replace the template parameter and the first constructor parameter by an EntitySet").
#endif
        : BaseT(backend,ordering_tag)
        , _es(gridview)
        , pfem(fem)
        , _pce(std::make_shared<CE>())
      {}
 
 
      // ****************************************************************************************************
      // Construct from EntitySet
      // ****************************************************************************************************
 
 
      GridFunctionSpace (const typename Traits::EntitySet& entitySet, const FEM& fem, const CE& ce, const B& backend = B(), const OrderingTag& ordering_tag = OrderingTag())
        : BaseT(backend,ordering_tag)
        , _es(entitySet)
        , pfem(stackobject_to_shared_ptr(fem))
        , _pce(stackobject_to_shared_ptr(ce))
      {
      }
 
      GridFunctionSpace (const typename Traits::EntitySet& entitySet, const std::shared_ptr<const FEM>& fem, const std::shared_ptr<const CE>& ce, const B& backend = B(), const OrderingTag& ordering_tag = OrderingTag())
        : BaseT(backend,ordering_tag)
        , _es(entitySet)
        , pfem(fem)
        , _pce(ce)
      {}
 
      GridFunctionSpace (const typename Traits::EntitySet& entitySet, const FEM& fem, const B& backend = B(), const OrderingTag& ordering_tag = OrderingTag())
        : BaseT(backend,ordering_tag)
        , _es(entitySet)
        , pfem(stackobject_to_shared_ptr(fem))
        , _pce(std::make_shared<CE>())
      {}
 
      GridFunctionSpace (const typename Traits::EntitySet& entitySet, const std::shared_ptr<const FEM>& fem, const B& backend = B(), const OrderingTag& ordering_tag = OrderingTag())
        : BaseT(backend,ordering_tag)
        , _es(entitySet)
        , pfem(fem)
        , _pce(std::make_shared<CE>())
      {}
 
 
      const typename Traits::GridView& gridView () const
      {
        return _es.gridView();
      }
 
      const typename Traits::EntitySet& entitySet () const
      {
        return _es;
      }
 
      const FEM& finiteElementMap () const
      {
        return *pfem;
      }
 
      std::shared_ptr<const FEM> finiteElementMapStorage () const
      {
        return pfem;
      }
 
      const typename Traits::ConstraintsType& constraints () const
      {
        return *_pce;
      }
 
      std::shared_ptr<const CE> constraintsStorage () const
      {
        return _pce;
      }
 
      //------------------------------
 
      const Ordering &ordering() const
      {
        if (!this->isRootSpace())
          {
            DUNE_THROW(GridFunctionSpaceHierarchyError,
                       "Ordering can only be obtained for root space in GridFunctionSpace tree.");
          }
        if (!_ordering)
          {
            create_ordering();
            this->update(*_ordering);
          }
        return *_ordering;
      }
 
      Ordering &ordering()
      {
        if (!this->isRootSpace())
          {
            DUNE_THROW(GridFunctionSpaceHierarchyError,
                       "Ordering can only be obtained for root space in GridFunctionSpace tree.");
          }
        if (!_ordering)
          {
            create_ordering();
            this->update(*_ordering);
          }
        return *_ordering;
      }
 
      std::shared_ptr<const Ordering> orderingStorage() const
      {
        if (!this->isRootSpace())
          {
            DUNE_THROW(GridFunctionSpaceHierarchyError,
                       "Ordering can only be obtained for root space in GridFunctionSpace tree.");
          }
        if (!_ordering)
          {
            create_ordering();
            this->update(*_ordering);
          }
        return _ordering;
      }
 
      std::shared_ptr<Ordering> orderingStorage()
      {
        if (!this->isRootSpace())
          {
            DUNE_THROW(GridFunctionSpaceHierarchyError,
                       "Ordering can only be obtained for root space in GridFunctionSpace tree.");
          }
        if (!_ordering)
          {
            create_ordering();
            this->update(*_ordering);
          }
        return _ordering;
      }
 
    private:
 
      // This method here is to avoid a double update of the Ordering when the user calls
      // GFS::update() before GFS::ordering().
      void create_ordering() const
      {
        _ordering = std::make_shared<Ordering>(ordering_transformation::transform(*this));
      }
 
      typename Traits::EntitySet _es;
      std::shared_ptr<FEM const> pfem;
      std::shared_ptr<CE const> _pce;
 
      mutable std::shared_ptr<Ordering> _ordering;
    };
 
 
  } // namespace PDELab
} // namespace Dune
 
#endif // DUNE_PDELAB_GRIDFUNCTIONSPACE_GRIDFUNCTIONSPACE_HH