1// -*- C++ -*-
2//===----------------------------------------------------------------------===//
3//
4// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
5// See https://llvm.org/LICENSE.txt for license information.
6// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
7//
8//===----------------------------------------------------------------------===//
9
10#ifndef _LIBCPP_CMATH
11#define _LIBCPP_CMATH
12
13/*
14 cmath synopsis
15
16Macros:
17
18 HUGE_VAL
19 HUGE_VALF // C99
20 HUGE_VALL // C99
21 INFINITY // C99
22 NAN // C99
23 FP_INFINITE // C99
24 FP_NAN // C99
25 FP_NORMAL // C99
26 FP_SUBNORMAL // C99
27 FP_ZERO // C99
28 FP_FAST_FMA // C99
29 FP_FAST_FMAF // C99
30 FP_FAST_FMAL // C99
31 FP_ILOGB0 // C99
32 FP_ILOGBNAN // C99
33 MATH_ERRNO // C99
34 MATH_ERREXCEPT // C99
35 math_errhandling // C99
36
37namespace std
38{
39
40Types:
41
42 float_t // C99
43 double_t // C99
44
45// C90
46
47floating_point abs(floating_point x);
48
49floating_point acos (arithmetic x);
50float acosf(float x);
51long double acosl(long double x);
52
53floating_point asin (arithmetic x);
54float asinf(float x);
55long double asinl(long double x);
56
57floating_point atan (arithmetic x);
58float atanf(float x);
59long double atanl(long double x);
60
61floating_point atan2 (arithmetic y, arithmetic x);
62float atan2f(float y, float x);
63long double atan2l(long double y, long double x);
64
65floating_point ceil (arithmetic x);
66float ceilf(float x);
67long double ceill(long double x);
68
69floating_point cos (arithmetic x);
70float cosf(float x);
71long double cosl(long double x);
72
73floating_point cosh (arithmetic x);
74float coshf(float x);
75long double coshl(long double x);
76
77floating_point exp (arithmetic x);
78float expf(float x);
79long double expl(long double x);
80
81floating_point fabs (arithmetic x);
82float fabsf(float x);
83long double fabsl(long double x);
84
85floating_point floor (arithmetic x);
86float floorf(float x);
87long double floorl(long double x);
88
89floating_point fmod (arithmetic x, arithmetic y);
90float fmodf(float x, float y);
91long double fmodl(long double x, long double y);
92
93floating_point frexp (arithmetic value, int* exp);
94float frexpf(float value, int* exp);
95long double frexpl(long double value, int* exp);
96
97floating_point ldexp (arithmetic value, int exp);
98float ldexpf(float value, int exp);
99long double ldexpl(long double value, int exp);
100
101floating_point log (arithmetic x);
102float logf(float x);
103long double logl(long double x);
104
105floating_point log10 (arithmetic x);
106float log10f(float x);
107long double log10l(long double x);
108
109floating_point modf (floating_point value, floating_point* iptr);
110float modff(float value, float* iptr);
111long double modfl(long double value, long double* iptr);
112
113floating_point pow (arithmetic x, arithmetic y);
114float powf(float x, float y);
115long double powl(long double x, long double y);
116
117floating_point sin (arithmetic x);
118float sinf(float x);
119long double sinl(long double x);
120
121floating_point sinh (arithmetic x);
122float sinhf(float x);
123long double sinhl(long double x);
124
125floating_point sqrt (arithmetic x);
126float sqrtf(float x);
127long double sqrtl(long double x);
128
129floating_point tan (arithmetic x);
130float tanf(float x);
131long double tanl(long double x);
132
133floating_point tanh (arithmetic x);
134float tanhf(float x);
135long double tanhl(long double x);
136
137// C99
138
139bool signbit(arithmetic x); // constexpr since C++23
140
141int fpclassify(arithmetic x); // constexpr since C++23
142
143bool isfinite(arithmetic x); // constexpr since C++23
144bool isinf(arithmetic x); // constexpr since C++23
145bool isnan(arithmetic x); // constexpr since C++23
146bool isnormal(arithmetic x); // constexpr since C++23
147
148bool isgreater(arithmetic x, arithmetic y); // constexpr since C++23
149bool isgreaterequal(arithmetic x, arithmetic y); // constexpr since C++23
150bool isless(arithmetic x, arithmetic y); // constexpr since C++23
151bool islessequal(arithmetic x, arithmetic y); // constexpr since C++23
152bool islessgreater(arithmetic x, arithmetic y); // constexpr since C++23
153bool isunordered(arithmetic x, arithmetic y); // constexpr since C++23
154
155floating_point acosh (arithmetic x);
156float acoshf(float x);
157long double acoshl(long double x);
158
159floating_point asinh (arithmetic x);
160float asinhf(float x);
161long double asinhl(long double x);
162
163floating_point atanh (arithmetic x);
164float atanhf(float x);
165long double atanhl(long double x);
166
167floating_point cbrt (arithmetic x);
168float cbrtf(float x);
169long double cbrtl(long double x);
170
171floating_point copysign (arithmetic x, arithmetic y);
172float copysignf(float x, float y);
173long double copysignl(long double x, long double y);
174
175floating_point erf (arithmetic x);
176float erff(float x);
177long double erfl(long double x);
178
179floating_point erfc (arithmetic x);
180float erfcf(float x);
181long double erfcl(long double x);
182
183floating_point exp2 (arithmetic x);
184float exp2f(float x);
185long double exp2l(long double x);
186
187floating_point expm1 (arithmetic x);
188float expm1f(float x);
189long double expm1l(long double x);
190
191floating_point fdim (arithmetic x, arithmetic y);
192float fdimf(float x, float y);
193long double fdiml(long double x, long double y);
194
195floating_point fma (arithmetic x, arithmetic y, arithmetic z);
196float fmaf(float x, float y, float z);
197long double fmal(long double x, long double y, long double z);
198
199floating_point fmax (arithmetic x, arithmetic y);
200float fmaxf(float x, float y);
201long double fmaxl(long double x, long double y);
202
203floating_point fmin (arithmetic x, arithmetic y);
204float fminf(float x, float y);
205long double fminl(long double x, long double y);
206
207double hermite(unsigned n, double x); // C++17
208float hermite(unsigned n, float x); // C++17
209long double hermite(unsigned n, long double x); // C++17
210float hermitef(unsigned n, float x); // C++17
211long double hermitel(unsigned n, long double x); // C++17
212template <class Integer>
213double hermite(unsigned n, Integer x); // C++17
214
215floating_point hypot (arithmetic x, arithmetic y);
216float hypotf(float x, float y);
217long double hypotl(long double x, long double y);
218
219double hypot(double x, double y, double z); // C++17
220float hypot(float x, float y, float z); // C++17
221long double hypot(long double x, long double y, long double z); // C++17
222
223int ilogb (arithmetic x);
224int ilogbf(float x);
225int ilogbl(long double x);
226
227floating_point lgamma (arithmetic x);
228float lgammaf(float x);
229long double lgammal(long double x);
230
231long long llrint (arithmetic x);
232long long llrintf(float x);
233long long llrintl(long double x);
234
235long long llround (arithmetic x);
236long long llroundf(float x);
237long long llroundl(long double x);
238
239floating_point log1p (arithmetic x);
240float log1pf(float x);
241long double log1pl(long double x);
242
243floating_point log2 (arithmetic x);
244float log2f(float x);
245long double log2l(long double x);
246
247floating_point logb (arithmetic x);
248float logbf(float x);
249long double logbl(long double x);
250
251long lrint (arithmetic x);
252long lrintf(float x);
253long lrintl(long double x);
254
255long lround (arithmetic x);
256long lroundf(float x);
257long lroundl(long double x);
258
259double nan (const char* str);
260float nanf(const char* str);
261long double nanl(const char* str);
262
263floating_point nearbyint (arithmetic x);
264float nearbyintf(float x);
265long double nearbyintl(long double x);
266
267floating_point nextafter (arithmetic x, arithmetic y);
268float nextafterf(float x, float y);
269long double nextafterl(long double x, long double y);
270
271floating_point nexttoward (arithmetic x, long double y);
272float nexttowardf(float x, long double y);
273long double nexttowardl(long double x, long double y);
274
275floating_point remainder (arithmetic x, arithmetic y);
276float remainderf(float x, float y);
277long double remainderl(long double x, long double y);
278
279floating_point remquo (arithmetic x, arithmetic y, int* pquo);
280float remquof(float x, float y, int* pquo);
281long double remquol(long double x, long double y, int* pquo);
282
283floating_point rint (arithmetic x);
284float rintf(float x);
285long double rintl(long double x);
286
287floating_point round (arithmetic x);
288float roundf(float x);
289long double roundl(long double x);
290
291floating_point scalbln (arithmetic x, long ex);
292float scalblnf(float x, long ex);
293long double scalblnl(long double x, long ex);
294
295floating_point scalbn (arithmetic x, int ex);
296float scalbnf(float x, int ex);
297long double scalbnl(long double x, int ex);
298
299floating_point tgamma (arithmetic x);
300float tgammaf(float x);
301long double tgammal(long double x);
302
303floating_point trunc (arithmetic x);
304float truncf(float x);
305long double truncl(long double x);
306
307constexpr float lerp(float a, float b, float t) noexcept; // C++20
308constexpr double lerp(double a, double b, double t) noexcept; // C++20
309constexpr long double lerp(long double a, long double b, long double t) noexcept; // C++20
310
311} // std
312
313*/
314
315#if __cplusplus < 201103L && defined(_LIBCPP_USE_FROZEN_CXX03_HEADERS)
316# include <__cxx03/cmath>
317#else
318# include <__config>
319# include <__type_traits/enable_if.h>
320# include <__type_traits/is_floating_point.h>
321# include <version>
322
323# include <__cmath/lerp.h>
324# include <__cmath/special_functions.h>
325# include <math.h>
326
327# ifndef _LIBCPP_MATH_H
328# error <cmath> tried including <math.h> but didn't find libc++'s <math.h> header. \
329 This usually means that your header search paths are not configured properly. \
330 The header search paths should contain the C++ Standard Library headers before \
331 any C Standard Library, and you are probably using compiler flags that make that \
332 not be the case.
333# endif
334
335# if !defined(_LIBCPP_HAS_NO_PRAGMA_SYSTEM_HEADER)
336# pragma GCC system_header
337# endif
338
339_LIBCPP_PUSH_MACROS
340# include <__undef_macros>
341
342_LIBCPP_BEGIN_NAMESPACE_STD
343
344using ::signbit _LIBCPP_USING_IF_EXISTS;
345using ::fpclassify _LIBCPP_USING_IF_EXISTS;
346using ::isfinite _LIBCPP_USING_IF_EXISTS;
347using ::isinf _LIBCPP_USING_IF_EXISTS;
348using ::isnan _LIBCPP_USING_IF_EXISTS;
349using ::isnormal _LIBCPP_USING_IF_EXISTS;
350using ::isgreater _LIBCPP_USING_IF_EXISTS;
351using ::isgreaterequal _LIBCPP_USING_IF_EXISTS;
352using ::isless _LIBCPP_USING_IF_EXISTS;
353using ::islessequal _LIBCPP_USING_IF_EXISTS;
354using ::islessgreater _LIBCPP_USING_IF_EXISTS;
355using ::isunordered _LIBCPP_USING_IF_EXISTS;
356using ::isunordered _LIBCPP_USING_IF_EXISTS;
357
358using ::float_t _LIBCPP_USING_IF_EXISTS;
359using ::double_t _LIBCPP_USING_IF_EXISTS;
360
361using ::abs _LIBCPP_USING_IF_EXISTS;
362
363using ::acos _LIBCPP_USING_IF_EXISTS;
364using ::acosf _LIBCPP_USING_IF_EXISTS;
365using ::asin _LIBCPP_USING_IF_EXISTS;
366using ::asinf _LIBCPP_USING_IF_EXISTS;
367using ::atan _LIBCPP_USING_IF_EXISTS;
368using ::atanf _LIBCPP_USING_IF_EXISTS;
369using ::atan2 _LIBCPP_USING_IF_EXISTS;
370using ::atan2f _LIBCPP_USING_IF_EXISTS;
371using ::ceil _LIBCPP_USING_IF_EXISTS;
372using ::ceilf _LIBCPP_USING_IF_EXISTS;
373using ::cos _LIBCPP_USING_IF_EXISTS;
374using ::cosf _LIBCPP_USING_IF_EXISTS;
375using ::cosh _LIBCPP_USING_IF_EXISTS;
376using ::coshf _LIBCPP_USING_IF_EXISTS;
377
378using ::exp _LIBCPP_USING_IF_EXISTS;
379using ::expf _LIBCPP_USING_IF_EXISTS;
380
381using ::fabs _LIBCPP_USING_IF_EXISTS;
382using ::fabsf _LIBCPP_USING_IF_EXISTS;
383using ::floor _LIBCPP_USING_IF_EXISTS;
384using ::floorf _LIBCPP_USING_IF_EXISTS;
385
386using ::fmod _LIBCPP_USING_IF_EXISTS;
387using ::fmodf _LIBCPP_USING_IF_EXISTS;
388
389using ::frexp _LIBCPP_USING_IF_EXISTS;
390using ::frexpf _LIBCPP_USING_IF_EXISTS;
391using ::ldexp _LIBCPP_USING_IF_EXISTS;
392using ::ldexpf _LIBCPP_USING_IF_EXISTS;
393
394using ::log _LIBCPP_USING_IF_EXISTS;
395using ::logf _LIBCPP_USING_IF_EXISTS;
396
397using ::log10 _LIBCPP_USING_IF_EXISTS;
398using ::log10f _LIBCPP_USING_IF_EXISTS;
399using ::modf _LIBCPP_USING_IF_EXISTS;
400using ::modff _LIBCPP_USING_IF_EXISTS;
401
402using ::pow _LIBCPP_USING_IF_EXISTS;
403using ::powf _LIBCPP_USING_IF_EXISTS;
404
405using ::sin _LIBCPP_USING_IF_EXISTS;
406using ::sinf _LIBCPP_USING_IF_EXISTS;
407using ::sinh _LIBCPP_USING_IF_EXISTS;
408using ::sinhf _LIBCPP_USING_IF_EXISTS;
409
410using ::sqrt _LIBCPP_USING_IF_EXISTS;
411using ::sqrtf _LIBCPP_USING_IF_EXISTS;
412using ::tan _LIBCPP_USING_IF_EXISTS;
413using ::tanf _LIBCPP_USING_IF_EXISTS;
414
415using ::tanh _LIBCPP_USING_IF_EXISTS;
416using ::tanhf _LIBCPP_USING_IF_EXISTS;
417
418using ::acosh _LIBCPP_USING_IF_EXISTS;
419using ::acoshf _LIBCPP_USING_IF_EXISTS;
420using ::asinh _LIBCPP_USING_IF_EXISTS;
421using ::asinhf _LIBCPP_USING_IF_EXISTS;
422using ::atanh _LIBCPP_USING_IF_EXISTS;
423using ::atanhf _LIBCPP_USING_IF_EXISTS;
424using ::cbrt _LIBCPP_USING_IF_EXISTS;
425using ::cbrtf _LIBCPP_USING_IF_EXISTS;
426
427using ::copysign _LIBCPP_USING_IF_EXISTS;
428using ::copysignf _LIBCPP_USING_IF_EXISTS;
429
430using ::erf _LIBCPP_USING_IF_EXISTS;
431using ::erff _LIBCPP_USING_IF_EXISTS;
432using ::erfc _LIBCPP_USING_IF_EXISTS;
433using ::erfcf _LIBCPP_USING_IF_EXISTS;
434using ::exp2 _LIBCPP_USING_IF_EXISTS;
435using ::exp2f _LIBCPP_USING_IF_EXISTS;
436using ::expm1 _LIBCPP_USING_IF_EXISTS;
437using ::expm1f _LIBCPP_USING_IF_EXISTS;
438using ::fdim _LIBCPP_USING_IF_EXISTS;
439using ::fdimf _LIBCPP_USING_IF_EXISTS;
440using ::fmaf _LIBCPP_USING_IF_EXISTS;
441using ::fma _LIBCPP_USING_IF_EXISTS;
442using ::fmax _LIBCPP_USING_IF_EXISTS;
443using ::fmaxf _LIBCPP_USING_IF_EXISTS;
444using ::fmin _LIBCPP_USING_IF_EXISTS;
445using ::fminf _LIBCPP_USING_IF_EXISTS;
446using ::hypot _LIBCPP_USING_IF_EXISTS;
447using ::hypotf _LIBCPP_USING_IF_EXISTS;
448using ::ilogb _LIBCPP_USING_IF_EXISTS;
449using ::ilogbf _LIBCPP_USING_IF_EXISTS;
450using ::lgamma _LIBCPP_USING_IF_EXISTS;
451using ::lgammaf _LIBCPP_USING_IF_EXISTS;
452using ::llrint _LIBCPP_USING_IF_EXISTS;
453using ::llrintf _LIBCPP_USING_IF_EXISTS;
454using ::llround _LIBCPP_USING_IF_EXISTS;
455using ::llroundf _LIBCPP_USING_IF_EXISTS;
456using ::log1p _LIBCPP_USING_IF_EXISTS;
457using ::log1pf _LIBCPP_USING_IF_EXISTS;
458using ::log2 _LIBCPP_USING_IF_EXISTS;
459using ::log2f _LIBCPP_USING_IF_EXISTS;
460using ::logb _LIBCPP_USING_IF_EXISTS;
461using ::logbf _LIBCPP_USING_IF_EXISTS;
462using ::lrint _LIBCPP_USING_IF_EXISTS;
463using ::lrintf _LIBCPP_USING_IF_EXISTS;
464using ::lround _LIBCPP_USING_IF_EXISTS;
465using ::lroundf _LIBCPP_USING_IF_EXISTS;
466
467using ::nan _LIBCPP_USING_IF_EXISTS;
468using ::nanf _LIBCPP_USING_IF_EXISTS;
469
470using ::nearbyint _LIBCPP_USING_IF_EXISTS;
471using ::nearbyintf _LIBCPP_USING_IF_EXISTS;
472using ::nextafter _LIBCPP_USING_IF_EXISTS;
473using ::nextafterf _LIBCPP_USING_IF_EXISTS;
474using ::nexttoward _LIBCPP_USING_IF_EXISTS;
475using ::nexttowardf _LIBCPP_USING_IF_EXISTS;
476using ::remainder _LIBCPP_USING_IF_EXISTS;
477using ::remainderf _LIBCPP_USING_IF_EXISTS;
478using ::remquo _LIBCPP_USING_IF_EXISTS;
479using ::remquof _LIBCPP_USING_IF_EXISTS;
480using ::rint _LIBCPP_USING_IF_EXISTS;
481using ::rintf _LIBCPP_USING_IF_EXISTS;
482using ::round _LIBCPP_USING_IF_EXISTS;
483using ::roundf _LIBCPP_USING_IF_EXISTS;
484using ::scalbln _LIBCPP_USING_IF_EXISTS;
485using ::scalblnf _LIBCPP_USING_IF_EXISTS;
486using ::scalbn _LIBCPP_USING_IF_EXISTS;
487using ::scalbnf _LIBCPP_USING_IF_EXISTS;
488using ::tgamma _LIBCPP_USING_IF_EXISTS;
489using ::tgammaf _LIBCPP_USING_IF_EXISTS;
490using ::trunc _LIBCPP_USING_IF_EXISTS;
491using ::truncf _LIBCPP_USING_IF_EXISTS;
492
493using ::acosl _LIBCPP_USING_IF_EXISTS;
494using ::asinl _LIBCPP_USING_IF_EXISTS;
495using ::atanl _LIBCPP_USING_IF_EXISTS;
496using ::atan2l _LIBCPP_USING_IF_EXISTS;
497using ::ceill _LIBCPP_USING_IF_EXISTS;
498using ::cosl _LIBCPP_USING_IF_EXISTS;
499using ::coshl _LIBCPP_USING_IF_EXISTS;
500using ::expl _LIBCPP_USING_IF_EXISTS;
501using ::fabsl _LIBCPP_USING_IF_EXISTS;
502using ::floorl _LIBCPP_USING_IF_EXISTS;
503using ::fmodl _LIBCPP_USING_IF_EXISTS;
504using ::frexpl _LIBCPP_USING_IF_EXISTS;
505using ::ldexpl _LIBCPP_USING_IF_EXISTS;
506using ::logl _LIBCPP_USING_IF_EXISTS;
507using ::log10l _LIBCPP_USING_IF_EXISTS;
508using ::modfl _LIBCPP_USING_IF_EXISTS;
509using ::powl _LIBCPP_USING_IF_EXISTS;
510using ::sinl _LIBCPP_USING_IF_EXISTS;
511using ::sinhl _LIBCPP_USING_IF_EXISTS;
512using ::sqrtl _LIBCPP_USING_IF_EXISTS;
513using ::tanl _LIBCPP_USING_IF_EXISTS;
514
515using ::tanhl _LIBCPP_USING_IF_EXISTS;
516using ::acoshl _LIBCPP_USING_IF_EXISTS;
517using ::asinhl _LIBCPP_USING_IF_EXISTS;
518using ::atanhl _LIBCPP_USING_IF_EXISTS;
519using ::cbrtl _LIBCPP_USING_IF_EXISTS;
520
521using ::copysignl _LIBCPP_USING_IF_EXISTS;
522
523using ::erfl _LIBCPP_USING_IF_EXISTS;
524using ::erfcl _LIBCPP_USING_IF_EXISTS;
525using ::exp2l _LIBCPP_USING_IF_EXISTS;
526using ::expm1l _LIBCPP_USING_IF_EXISTS;
527using ::fdiml _LIBCPP_USING_IF_EXISTS;
528using ::fmal _LIBCPP_USING_IF_EXISTS;
529using ::fmaxl _LIBCPP_USING_IF_EXISTS;
530using ::fminl _LIBCPP_USING_IF_EXISTS;
531using ::hypotl _LIBCPP_USING_IF_EXISTS;
532using ::ilogbl _LIBCPP_USING_IF_EXISTS;
533using ::lgammal _LIBCPP_USING_IF_EXISTS;
534using ::llrintl _LIBCPP_USING_IF_EXISTS;
535using ::llroundl _LIBCPP_USING_IF_EXISTS;
536using ::log1pl _LIBCPP_USING_IF_EXISTS;
537using ::log2l _LIBCPP_USING_IF_EXISTS;
538using ::logbl _LIBCPP_USING_IF_EXISTS;
539using ::lrintl _LIBCPP_USING_IF_EXISTS;
540using ::lroundl _LIBCPP_USING_IF_EXISTS;
541using ::nanl _LIBCPP_USING_IF_EXISTS;
542using ::nearbyintl _LIBCPP_USING_IF_EXISTS;
543using ::nextafterl _LIBCPP_USING_IF_EXISTS;
544using ::nexttowardl _LIBCPP_USING_IF_EXISTS;
545using ::remainderl _LIBCPP_USING_IF_EXISTS;
546using ::remquol _LIBCPP_USING_IF_EXISTS;
547using ::rintl _LIBCPP_USING_IF_EXISTS;
548using ::roundl _LIBCPP_USING_IF_EXISTS;
549using ::scalblnl _LIBCPP_USING_IF_EXISTS;
550using ::scalbnl _LIBCPP_USING_IF_EXISTS;
551using ::tgammal _LIBCPP_USING_IF_EXISTS;
552using ::truncl _LIBCPP_USING_IF_EXISTS;
553
554template <class _A1, __enable_if_t<is_floating_point<_A1>::value, int> = 0>
555_LIBCPP_HIDE_FROM_ABI _LIBCPP_CONSTEXPR bool __constexpr_isinf(_A1 __lcpp_x) _NOEXCEPT {
556# if __has_builtin(__builtin_isinf)
557 return __builtin_isinf(__lcpp_x);
558# else
559 return isinf(__lcpp_x);
560# endif
561}
562
563template <class _A1, __enable_if_t<!is_floating_point<_A1>::value, int> = 0>
564_LIBCPP_HIDE_FROM_ABI _LIBCPP_CONSTEXPR bool __constexpr_isinf(_A1 __lcpp_x) _NOEXCEPT {
565 return std::isinf(__lcpp_x);
566}
567
568_LIBCPP_END_NAMESPACE_STD
569
570_LIBCPP_POP_MACROS
571
572#endif // __cplusplus < 201103L && defined(_LIBCPP_USE_FROZEN_CXX03_HEADERS)
573
574#endif // _LIBCPP_CMATH
575