domains.hh

Go to the documentation of this file.

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
 
// SPDX-FileCopyrightText: Copyright © DUNE-FUFEM Project contributors, see file AUTHORS.md
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception OR LGPL-3.0-or-later
 
#ifndef DUNE_FUFEM_FORMS_DOMAINS_HH
#define DUNE_FUFEM_FORMS_DOMAINS_HH
 
#include <type_traits>
#include <utility>
 
#include <dune/common/version.hh>
 
#if DUNE_VERSION_GTE(DUNE_FUNCTIONS, 2, 11)
#include <dune/functions/common/subdomain.hh>
#endif
 
#include <dune/fufem/boundarypatch.hh>
 
 
 
namespace Dune::Fufem::Forms {
 
 
 
  template<class T>
  struct IsBulkDomain : public std::false_type {};
 
 
 
  template<class T>
  struct IsBoundaryDomain : public std::false_type {};
 
 
 
  template<class T>
  struct IsSkeletonDomain : public std::false_type {};
 
 
 
  template<class T>
  struct IsDomain : public std::disjunction<IsBulkDomain<T>, IsBoundaryDomain<T>, IsSkeletonDomain<T>> {};
 
 
 
  class Bulk
  {
  public:
    template<class Element>
    static constexpr bool contains(const Element&)
    {
      return true;
    }
  };
 
 
 
  template<>
  struct IsBulkDomain<Bulk> : public std::true_type {};
 
 
 
  class Boundary
  {
  public:
    template<class Intersection>
    static constexpr bool contains(const Intersection& intersection)
    {
      return true;
    }
 
    template<class Element>
    static constexpr bool containsFaceOf(const Element& element)
    {
      return element.hasBoundaryIntersections();
    }
  };
 
 
 
  template<>
  struct IsBoundaryDomain<Boundary> : public std::true_type {};
 
 
 
  class Skeleton
  {
  public:
    template<class Intersection>
    static constexpr bool contains(const Intersection&)
    {
      return true;
    }
  };
 
 
 
  template<>
  struct IsSkeletonDomain<Skeleton> : public std::true_type {};
 
 
 
  template<class GridView>
  struct IsBoundaryDomain<Dune::Fufem::Boundary<GridView>> : public std::true_type {};
 
 
 
  template<class GridView>
  struct IsBoundaryDomain<Dune::Fufem::BoundaryPatch<GridView>> : public std::true_type {};
 
 
 
  template<class GridView>
  struct IsBoundaryDomain<::BoundaryPatch<GridView>> : public std::true_type {};
 
 
 
#if DUNE_VERSION_GTE(DUNE_FUNCTIONS, 2, 11)
  template<class GridView>
  struct IsBulkDomain<Dune::Functions::Experimental::SubDomain<GridView>> : public std::true_type {};
 
 
 
  template<class SubDomainA, class SubDomainB>
  struct IsSkeletonDomain<Dune::Functions::Experimental::SubDomainInterface<SubDomainA, SubDomainB>> : public std::true_type {};
 
 
 
  template<class SubDomain>
  struct IsSkeletonDomain<Dune::Functions::Experimental::SubDomainSkeleton<SubDomain>> : public std::true_type {};
#endif
 
 
 
  inline constexpr Bulk bulk = {};
 
 
 
  inline constexpr Boundary boundary = {};
 
 
 
  inline constexpr Skeleton skeleton = {};
 
 
 
} // namespace Dune::Fufem::Forms
 
 
#endif // DUNE_FUFEM_FORMS_DOMAINS_HH