decay.hh

#ifndef __DUNE_ACFEM_EXPRESSIONS_DECAY_HH__
#define __DUNE_ACFEM_EXPRESSIONS_DECAY_HH__
 
#include "transform.hh"
#include "subexpression.hh"
 
namespace Dune::ACFem::Expressions {
 
  template<class Optimize = OptimizeTop>
  struct OperateDecay
  {
    template<class F, class... Arg>
    constexpr auto operator()(F&& f, Arg&&... arg) const
    {
      return subExpressionExpression(operate(Optimize{}, std::decay_t<F>(f), subExpressionExpression(std::decay_t<Arg>(arg))...));
    }
  };
 
  template<class E, class Optimize = OptimizeTop, class TreePos = IndexSequence<> >
  auto decayCopy(E&& e, Optimize = Optimize{}, TreePos = TreePos{})
  {
    return transform(std::forward<E>(e), FalseType{}, OperateDecay<Optimize>{}, TreePos{});
  }
 
  template<class T, class TreePos = IndexSequence<> >
  using DecayExpression = std::decay_t<decltype(decayCopy(std::declval<T>(), OptimizeTop{}, TreePos{}))>;
 
}
 
#endif // __DUNE_ACFEM_EXPRESSIONS_DECAY_HH__