15#ifndef dealii_numbers_h
16#define dealii_numbers_h
23#include <Kokkos_Macros.hpp>
32#ifdef DEAL_II_WITH_ADOLC
60 template <
typename Number>
82#if DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 512
84#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 256
86#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128
106#if DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128 && defined(__ALTIVEC__)
108#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 512 && defined(__AVX512F__)
110#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 256 && defined(__AVX__)
112#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128 && defined(__SSE2__)
114#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128 && defined(__ARM_NEON)
126template <
typename Number,
134#ifdef DEAL_II_WITH_ADOLC
152abs(
const adtl::adouble &
x);
160 template <
typename Number, std::
size_t w
idth>
163 template <
typename Number, std::
size_t w
idth>
166 template <
typename Number, std::
size_t w
idth>
170 template <
typename Number, std::
size_t w
idth>
174 template <
typename Number,
size_t w
idth>
177 template <
typename Number,
size_t w
idth>
180 template <
typename Number,
size_t w
idth>
183 template <
typename Number,
size_t w
idth>
186 template <
typename Number,
size_t w
idth>
189 template <
typename Number,
size_t w
idth>
216 constexpr double E = 2.7182818284590452354;
221 constexpr double LOG2E = 1.4426950408889634074;
226 constexpr double LOG10E = 0.43429448190325182765;
231 constexpr double LN2 = 0.69314718055994530942;
236 constexpr double LN10 = 2.30258509299404568402;
241 constexpr double PI = 3.14159265358979323846;
246 constexpr double PI_2 = 1.57079632679489661923;
251 constexpr double PI_4 = 0.78539816339744830962;
256 constexpr double SQRT2 = 1.41421356237309504880;
261 constexpr double SQRT1_2 = 0.70710678118654752440;
298 is_finite(
const std::complex<long double> &
x);
310 template <
typename Number1,
typename Number2>
324 template <
typename Number1,
typename Number2>
335 template <
typename Number>
349 template <
typename Number1,
typename Number2>
363 template <
typename Number1,
typename Number2>
380 template <
typename Number1,
typename Number2>
394 template <
typename Number1,
typename Number2>
407 template <
typename number>
463 template <
typename number>
489 static constexpr std::complex<number>
506 abs(
const std::complex<number> &
x);
514 return std::isnan(
x);
522 return std::isfinite(
x);
555 template <
typename number>
564 template <
typename number>
573 template <
typename number>
579#ifdef DEAL_II_WITH_ADOLC
591 template <
typename number>
592 constexpr std::complex<number>
600 template <
typename number>
606#ifdef DEAL_II_WITH_ADOLC
615 template <
typename number>
633 template <
typename T>
638 template <
typename NumberType>
650 template <
typename From,
typename To>
655 template <
typename T>
658 template <
typename F,
typename T>
659 static constexpr auto
660 test(
int) ->
decltype(
f(
static_cast<T
>(std::declval<F>())),
true)
665 template <
typename F,
typename T>
666 static constexpr auto
680 template <
typename T>
697 template <
typename F>
700 std::enable_if_t<!std::is_same_v<std::decay_t<T>, std::decay_t<F>> &&
701 std::is_constructible_v<T, F>> * =
nullptr)
707 template <
typename F>
710 std::enable_if_t<!std::is_same_v<std::decay_t<T>, std::decay_t<F>> &&
711 !std::is_constructible_v<T, F> &&
715 return static_cast<T
>(f);
722 template <
typename F>
726 std::enable_if_t<!std::is_same_v<std::decay_t<T>, std::decay_t<F>> &&
727 !std::is_constructible_v<T, F> &&
735 template <
typename T>
738 static constexpr const std::complex<T> &
744 static constexpr std::complex<T>
747 return std::complex<T>(t);
751 template <
typename U>
752 static constexpr std::complex<T>
764#ifdef DEAL_II_ADOLC_WITH_ADVANCED_BRANCHING
791 template <
typename Number>
811 template <
typename Number>
846 template <
typename Number>
867 template <
typename Number>
879 template <
typename Number1,
typename Number2>
887 template <
typename Number1,
typename Number2>
888 inline constexpr bool
895 template <
typename Number>
903 template <
typename Number1,
typename Number2>
911 template <
typename Number1,
typename Number2>
920 template <
typename Number1,
typename Number2>
928 template <
typename Number1,
typename Number2>
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_HOST_DEVICE
#define DEAL_II_NAMESPACE_CLOSE
#define DEAL_II_HOST_DEVICE_ALWAYS_INLINE
bool value_is_less_than_or_equal_to(const Number1 &value_1, const Number2 &value_2)
bool value_is_greater_than(const Number1 &value_1, const Number2 &value_2)
constexpr bool values_are_not_equal(const Number1 &value_1, const Number2 &value_2)
constexpr bool value_is_zero(const Number &value)
constexpr bool values_are_equal(const Number1 &value_1, const Number2 &value_2)
bool value_is_less_than(const Number1 &value_1, const Number2 &value_2)
bool is_finite(const double x)
bool value_is_greater_than_or_equal_to(const Number1 &value_1, const Number2 &value_2)
bool is_nan(const double x)
::VectorizedArray< Number, width > log(const ::VectorizedArray< Number, width > &)
::VectorizedArray< Number, width > exp(const ::VectorizedArray< Number, width > &)
::VectorizedArray< Number, width > tan(const ::VectorizedArray< Number, width > &)
::VectorizedArray< Number, width > min(const ::VectorizedArray< Number, width > &, const ::VectorizedArray< Number, width > &)
::VectorizedArray< Number, width > max(const ::VectorizedArray< Number, width > &, const ::VectorizedArray< Number, width > &)
::VectorizedArray< Number, width > cos(const ::VectorizedArray< Number, width > &)
::VectorizedArray< Number, width > sin(const ::VectorizedArray< Number, width > &)
::VectorizedArray< Number, width > sqrt(const ::VectorizedArray< Number, width > &)
::VectorizedArray< Number, width > pow(const ::VectorizedArray< Number, width > &, const Number p)
::VectorizedArray< Number, width > abs(const ::VectorizedArray< Number, width > &)
static constexpr std::complex< T > value(const std::complex< U > &t)
static constexpr std::complex< T > value(const T &t)
static constexpr const std::complex< T > & value(const std::complex< T > &t)
static constexpr T value(const F &f, std::enable_if_t<!std::is_same_v< std::decay_t< T >, std::decay_t< F > > &&std::is_constructible_v< T, F > > *=nullptr)
static constexpr const T & value(const T &t)
static constexpr T value(const F &f, std::enable_if_t<!std::is_same_v< std::decay_t< T >, std::decay_t< F > > &&!std::is_constructible_v< T, F > &&is_explicitly_convertible< const F, T >::value > *=nullptr)
static T value(const F &f, std::enable_if_t<!std::is_same_v< std::decay_t< T >, std::decay_t< F > > &&!std::is_constructible_v< T, F > &&!is_explicitly_convertible< const F, T >::value &&Differentiation::AD::is_ad_number< F >::value > *=nullptr)
static constexpr unsigned int max_width
static constexpr auto test(...) -> bool
static constexpr auto test(int) -> decltype(f(static_cast< T >(std::declval< F >())), true)
std::complex< double > double_type
static constexpr const number & conjugate(const number &x)
static constexpr bool is_complex
static real_type abs(const number &x)
static constexpr real_type abs_square(const number &x)