1/*===-- __clang_openmp_device_functions.h - OpenMP math declares ------ c++ -===
2 *
3 * Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
4 * See https://llvm.org/LICENSE.txt for license information.
5 * SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
6 *
7 *===-----------------------------------------------------------------------===
8 */
9
10#ifndef __CLANG_OPENMP_CMATH_H__
11#define __CLANG_OPENMP_CMATH_H__
12
13#ifndef _OPENMP
14#error "This file is for OpenMP compilation only."
15#endif
16
17#include_next <cmath>
18
19// Make sure we include our math.h overlay, it probably happend already but we
20// need to be sure.
21#include <math.h>
22
23// We (might) need cstdlib because __clang_cuda_cmath.h below declares `abs`
24// which might live in cstdlib.
25#include <cstdlib>
26
27// We need limits because __clang_cuda_cmath.h below uses `std::numeric_limit`.
28#include <limits>
29
30#pragma omp begin declare variant match(                                       \
31    device = {arch(nvptx, nvptx64)}, implementation = {extension(match_any, allow_templates)})
32
33#define __CUDA__
34#define __OPENMP_NVPTX__
35#include <__clang_cuda_cmath.h>
36#undef __OPENMP_NVPTX__
37#undef __CUDA__
38
39// Overloads not provided by the CUDA wrappers but by the CUDA system headers.
40// Since we do not include the latter we define them ourselves.
41#define __DEVICE__ static constexpr __attribute__((always_inline, nothrow))
42
43__DEVICE__ float acosh(float __x) { return ::acoshf(__x); }
44__DEVICE__ float asinh(float __x) { return ::asinhf(__x); }
45__DEVICE__ float atanh(float __x) { return ::atanhf(__x); }
46__DEVICE__ float cbrt(float __x) { return ::cbrtf(__x); }
47__DEVICE__ float erf(float __x) { return ::erff(__x); }
48__DEVICE__ float erfc(float __x) { return ::erfcf(__x); }
49__DEVICE__ float exp2(float __x) { return ::exp2f(__x); }
50__DEVICE__ float expm1(float __x) { return ::expm1f(__x); }
51__DEVICE__ float fdim(float __x, float __y) { return ::fdimf(__x, __y); }
52__DEVICE__ float hypot(float __x, float __y) { return ::hypotf(__x, __y); }
53__DEVICE__ int ilogb(float __x) { return ::ilogbf(__x); }
54__DEVICE__ float lgamma(float __x) { return ::lgammaf(__x); }
55__DEVICE__ long long int llrint(float __x) { return ::llrintf(__x); }
56__DEVICE__ long long int llround(float __x) { return ::llroundf(__x); }
57__DEVICE__ float log1p(float __x) { return ::log1pf(__x); }
58__DEVICE__ float log2(float __x) { return ::log2f(__x); }
59__DEVICE__ float logb(float __x) { return ::logbf(__x); }
60__DEVICE__ long int lrint(float __x) { return ::lrintf(__x); }
61__DEVICE__ long int lround(float __x) { return ::lroundf(__x); }
62__DEVICE__ float nextafter(float __x, float __y) {
63  return ::nextafterf(__x, __y);
64}
65__DEVICE__ float remainder(float __x, float __y) {
66  return ::remainderf(__x, __y);
67}
68__DEVICE__ float scalbln(float __x, long int __y) {
69  return ::scalblnf(__x, __y);
70}
71__DEVICE__ float scalbn(float __x, int __y) { return ::scalbnf(__x, __y); }
72__DEVICE__ float tgamma(float __x) { return ::tgammaf(__x); }
73
74#undef __DEVICE__
75
76#pragma omp end declare variant
77
78#ifdef __AMDGCN__
79#pragma omp begin declare variant match(device = {arch(amdgcn)})
80
81#pragma push_macro("__constant__")
82#define __constant__ __attribute__((constant))
83#define __OPENMP_AMDGCN__
84
85#include <__clang_hip_cmath.h>
86
87#pragma pop_macro("__constant__")
88#undef __OPENMP_AMDGCN__
89
90// Define overloads otherwise which are absent
91#define __DEVICE__ static constexpr __attribute__((always_inline, nothrow))
92
93__DEVICE__ float acos(float __x) { return ::acosf(__x); }
94__DEVICE__ float acosh(float __x) { return ::acoshf(__x); }
95__DEVICE__ float asin(float __x) { return ::asinf(__x); }
96__DEVICE__ float asinh(float __x) { return ::asinhf(__x); }
97__DEVICE__ float atan(float __x) { return ::atanf(__x); }
98__DEVICE__ float atan2(float __x, float __y) { return ::atan2f(__x, __y); }
99__DEVICE__ float atanh(float __x) { return ::atanhf(__x); }
100__DEVICE__ float cbrt(float __x) { return ::cbrtf(__x); }
101__DEVICE__ float cosh(float __x) { return ::coshf(__x); }
102__DEVICE__ float erf(float __x) { return ::erff(__x); }
103__DEVICE__ float erfc(float __x) { return ::erfcf(__x); }
104__DEVICE__ float exp2(float __x) { return ::exp2f(__x); }
105__DEVICE__ float expm1(float __x) { return ::expm1f(__x); }
106__DEVICE__ float fdim(float __x, float __y) { return ::fdimf(__x, __y); }
107__DEVICE__ float hypot(float __x, float __y) { return ::hypotf(__x, __y); }
108__DEVICE__ int ilogb(float __x) { return ::ilogbf(__x); }
109__DEVICE__ float ldexp(float __arg, int __exp) {
110  return ::ldexpf(__arg, __exp);
111}
112__DEVICE__ float lgamma(float __x) { return ::lgammaf(__x); }
113__DEVICE__ float log1p(float __x) { return ::log1pf(__x); }
114__DEVICE__ float logb(float __x) { return ::logbf(__x); }
115__DEVICE__ float nextafter(float __x, float __y) {
116  return ::nextafterf(__x, __y);
117}
118__DEVICE__ float remainder(float __x, float __y) {
119  return ::remainderf(__x, __y);
120}
121__DEVICE__ float scalbn(float __x, int __y) { return ::scalbnf(__x, __y); }
122__DEVICE__ float sinh(float __x) { return ::sinhf(__x); }
123__DEVICE__ float tan(float __x) { return ::tanf(__x); }
124__DEVICE__ float tanh(float __x) { return ::tanhf(__x); }
125__DEVICE__ float tgamma(float __x) { return ::tgammaf(__x); }
126
127#undef __DEVICE__
128
129#pragma omp end declare variant
130#endif // __AMDGCN__
131
132#endif
133