15#ifndef dealii_numbers_h
16#define dealii_numbers_h
31#ifdef DEAL_II_WITH_ADOLC
59 template <
typename Number>
81#if DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 512
83#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 256
85#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128
105#if DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128 && defined(__ALTIVEC__)
107#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 512 && defined(__AVX512F__)
109#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 256 && defined(__AVX__)
111#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128 && defined(__SSE2__)
113#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128 && defined(__ARM_NEON)
125template <
typename Number,
133#ifdef DEAL_II_WITH_ADOLC
148abs(
const adouble &x);
151abs(
const adtl::adouble &x);
159 template <
typename Number, std::
size_t w
idth>
160 DEAL_II_ALWAYS_INLINE ::VectorizedArray<Number, width>
161 sqrt(const ::VectorizedArray<Number, width> &);
162 template <
typename Number, std::
size_t w
idth>
163 DEAL_II_ALWAYS_INLINE ::VectorizedArray<Number, width>
164 abs(const ::VectorizedArray<Number, width> &);
165 template <
typename Number, std::
size_t w
idth>
166 DEAL_II_ALWAYS_INLINE ::VectorizedArray<Number, width>
167 max(const ::VectorizedArray<Number, width> &,
168 const ::VectorizedArray<Number, width> &);
169 template <
typename Number, std::
size_t w
idth>
170 DEAL_II_ALWAYS_INLINE ::VectorizedArray<Number, width>
171 min(const ::VectorizedArray<Number, width> &,
172 const ::VectorizedArray<Number, width> &);
173 template <
typename Number,
size_t w
idth>
175 pow(const ::VectorizedArray<Number, width> &,
const Number p);
176 template <
typename Number,
size_t w
idth>
178 sin(const ::VectorizedArray<Number, width> &);
179 template <
typename Number,
size_t w
idth>
181 cos(const ::VectorizedArray<Number, width> &);
182 template <
typename Number,
size_t w
idth>
184 tan(const ::VectorizedArray<Number, width> &);
185 template <
typename Number,
size_t w
idth>
187 exp(const ::VectorizedArray<Number, width> &);
188 template <
typename Number,
size_t w
idth>
190 log(const ::VectorizedArray<Number, width> &);
215 constexpr double E = 2.7182818284590452354;
220 constexpr double LOG2E = 1.4426950408889634074;
225 constexpr double LOG10E = 0.43429448190325182765;
230 constexpr double LN2 = 0.69314718055994530942;
235 constexpr double LN10 = 2.30258509299404568402;
240 constexpr double PI = 3.14159265358979323846;
245 constexpr double PI_2 = 1.57079632679489661923;
250 constexpr double PI_4 = 0.78539816339744830962;
255 constexpr double SQRT2 = 1.41421356237309504880;
260 constexpr double SQRT1_2 = 0.70710678118654752440;
279 is_finite(
const std::complex<double> &x);
297 is_finite(
const std::complex<long double> &x);
309 template <
typename Number1,
typename Number2>
323 template <
typename Number1,
typename Number2>
334 template <
typename Number>
348 template <
typename Number1,
typename Number2>
362 template <
typename Number1,
typename Number2>
365 const Number2 &value_2);
379 template <
typename Number1,
typename Number2>
393 template <
typename Number1,
typename Number2>
396 const Number2 &value_2);
406 template <
typename number>
462 template <
typename number>
488 static constexpr std::complex<number>
489 conjugate(
const std::complex<number> &x);
505 abs(
const std::complex<number> &x);
513 return std::isnan(x);
521 return std::isfinite(x);
554 template <
typename number>
563 template <
typename number>
572 template <
typename number>
578#ifdef DEAL_II_WITH_ADOLC
590 template <
typename number>
591 constexpr std::complex<number>
599 template <
typename number>
605#ifdef DEAL_II_WITH_ADOLC
614 template <
typename number>
632 template <
typename T>
637 template <
typename NumberType>
649 template <
typename From,
typename To>
654 template <
typename T>
657 template <
typename F,
typename T>
658 static constexpr auto
659 test(
int) ->
decltype(
f(
static_cast<T
>(std::declval<F>())),
true)
659 test(
int) ->
decltype(
f(
static_cast<T
>(std::declval<F>())),
true) {
…}
664 template <
typename F,
typename T>
665 static constexpr auto
672 static const bool value = test<From, To>(0);
679 template <
typename T>
696 template <
typename F>
699 std::enable_if_t<!std::is_same_v<std::decay_t<T>, std::decay_t<F>> &&
700 std::is_constructible_v<T, F>> * =
nullptr)
706 template <
typename F>
709 std::enable_if_t<!std::is_same_v<std::decay_t<T>, std::decay_t<F>> &&
710 !std::is_constructible_v<T, F> &&
714 return static_cast<T
>(f);
721 template <
typename F>
725 std::enable_if_t<!std::is_same_v<std::decay_t<T>, std::decay_t<F>> &&
726 !std::is_constructible_v<T, F> &&
734 template <
typename T>
737 static constexpr const std::complex<T> &
743 static constexpr std::complex<T>
746 return std::complex<T>(t);
750 template <
typename U>
751 static constexpr std::complex<T>
763#ifdef DEAL_II_ADOLC_WITH_ADVANCED_BRANCHING
790 template <
typename Number>
810 template <
typename Number>
845 template <
typename Number>
866 template <
typename Number>
878 template <
typename Number1,
typename Number2>
886 template <
typename Number1,
typename Number2>
887 inline constexpr bool
894 template <
typename Number>
902 template <
typename Number1,
typename Number2>
910 template <
typename Number1,
typename Number2>
919 template <
typename Number1,
typename Number2>
927 template <
typename Number1,
typename Number2>
930 const Number2 &value_2)
#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)