15#ifndef dealii_numbers_h
16#define dealii_numbers_h
30#ifdef DEAL_II_WITH_ADOLC
58 template <
typename Number>
80#if DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 512
82#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 256
84#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128
104#if DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128 && defined(__ALTIVEC__)
106#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 512 && defined(__AVX512F__)
108#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 256 && defined(__AVX__)
110#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128 && defined(__SSE2__)
112#elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128 && defined(__ARM_NEON)
124template <
typename Number,
132#ifdef DEAL_II_WITH_ADOLC
147abs(
const adouble &x);
150abs(
const adtl::adouble &x);
158 template <
typename Number, std::
size_t w
idth>
159 DEAL_II_ALWAYS_INLINE ::VectorizedArray<Number, width>
160 sqrt(const ::VectorizedArray<Number, width> &);
161 template <
typename Number, std::
size_t w
idth>
162 DEAL_II_ALWAYS_INLINE ::VectorizedArray<Number, width>
163 abs(const ::VectorizedArray<Number, width> &);
164 template <
typename Number, std::
size_t w
idth>
165 DEAL_II_ALWAYS_INLINE ::VectorizedArray<Number, width>
166 max(const ::VectorizedArray<Number, width> &,
167 const ::VectorizedArray<Number, width> &);
168 template <
typename Number, std::
size_t w
idth>
169 DEAL_II_ALWAYS_INLINE ::VectorizedArray<Number, width>
170 min(const ::VectorizedArray<Number, width> &,
171 const ::VectorizedArray<Number, width> &);
172 template <
typename Number,
size_t w
idth>
174 pow(const ::VectorizedArray<Number, width> &,
const Number p);
175 template <
typename Number,
size_t w
idth>
177 sin(const ::VectorizedArray<Number, width> &);
178 template <
typename Number,
size_t w
idth>
180 cos(const ::VectorizedArray<Number, width> &);
181 template <
typename Number,
size_t w
idth>
183 tan(const ::VectorizedArray<Number, width> &);
184 template <
typename Number,
size_t w
idth>
186 exp(const ::VectorizedArray<Number, width> &);
187 template <
typename Number,
size_t w
idth>
189 log(const ::VectorizedArray<Number, width> &);
214 constexpr double E = 2.7182818284590452354;
219 constexpr double LOG2E = 1.4426950408889634074;
224 constexpr double LOG10E = 0.43429448190325182765;
229 constexpr double LN2 = 0.69314718055994530942;
234 constexpr double LN10 = 2.30258509299404568402;
239 constexpr double PI = 3.14159265358979323846;
244 constexpr double PI_2 = 1.57079632679489661923;
249 constexpr double PI_4 = 0.78539816339744830962;
254 constexpr double SQRT2 = 1.41421356237309504880;
259 constexpr double SQRT1_2 = 0.70710678118654752440;
278 is_finite(
const std::complex<double> &x);
296 is_finite(
const std::complex<long double> &x);
308 template <
typename Number1,
typename Number2>
322 template <
typename Number1,
typename Number2>
333 template <
typename Number>
347 template <
typename Number1,
typename Number2>
361 template <
typename Number1,
typename Number2>
364 const Number2 &value_2);
378 template <
typename Number1,
typename Number2>
392 template <
typename Number1,
typename Number2>
395 const Number2 &value_2);
405 template <
typename number>
461 template <
typename number>
487 static constexpr std::complex<number>
488 conjugate(
const std::complex<number> &x);
504 abs(
const std::complex<number> &x);
512 return std::isnan(x);
520 return std::isfinite(x);
553 template <
typename number>
562 template <
typename number>
571 template <
typename number>
577#ifdef DEAL_II_WITH_ADOLC
589 template <
typename number>
590 constexpr std::complex<number>
598 template <
typename number>
604#ifdef DEAL_II_WITH_ADOLC
613 template <
typename number>
631 template <
typename T>
636 template <
typename NumberType>
648 template <
typename From,
typename To>
653 template <
typename T>
656 template <
typename F,
typename T>
657 static constexpr auto
658 test(
int) ->
decltype(
f(
static_cast<T
>(std::declval<F>())),
true)
658 test(
int) ->
decltype(
f(
static_cast<T
>(std::declval<F>())),
true) {
…}
663 template <
typename F,
typename T>
664 static constexpr auto
671 static const bool value = test<From, To>(0);
678 template <
typename T>
695 template <
typename F>
698 std::enable_if_t<!std::is_same_v<std::decay_t<T>, std::decay_t<F>> &&
699 std::is_constructible_v<T, F>> * =
nullptr)
705 template <
typename F>
708 std::enable_if_t<!std::is_same_v<std::decay_t<T>, std::decay_t<F>> &&
709 !std::is_constructible_v<T, F> &&
713 return static_cast<T
>(f);
720 template <
typename F>
724 std::enable_if_t<!std::is_same_v<std::decay_t<T>, std::decay_t<F>> &&
725 !std::is_constructible_v<T, F> &&
733 template <
typename T>
736 static constexpr const std::complex<T> &
742 static constexpr std::complex<T>
745 return std::complex<T>(t);
749 template <
typename U>
750 static constexpr std::complex<T>
762#ifdef DEAL_II_ADOLC_WITH_ADVANCED_BRANCHING
789 template <
typename Number>
809 template <
typename Number>
844 template <
typename Number>
865 template <
typename Number>
877 template <
typename Number1,
typename Number2>
885 template <
typename Number1,
typename Number2>
886 inline constexpr bool
893 template <
typename Number>
901 template <
typename Number1,
typename Number2>
909 template <
typename Number1,
typename Number2>
918 template <
typename Number1,
typename Number2>
926 template <
typename Number1,
typename Number2>
929 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)