quadraturehints.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_QUADRATUREHINTS_HH
#define DUNE_FUFEM_FORMS_QUADRATUREHINTS_HH
 
#include <type_traits>
#include <utility>
 
#include <dune/common/typetraits.hh>
#include <dune/common/std/functional.hh>
 
#include <dune/fufem/quadraturerules/quadraturerulecache.hh>
 
 
 
namespace Dune::Fufem::Forms::QuadratureHints {
 
 
 
  template<class T>
  struct IsQuadratureHint : public std::false_type {};
 
 
 
  class LiftOrder
  {
  public:
 
    LiftOrder(int increment)
      : increment_(increment)
    {}
 
    QuadratureRuleKey operator()(QuadratureRuleKey key) const
    {
      key.setOrder(key.order() + increment_);
      return key;
    }
 
  private:
    int increment_;
  };
 
  template<>
  struct IsQuadratureHint<LiftOrder> : public std::true_type {};
 
 
 
  class FixedOrder
  {
  public:
 
    FixedOrder(int order)
      : order_(order)
    {}
 
    QuadratureRuleKey operator()(QuadratureRuleKey key) const
    {
      key.setOrder(order_);
      return key;
    }
 
  private:
    int order_;
  };
 
  template<>
  struct IsQuadratureHint<FixedOrder> : public std::true_type {};
 
 
 
  template<class T>
  class QuadratureHint
  {
  public:
 
    QuadratureHint(T transformation)
      : transformation_(std::move(transformation))
    {}
 
    QuadratureHint()
    requires std::is_default_constructible_v<T>
      : transformation_()
    {}
 
    decltype(auto) operator()(QuadratureRuleKey&& key) const
    {
      return transformation_(std::move(key));
    }
 
    decltype(auto) operator()(const QuadratureRuleKey& key) const
    {
      return transformation_(key);
    }
 
  private:
    T transformation_;
  };
 
  template<class T>
  struct IsQuadratureHint<QuadratureHint<T>> : public std::true_type {};
 
 
 
  using KeepOrder = QuadratureHint<Dune::Std::identity>;
 
 
 
} // namespace Dune::Fufem::Forms::QuadratureHints
 
 
 
#endif // DUNE_FUFEM_FORMS_QUADRATUREHINTS_HH