baseclass.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_BASECLASS_HH
#define DUNE_FUFEM_FORMS_BASECLASS_HH
 
#include <array>
#include <cstddef>
#include <type_traits>
#include <tuple>
#include <utility>
 
#include <dune/common/typetraits.hh>
#include <dune/common/indices.hh>
#include <dune/common/rangeutilities.hh>
 
 
 
namespace Dune::Fufem::Forms::Impl {
 
  template<class Tuple, class F>
  void forEachTupleEntry(Tuple&& tuple, F&& f){
    std::apply([&](auto&&... entry){
      (f(std::forward<decltype(entry)>(entry)), ...);
    }, tuple);
  }
 
 
 
} // namespace Dune::Fufem::Forms::Impl
 
 
 
namespace Dune::Fufem::Forms {
 
 
 
  template<std::size_t k>
  class MultilinearOperator {
  public:
    static constexpr std::size_t arity = k;
  };
 
  template<std::size_t k>
  class UnaryOperator : public MultilinearOperator<1>
  {
  public:
    static constexpr std::size_t argIndex = k;
  };
 
  // Forward declaration
  template<class Operator, class... Operators>
  class SumOperator;
 
 
 
namespace Impl {
 
 
  std::false_type isOperatorHelper(const void*)
  {
    return {};
  }
 
  template<std::size_t k>
  std::true_type isOperatorHelper(const MultilinearOperator<k>*)
  {
    return {};
  }
 
  template<class T>
  struct IsSumOperatorHelper : std::false_type {};
 
  template<class... Op>
  struct IsSumOperatorHelper<SumOperator<Op...>> : std::true_type {};
 
}
 
 
 
  template<class Op>
  using IsOperator = decltype(Impl::isOperatorHelper(std::declval<const std::decay_t<Op>*>()));
 
  template<class Op>
  inline constexpr bool isOperator_v = IsOperator<Op>::value;
 
 
 
  template<class Op>
  using IsSumOperator = Impl::IsSumOperatorHelper<std::decay_t<Op>>;
 
  template<class Op>
  inline constexpr bool isSumOperator_v = IsSumOperator<Op>::value;
 
 
 
  template<class Op>
  using IsOperatorOrSumOperator = std::disjunction<IsOperator<Op>, IsSumOperator<Op>>;
 
  template<class Op>
  inline constexpr bool isOperatorOrSumOperator_v = IsOperatorOrSumOperator<Op>::value;
 
} // namespace Dune::Fufem::Forms
 
 
 
 
namespace Dune::Fufem::Forms::Impl {
 
  // Execute callback on first LocalView in list, that belongs to the basis
  // and return true if any matching LocalView was found.
  template<class F, class Basis, class... LocalViews>
  bool visitMatchingLocalView(F&& f, const Basis& basis, LocalViews&&... localViews)
  {
    auto visitIfMatching = [&](const auto& localView) {
      if constexpr (std::is_same_v<Basis, typename std::decay_t<decltype(localView)>::GlobalBasis>)
        if (&localView.globalBasis() == &basis)
        {
          f(localView);
          return true;
        }
      return false;
    };
    return (false || ... || visitIfMatching(localViews));
  }
 
} // namespace Dune::Fufem::Forms::Impl
 
 
 
#endif // DUNE_FUFEM_FORMS_BASECLASS_HH