16 #ifndef dealii_numbers_h
17 #define dealii_numbers_h
24 #ifdef DEAL_II_WITH_CUDA
25 # include <cuComplex.h>
28 #include <Kokkos_Macros.hpp>
33 #include <type_traits>
35 #define DEAL_II_HOST_DEVICE KOKKOS_FUNCTION
36 #define DEAL_II_CUDA_HOST_DEV DEAL_II_HOST_DEVICE
37 #define DEAL_II_HOST_DEVICE_ALWAYS_INLINE KOKKOS_FORCEINLINE_FUNCTION
43 #if (defined(__clang__) && defined(__CUDA__)) || defined(KOKKOS_ENABLE_HIP)
44 # define DEAL_II_HOST __host__
51 #ifdef DEAL_II_WITH_ADOLC
79 template <
typename Number>
101 #if DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 512
103 #elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 256
105 #elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128
125 #if DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128 && defined(__ALTIVEC__)
127 #elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 512 && defined(__AVX512F__)
129 #elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 256 && defined(__AVX__)
131 #elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128 && defined(__SSE2__)
133 #elif DEAL_II_VECTORIZATION_WIDTH_IN_BITS >= 128 && defined(__ARM_NEON)
145 template <
typename Number,
149 template <
typename T>
153 #ifdef DEAL_II_WITH_ADOLC
168 abs(
const adouble &x);
171 abs(
const adtl::adouble &x);
179 template <
typename Number, std::
size_t w
idth>
180 DEAL_II_ALWAYS_INLINE ::VectorizedArray<Number, width>
181 sqrt(const ::VectorizedArray<Number, width> &);
182 template <
typename Number, std::
size_t w
idth>
183 DEAL_II_ALWAYS_INLINE ::VectorizedArray<Number, width>
184 abs(const ::VectorizedArray<Number, width> &);
185 template <
typename Number, std::
size_t w
idth>
186 DEAL_II_ALWAYS_INLINE ::VectorizedArray<Number, width>
187 max(const ::VectorizedArray<Number, width> &,
188 const ::VectorizedArray<Number, width> &);
189 template <
typename Number, std::
size_t w
idth>
190 DEAL_II_ALWAYS_INLINE ::VectorizedArray<Number, width>
191 min(const ::VectorizedArray<Number, width> &,
192 const ::VectorizedArray<Number, width> &);
193 template <
typename Number,
size_t w
idth>
195 pow(const ::VectorizedArray<Number, width> &,
const Number p);
196 template <
typename Number,
size_t w
idth>
198 sin(const ::VectorizedArray<Number, width> &);
199 template <
typename Number,
size_t w
idth>
201 cos(const ::VectorizedArray<Number, width> &);
202 template <
typename Number,
size_t w
idth>
204 tan(const ::VectorizedArray<Number, width> &);
205 template <
typename Number,
size_t w
idth>
207 exp(const ::VectorizedArray<Number, width> &);
208 template <
typename Number,
size_t w
idth>
210 log(const ::VectorizedArray<Number, width> &);
235 static constexpr
double E = 2.7182818284590452354;
240 static constexpr
double LOG2E = 1.4426950408889634074;
245 static constexpr
double LOG10E = 0.43429448190325182765;
250 static constexpr
double LN2 = 0.69314718055994530942;
255 static constexpr
double LN10 = 2.30258509299404568402;
260 static constexpr
double PI = 3.14159265358979323846;
265 static constexpr
double PI_2 = 1.57079632679489661923;
270 static constexpr
double PI_4 = 0.78539816339744830962;
275 static constexpr
double SQRT2 = 1.41421356237309504880;
280 static constexpr
double SQRT1_2 = 0.70710678118654752440;
299 is_finite(
const std::complex<double> &x);
317 is_finite(
const std::complex<long double> &x);
329 template <
typename Number1,
typename Number2>
343 template <
typename Number1,
typename Number2>
354 template <
typename Number>
368 template <
typename Number1,
typename Number2>
382 template <
typename Number1,
typename Number2>
385 const Number2 &value_2);
399 template <
typename Number1,
typename Number2>
413 template <
typename Number1,
typename Number2>
416 const Number2 &value_2);
426 template <
typename number>
474 abs(
const number &x);
482 template <
typename number>
508 static constexpr std::complex<number>
509 conjugate(
const std::complex<number> &x);
525 abs(
const std::complex<number> &x);
533 return std::isnan(x);
541 return std::isfinite(x);
574 template <
typename number>
583 template <
typename number>
592 template <
typename number>
598 #ifdef DEAL_II_WITH_ADOLC
610 template <
typename number>
611 constexpr std::complex<number>
619 template <
typename number>
625 #ifdef DEAL_II_WITH_ADOLC
634 template <
typename number>
652 template <
typename T>
657 template <
typename NumberType>
669 template <
typename From,
typename To>
674 template <
typename T>
677 template <
typename F,
typename T>
678 static constexpr
auto
679 test(
int) -> decltype(
f(
static_cast<T>(std::declval<F>())),
true)
684 template <
typename F,
typename T>
685 static constexpr
auto
692 static const bool value = test<From, To>(0);
699 template <
typename T>
716 template <
typename F>
719 std::enable_if_t<!std::is_same_v<std::decay_t<T>, std::decay_t<F>> &&
720 std::is_constructible_v<T, F>> * =
nullptr)
726 template <
typename F>
729 std::enable_if_t<!std::is_same_v<std::decay_t<T>, std::decay_t<F>> &&
730 !std::is_constructible_v<T, F> &&
734 return static_cast<T>(f);
741 template <
typename F>
745 std::enable_if_t<!std::is_same_v<std::decay_t<T>, std::decay_t<F>> &&
746 !std::is_constructible_v<T, F> &&
754 template <
typename T>
757 static constexpr
const std::complex<T> &
763 static constexpr std::complex<T>
766 return std::complex<T>(t);
770 template <
typename U>
771 static constexpr std::complex<T>
779 #ifdef DEAL_II_WITH_CUDA
786 return make_cuComplex(t, 0.f);
793 static cuDoubleComplex
796 return make_cuDoubleComplex(t, 0.);
804 #ifdef DEAL_II_ADOLC_WITH_ADVANCED_BRANCHING
831 template <
typename Number>
851 template <
typename Number>
886 template <
typename Number>
907 template <
typename Number>
919 template <
typename Number1,
typename Number2>
927 template <
typename Number1,
typename Number2>
935 template <
typename Number>
943 template <
typename Number1,
typename Number2>
951 template <
typename Number1,
typename Number2>
960 template <
typename Number1,
typename Number2>
968 template <
typename Number1,
typename Number2>
971 const Number2 &value_2)
VectorizedArray< Number, width > abs(const ::VectorizedArray< Number, width > &x)
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
double norm(const FEValuesBase< dim > &fe, const ArrayView< const std::vector< Tensor< 1, dim >>> &Du)
Tensor< 2, dim, Number > F(const Tensor< 2, dim, Number > &Grad_u)
static constexpr double LOG10E
static constexpr double PI_2
bool value_is_less_than_or_equal_to(const Number1 &value_1, const Number2 &value_2)
static constexpr double E
static constexpr double PI
bool value_is_greater_than(const Number1 &value_1, const Number2 &value_2)
static constexpr double SQRT2
constexpr bool value_is_zero(const Number &value)
bool values_are_not_equal(const Number1 &value_1, const Number2 &value_2)
static constexpr double SQRT1_2
constexpr bool values_are_equal(const Number1 &value_1, const Number2 &value_2)
static constexpr double PI_4
static constexpr double LN10
static constexpr double LN2
bool value_is_less_than(const Number1 &value_1, const Number2 &value_2)
bool is_finite(const double x)
static constexpr double LOG2E
bool value_is_greater_than_or_equal_to(const Number1 &value_1, const Number2 &value_2)
bool is_nan(const double x)
#define DEAL_II_HOST_DEVICE
#define DEAL_II_HOST_DEVICE_ALWAYS_INLINE
static cuComplex value(const float t)
static cuDoubleComplex value(const double t)
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 DEAL_II_HOST_DEVICE_ALWAYS_INLINE 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 DEAL_II_HOST_DEVICE_ALWAYS_INLINE 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 DEAL_II_HOST_DEVICE_ALWAYS_INLINE const T & value(const T &t)
constexpr static 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)