fieldpromotion.hh

Go to the documentation of this file.

#ifndef __DUNE_ACFEM_EXPRESSIONS_FIELDPROMOTION_HH__
#define __DUNE_ACFEM_EXPRESSIONS_FIELDPROMOTION_HH__
 
#include <limits>
#include <complex>
 
#include "../common/types.hh"
 
#include "policy.hh"
#include "typetraits.hh"
 
namespace Dune
{
 
  namespace ACFem
  {
 
    namespace {
 
      template<class F, class SFINAE = void>
      struct FloatingPointClosureHelper
      {
        using Type = Expressions::Policy::FloatingPointDefault;
      };
 
      template<class F>
      struct FloatingPointClosureHelper<F, std::enable_if_t<!IsDecay<F>::value> >
      {
        using Type = copy_cv_reference_t<F, typename FloatingPointClosureHelper<std::decay_t<F> >::Type>;
      };
 
      template<class F>
      struct FloatingPointClosureHelper<
        F,
        std::enable_if_t<IsDecay<F>::value
                         &&
                         std::numeric_limits<typename FieldTraits<F>::real_type>::is_iec559>
        >
      {
        using Type = typename FieldTraits<F>::real_type;
      };
 
      template<class T>
      struct FloatingPointClosureHelper<std::complex<T> >
      {
        using Type = std::complex<typename FloatingPointClosureHelper<T>::Type>;
      };
 
    } // NS anonymous
 
    template<class T>
    using FloatingPointClosure = typename FloatingPointClosureHelper<T>::Type;
 
    template<class T, std::enable_if_t<!std::is_same<FloatingPointClosure<std::decay_t<T> >, std::decay_t<T> >::value, int> = 0>
    constexpr auto floatingPointClosure(T&& t)
    {
      return FloatingPointClosure<std::decay_t<T> >(t);
    }
 
    template<class T, std::enable_if_t<std::is_same<FloatingPointClosure<std::decay_t<T> >, std::decay_t<T> >::value, int> = 0>
    constexpr decltype(auto) floatingPointClosure(T&& t)
    {
      return std::forward<T>(t);
    }
 
    template<class T1, class T2, class SFINAE = void>
    struct FieldPromotion
    {
      using Type = std::decay_t<decltype(std::declval<typename FieldTraits<std::decay_t<T1> >::field_type>()
                                         *
                                         std::declval<typename FieldTraits<std::decay_t<T2> >::field_type>())>;
    };
 
    template<class T1, class T2>
    using FieldPromotionType = typename FieldPromotion<T1, T2>::Type;
 
 
  } // ACFem
 
} // Dune
 
#endif // __DUNE_ACFEM_EXPRESSIONS_FIELDPROMOTION_HH__