math.h@ 6079

Last change on this file since 6079 was 5999, checked in by vboxsync, 17 years ago
The Giant CDDL Dual-License Header Change.
Property svn:eol-style set to `native`
File size: 27.2 KB

Line
1	/** @file
2	* innotek Portable Runtime / No-CRT - math.h.
3	*/
4
5	/*
6	* Copyright (C) 2006-2007 innotek GmbH
7	*
8	* This file is part of VirtualBox Open Source Edition (OSE), as
9	* available from http://www.virtualbox.org. This file is free software;
10	* you can redistribute it and/or modify it under the terms of the GNU
11	* General Public License (GPL) as published by the Free Software
12	* Foundation, in version 2 as it comes in the "COPYING" file of the
13	* VirtualBox OSE distribution. VirtualBox OSE is distributed in the
14	* hope that it will be useful, but WITHOUT ANY WARRANTY of any kind.
15	*
16	* The contents of this file may alternatively be used under the terms
17	* of the Common Development and Distribution License Version 1.0
18	* (CDDL) only, as it comes in the "COPYING.CDDL" file of the
19	* VirtualBox OSE distribution, in which case the provisions of the
20	* CDDL are applicable instead of those of the GPL.
21	*
22	* You may elect to license modified versions of this file under the
23	* terms and conditions of either the GPL or the CDDL or both.
24	* --------------------------------------------------------------------
25	*
26	* This code is based on:
27	*
28	* from: @(#)fdlibm.h 5.1 93/09/24
29	* $FreeBSD: src/lib/msun/src/math.h,v 1.61 2005/04/16 21:12:47 das Exp $
30	* FreeBSD HEAD 2005-06-xx
31	*
32	* ====================================================
33	* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
34	*
35	* Developed at SunPro, a Sun Microsystems, Inc. business.
36	* Permission to use, copy, modify, and distribute this
37	* software is freely granted, provided that this notice
38	* is preserved.
39	* ====================================================
40	*/
41
42	#ifndef ___iprt_nocrt_math_h
43	#define ___iprt_nocrt_math_h
44
45	#if !defined(__GNUC__) /* && !defined(__YOUR_COMPILER__) */
46	# error "IPRT: Adjust this header for your compiler"
47	#endif
48
49	#include <iprt/types.h>
50	/#include <machine/_limits.h>/
51
52	/* from sys/cdefs.h */
53	#if defined(__GNUC__) && !defined(__INTEL_COMPILER)
54	#define __GNUC_PREREQ__(ma, mi) \
55	(__GNUC__ > (ma) \|\| __GNUC__ == (ma) && __GNUC_MINOR__ >= (mi))
56	#else
57	#define __GNUC_PREREQ__(ma, mi) 0
58	#endif
59	#define __pure2
60
61
62	/*
63	* ANSI/POSIX
64	*/
65	extern const union __infinity_un {
66	unsigned char __uc[8];
67	double __ud;
68	} RT_NOCRT(__infinity);
69
70	extern const union __nan_un {
71	unsigned char __uc[sizeof(float)];
72	float __uf;
73	} RT_NOCRT(__nan);
74
75	#if __GNUC_PREREQ__(3, 3) \|\| (defined(__INTEL_COMPILER) && __INTEL_COMPILER >= 800)
76	#define __MATH_BUILTIN_CONSTANTS
77	#endif
78
79	#if __GNUC_PREREQ__(3, 0) && !defined(__INTEL_COMPILER)
80	#define __MATH_BUILTIN_RELOPS
81	#endif
82
83	#ifdef __MATH_BUILTIN_CONSTANTS
84	#define HUGE_VAL __builtin_huge_val()
85	#else
86	#define HUGE_VAL (RT_NOCRT(__infinity).__ud)
87	#endif
88
89	#if 1/* __ISO_C_VISIBLE >= 1999*/
90	#define FP_ILOGB0 (-__INT_MAX)
91	#define FP_ILOGBNAN __INT_MAX
92
93	#ifdef __MATH_BUILTIN_CONSTANTS
94	#define HUGE_VALF __builtin_huge_valf()
95	#define HUGE_VALL __builtin_huge_vall()
96	#define INFINITY __builtin_inf()
97	#define NAN __builtin_nan("")
98	#else
99	#define HUGE_VALF (float)HUGE_VAL
100	#define HUGE_VALL (long double)HUGE_VAL
101	#define INFINITY HUGE_VALF
102	#define NAN (__nan.__uf)
103	#endif /* __MATH_BUILTIN_CONSTANTS */
104
105	#define MATH_ERRNO 1
106	#define MATH_ERREXCEPT 2
107	#define math_errhandling MATH_ERREXCEPT
108
109	/* XXX We need a <machine/math.h>. */
110	#if defined(__ia64__) \|\| defined(__sparc64__)
111	#define FP_FAST_FMA
112	#endif
113	#ifdef __ia64__
114	#define FP_FAST_FMAL
115	#endif
116	#define FP_FAST_FMAF
117
118	/* Symbolic constants to classify floating point numbers. */
119	#define FP_INFINITE 0x01
120	#define FP_NAN 0x02
121	#define FP_NORMAL 0x04
122	#define FP_SUBNORMAL 0x08
123	#define FP_ZERO 0x10
124	#define fpclassify(x) \
125	((sizeof (x) == sizeof (float)) ? RT_NOCRT(__fpclassifyf)(x) \
126	: (sizeof (x) == sizeof (double)) ? RT_NOCRT(__fpclassifyd)(x) \
127	: RT_NOCRT(__fpclassifyl)(x))
128
129	#define isfinite(x) \
130	((sizeof (x) == sizeof (float)) ? RT_NOCRT(__isfinitef)(x) \
131	: (sizeof (x) == sizeof (double)) ? RT_NOCRT(__isfinite)(x) \
132	: RT_NOCRT(__isfinitel)(x))
133	#define isinf(x) \
134	((sizeof (x) == sizeof (float)) ? RT_NOCRT(__isinff)(x) \
135	: (sizeof (x) == sizeof (double)) ? isinf(x) \
136	: RT_NOCRT(__isinfl)(x))
137	#define isnan(x) \
138	((sizeof (x) == sizeof (float)) ? isnanf(x) \
139	: (sizeof (x) == sizeof (double)) ? isnan(x) \
140	: RT_NOCRT(__isnanl)(x))
141	#define isnormal(x) \
142	((sizeof (x) == sizeof (float)) ? RT_NOCRT(__isnormalf)(x) \
143	: (sizeof (x) == sizeof (double)) ? RT_NOCRT(__isnormal)(x) \
144	: RT_NOCRT(__isnormall)(x))
145
146	#ifdef __MATH_BUILTIN_RELOPS
147	#define isgreater(x, y) __builtin_isgreater((x), (y))
148	#define isgreaterequal(x, y) __builtin_isgreaterequal((x), (y))
149	#define isless(x, y) __builtin_isless((x), (y))
150	#define islessequal(x, y) __builtin_islessequal((x), (y))
151	#define islessgreater(x, y) __builtin_islessgreater((x), (y))
152	#define isunordered(x, y) __builtin_isunordered((x), (y))
153	#else
154	#define isgreater(x, y) (!isunordered((x), (y)) && (x) > (y))
155	#define isgreaterequal(x, y) (!isunordered((x), (y)) && (x) >= (y))
156	#define isless(x, y) (!isunordered((x), (y)) && (x) < (y))
157	#define islessequal(x, y) (!isunordered((x), (y)) && (x) <= (y))
158	#define islessgreater(x, y) (!isunordered((x), (y)) && \
159	((x) > (y) \|\| (y) > (x)))
160	#define isunordered(x, y) (isnan(x) \|\| isnan(y))
161	#endif /* __MATH_BUILTIN_RELOPS */
162
163	#define signbit(x) \
164	((sizeof (x) == sizeof (float)) ? RT_NOCRT(__signbitf)(x) \
165	: (sizeof (x) == sizeof (double)) ? RT_NOCRT(__signbit)(x) \
166	: RT_NOCRT(__signbitl)(x))
167
168	typedef double double_t;
169	typedef float float_t;
170	#endif /* __ISO_C_VISIBLE >= 1999 */
171
172	/*
173	* XOPEN/SVID
174	*/
175	#if 1/* __BSD_VISIBLE \|\| __XSI_VISIBLE*/
176	#define M_E 2.7182818284590452354 /* e */
177	#define M_LOG2E 1.4426950408889634074 /* log 2e */
178	#define M_LOG10E 0.43429448190325182765 /* log 10e */
179	#define M_LN2 0.69314718055994530942 /* log e2 */
180	#define M_LN10 2.30258509299404568402 /* log e10 */
181	#define M_PI 3.14159265358979323846 /* pi */
182	#define M_PI_2 1.57079632679489661923 /* pi/2 */
183	#define M_PI_4 0.78539816339744830962 /* pi/4 */
184	#define M_1_PI 0.31830988618379067154 /* 1/pi */
185	#define M_2_PI 0.63661977236758134308 /* 2/pi */
186	#define M_2_SQRTPI 1.12837916709551257390 /* 2/sqrt(pi) */
187	#define M_SQRT2 1.41421356237309504880 /* sqrt(2) */
188	#define M_SQRT1_2 0.70710678118654752440 /* 1/sqrt(2) */
189
190	#define MAXFLOAT ((float)3.40282346638528860e+38)
191	extern int RT_NOCRT(signgam);
192	#endif /* __BSD_VISIBLE \|\| __XSI_VISIBLE */
193
194	#if 1/* __BSD_VISIBLE*/
195	#if 0
196	/* Old value from 4.4BSD-Lite math.h; this is probably better. */
197	#define HUGE HUGE_VAL
198	#else
199	#define HUGE MAXFLOAT
200	#endif
201	#endif /* __BSD_VISIBLE */
202
203	/*
204	* Most of these functions depend on the rounding mode and have the side
205	* effect of raising floating-point exceptions, so they are not declared
206	* as __pure2. In C99, FENV_ACCESS affects the purity of these functions.
207	*/
208	__BEGIN_DECLS
209	/*
210	* ANSI/POSIX
211	*/
212	int RT_NOCRT(__fpclassifyd)(double) __pure2;
213	int RT_NOCRT(__fpclassifyf)(float) __pure2;
214	int RT_NOCRT(__fpclassifyl)(long double) __pure2;
215	int RT_NOCRT(__isfinitef)(float) __pure2;
216	int RT_NOCRT(__isfinite)(double) __pure2;
217	int RT_NOCRT(__isfinitel)(long double) __pure2;
218	int RT_NOCRT(__isinff)(float) __pure2;
219	int RT_NOCRT(__isinfl)(long double) __pure2;
220	int RT_NOCRT(__isnanl)(long double) __pure2;
221	int RT_NOCRT(__isnormalf)(float) __pure2;
222	int RT_NOCRT(__isnormal)(double) __pure2;
223	int RT_NOCRT(__isnormall)(long double) __pure2;
224	int RT_NOCRT(__signbit)(double) __pure2;
225	int RT_NOCRT(__signbitf)(float) __pure2;
226	int RT_NOCRT(__signbitl)(long double) __pure2;
227
228	double RT_NOCRT(acos)(double);
229	double RT_NOCRT(asin)(double);
230	double RT_NOCRT(atan)(double);
231	double RT_NOCRT(atan2)(double, double);
232	double RT_NOCRT(cos)(double);
233	double RT_NOCRT(sin)(double);
234	double RT_NOCRT(tan)(double);
235
236	double RT_NOCRT(cosh)(double);
237	double RT_NOCRT(sinh)(double);
238	double RT_NOCRT(tanh)(double);
239
240	double RT_NOCRT(exp)(double);
241	double RT_NOCRT(frexp)(double, int ); / fundamentally !__pure2 */
242	double RT_NOCRT(ldexp)(double, int);
243	double RT_NOCRT(log)(double);
244	double RT_NOCRT(log10)(double);
245	double RT_NOCRT(modf)(double, double ); / fundamentally !__pure2 */
246
247	double RT_NOCRT(pow)(double, double);
248	double RT_NOCRT(sqrt)(double);
249
250	double RT_NOCRT(ceil)(double);
251	double RT_NOCRT(fabs)(double) __pure2;
252	double RT_NOCRT(floor)(double);
253	double RT_NOCRT(fmod)(double, double);
254
255	/*
256	* These functions are not in C90.
257	*/
258	#if 1 /__BSD_VISIBLE \|\| __ISO_C_VISIBLE >= 1999 \|\| __XSI_VISIBLE/
259	double RT_NOCRT(acosh)(double);
260	double RT_NOCRT(asinh)(double);
261	double RT_NOCRT(atanh)(double);
262	double RT_NOCRT(cbrt)(double);
263	double RT_NOCRT(erf)(double);
264	double RT_NOCRT(erfc)(double);
265	double RT_NOCRT(exp2)(double);
266	double RT_NOCRT(expm1)(double);
267	double RT_NOCRT(fma)(double, double, double);
268	double RT_NOCRT(hypot)(double, double);
269	int RT_NOCRT(ilogb)(double) __pure2;
270	/int isinf(double) __pure2;/
271	/int isnan(double) __pure2;/
272	double RT_NOCRT(lgamma)(double);
273	long long RT_NOCRT(llrint)(double);
274	long long RT_NOCRT(llround)(double);
275	double RT_NOCRT(log1p)(double);
276	double RT_NOCRT(logb)(double);
277	long RT_NOCRT(lrint)(double);
278	long RT_NOCRT(lround)(double);
279	double RT_NOCRT(nextafter)(double, double);
280	double RT_NOCRT(remainder)(double, double);
281	double RT_NOCRT(remquo)(double, double, int *);
282	double RT_NOCRT(rint)(double);
283	#endif /* __BSD_VISIBLE \|\| __ISO_C_VISIBLE >= 1999 \|\| __XSI_VISIBLE */
284
285	#if 1/* __BSD_VISIBLE \|\| __XSI_VISIBLE*/
286	double RT_NOCRT(j0)(double);
287	double RT_NOCRT(j1)(double);
288	double RT_NOCRT(jn)(int, double);
289	double RT_NOCRT(scalb)(double, double);
290	double RT_NOCRT(y0)(double);
291	double RT_NOCRT(y1)(double);
292	double RT_NOCRT(yn)(int, double);
293
294	#if 1/* __XSI_VISIBLE <= 500 \|\| __BSD_VISIBLE*/
295	double RT_NOCRT(gamma)(double);
296	#endif
297	#endif /* __BSD_VISIBLE \|\| __XSI_VISIBLE */
298
299	#if 1/* __BSD_VISIBLE \|\| __ISO_C_VISIBLE >= 1999*/
300	double RT_NOCRT(copysign)(double, double) __pure2;
301	double RT_NOCRT(fdim)(double, double);
302	double RT_NOCRT(fmax)(double, double) __pure2;
303	double RT_NOCRT(fmin)(double, double) __pure2;
304	double RT_NOCRT(nearbyint)(double);
305	double RT_NOCRT(round)(double);
306	double RT_NOCRT(scalbln)(double, long);
307	double RT_NOCRT(scalbn)(double, int);
308	double RT_NOCRT(tgamma)(double);
309	double RT_NOCRT(trunc)(double);
310	#endif
311
312	/*
313	* BSD math library entry points
314	*/
315	#if 1/* __BSD_VISIBLE*/
316	double RT_NOCRT(drem)(double, double);
317	int RT_NOCRT(finite)(double) __pure2;
318	int RT_NOCRT(isnanf)(float) __pure2;
319
320	/*
321	* Reentrant version of gamma & lgamma; passes signgam back by reference
322	* as the second argument; user must allocate space for signgam.
323	*/
324	double RT_NOCRT(gamma_r)(double, int *);
325	double RT_NOCRT(lgamma_r)(double, int *);
326
327	/*
328	* IEEE Test Vector
329	*/
330	double RT_NOCRT(significand)(double);
331	#endif /* __BSD_VISIBLE */
332
333	/* float versions of ANSI/POSIX functions */
334	#if 1/* __ISO_C_VISIBLE >= 1999*/
335	float RT_NOCRT(acosf)(float);
336	float RT_NOCRT(asinf)(float);
337	float RT_NOCRT(atanf)(float);
338	float RT_NOCRT(atan2f)(float, float);
339	float RT_NOCRT(cosf)(float);
340	float RT_NOCRT(sinf)(float);
341	float RT_NOCRT(tanf)(float);
342
343	float RT_NOCRT(coshf)(float);
344	float RT_NOCRT(sinhf)(float);
345	float RT_NOCRT(tanhf)(float);
346
347	float RT_NOCRT(exp2f)(float);
348	float RT_NOCRT(expf)(float);
349	float RT_NOCRT(expm1f)(float);
350	float RT_NOCRT(frexpf)(float, int ); / fundamentally !__pure2 */
351	int RT_NOCRT(ilogbf)(float) __pure2;
352	float RT_NOCRT(ldexpf)(float, int);
353	float RT_NOCRT(log10f)(float);
354	float RT_NOCRT(log1pf)(float);
355	float RT_NOCRT(logf)(float);
356	float RT_NOCRT(modff)(float, float ); / fundamentally !__pure2 */
357
358	float RT_NOCRT(powf)(float, float);
359	float RT_NOCRT(sqrtf)(float);
360
361	float RT_NOCRT(ceilf)(float);
362	float RT_NOCRT(fabsf)(float) __pure2;
363	float RT_NOCRT(floorf)(float);
364	float RT_NOCRT(fmodf)(float, float);
365	float RT_NOCRT(roundf)(float);
366
367	float RT_NOCRT(erff)(float);
368	float RT_NOCRT(erfcf)(float);
369	float RT_NOCRT(hypotf)(float, float);
370	float RT_NOCRT(lgammaf)(float);
371
372	float RT_NOCRT(acoshf)(float);
373	float RT_NOCRT(asinhf)(float);
374	float RT_NOCRT(atanhf)(float);
375	float RT_NOCRT(cbrtf)(float);
376	float RT_NOCRT(logbf)(float);
377	float RT_NOCRT(copysignf)(float, float) __pure2;
378	long long RT_NOCRT(llrintf)(float);
379	long long RT_NOCRT(llroundf)(float);
380	long RT_NOCRT(lrintf)(float);
381	long RT_NOCRT(lroundf)(float);
382	float RT_NOCRT(nearbyintf)(float);
383	float RT_NOCRT(nextafterf)(float, float);
384	float RT_NOCRT(remainderf)(float, float);
385	float RT_NOCRT(remquof)(float, float, int *);
386	float RT_NOCRT(rintf)(float);
387	float RT_NOCRT(scalblnf)(float, long);
388	float RT_NOCRT(scalbnf)(float, int);
389	float RT_NOCRT(truncf)(float);
390
391	float RT_NOCRT(fdimf)(float, float);
392	float RT_NOCRT(fmaf)(float, float, float);
393	float RT_NOCRT(fmaxf)(float, float) __pure2;
394	float RT_NOCRT(fminf)(float, float) __pure2;
395	#endif
396
397	/*
398	* float versions of BSD math library entry points
399	*/
400	#if 1/* __BSD_VISIBLE*/
401	float RT_NOCRT(dremf)(float, float);
402	int RT_NOCRT(finitef)(float) __pure2;
403	float RT_NOCRT(gammaf)(float);
404	float RT_NOCRT(j0f)(float);
405	float RT_NOCRT(j1f)(float);
406	float RT_NOCRT(jnf)(int, float);
407	float RT_NOCRT(scalbf)(float, float);
408	float RT_NOCRT(y0f)(float);
409	float RT_NOCRT(y1f)(float);
410	float RT_NOCRT(ynf)(int, float);
411
412	/*
413	* Float versions of reentrant version of gamma & lgamma; passes
414	* signgam back by reference as the second argument; user must
415	* allocate space for signgam.
416	*/
417	float RT_NOCRT(gammaf_r)(float, int *);
418	float RT_NOCRT(lgammaf_r)(float, int *);
419
420	/*
421	* float version of IEEE Test Vector
422	*/
423	float RT_NOCRT(significandf)(float);
424	#endif /* __BSD_VISIBLE */
425
426	/*
427	* long double versions of ISO/POSIX math functions
428	*/
429	#if 1/* __ISO_C_VISIBLE >= 1999*/
430	#if 1 /* bird: we've got these */
431	long double RT_NOCRT(acoshl)(long double);
432	long double RT_NOCRT(acosl)(long double);
433	long double RT_NOCRT(asinhl)(long double);
434	long double RT_NOCRT(asinl)(long double);
435	long double RT_NOCRT(atan2l)(long double, long double);
436	long double RT_NOCRT(atanhl)(long double);
437	long double RT_NOCRT(atanl)(long double);
438	long double RT_NOCRT(cbrtl)(long double);
439	#endif
440	long double RT_NOCRT(ceill)(long double);
441	long double RT_NOCRT(copysignl)(long double, long double) __pure2;
442	#if 1 /* bird */
443	long double RT_NOCRT(coshl)(long double);
444	long double RT_NOCRT(cosl)(long double);
445	long double RT_NOCRT(erfcl)(long double);
446	long double RT_NOCRT(erfl)(long double);
447	long double RT_NOCRT(exp2l)(long double);
448	long double RT_NOCRT(expl)(long double);
449	long double RT_NOCRT(expm1l)(long double);
450	#endif
451	long double RT_NOCRT(fabsl)(long double) __pure2;
452	long double RT_NOCRT(fdiml)(long double, long double);
453	long double RT_NOCRT(floorl)(long double);
454	long double RT_NOCRT(fmal)(long double, long double, long double);
455	long double RT_NOCRT(fmaxl)(long double, long double) __pure2;
456	long double RT_NOCRT(fminl)(long double, long double) __pure2;
457	#if 1 /* bird */
458	long double RT_NOCRT(fmodl)(long double, long double);
459	#endif
460	long double RT_NOCRT(frexpl)(long double value, int ); / fundamentally !__pure2 */
461	#if 1 /* bird */
462	long double RT_NOCRT(hypotl)(long double, long double);
463	#endif
464	int RT_NOCRT(ilogbl)(long double) __pure2;
465	long double RT_NOCRT(ldexpl)(long double, int);
466	#if 1 /* bird */
467	long double RT_NOCRT(lgammal)(long double);
468	long long RT_NOCRT(llrintl)(long double);
469	#endif
470	long long RT_NOCRT(llroundl)(long double);
471	#if 1 /* bird */
472	long double RT_NOCRT(log10l)(long double);
473	long double RT_NOCRT(log1pl)(long double);
474	long double RT_NOCRT(log2l)(long double);
475	long double RT_NOCRT(logbl)(long double);
476	long double RT_NOCRT(logl)(long double);
477	long RT_NOCRT(lrintl)(long double);
478	#endif
479	long RT_NOCRT(lroundl)(long double);
480	#if 1 /* bird */
481	long double RT_NOCRT(modfl)(long double, long double ); / fundamentally !__pure2 */
482	long double RT_NOCRT(nanl)(const char *) __pure2;
483	long double RT_NOCRT(nearbyintl)(long double);
484	#endif
485	long double RT_NOCRT(nextafterl)(long double, long double);
486	double RT_NOCRT(nexttoward)(double, long double);
487	float RT_NOCRT(nexttowardf)(float, long double);
488	long double RT_NOCRT(nexttowardl)(long double, long double);
489	#if 1 /* bird */
490	long double RT_NOCRT(powl)(long double, long double);
491	long double RT_NOCRT(remainderl)(long double, long double);
492	long double RT_NOCRT(remquol)(long double, long double, int *);
493	long double RT_NOCRT(rintl)(long double);
494	#endif
495	long double RT_NOCRT(roundl)(long double);
496	long double RT_NOCRT(scalblnl)(long double, long);
497	long double RT_NOCRT(scalbnl)(long double, int);
498	#if 1 /* bird: we 've got most of these. */
499	long double RT_NOCRT(sinhl)(long double);
500	long double RT_NOCRT(sinl)(long double);
501	long double RT_NOCRT(sqrtl)(long double);
502	long double RT_NOCRT(tanhl)(long double);
503	long double RT_NOCRT(tanl)(long double);
504	long double RT_NOCRT(tgammal)(long double);
505	#endif
506	long double RT_NOCRT(truncl)(long double);
507
508	/* bird: these were missing, gcc apparently inlines them. */
509	double RT_NOCRT(nan)(const char *);
510	float RT_NOCRT(nanf)(const char *);
511
512	#endif /* __ISO_C_VISIBLE >= 1999 */
513
514	#if 1/def __USE_GNU/
515	/*
516	* In GLIBC there are long variants of the XOPEN/SVID constant
517	* block some pages ago. We need this to get the math tests going.
518	*/
519	#define M_El 2.7182818284590452353602874713526625L
520	#define M_LOG2El 1.4426950408889634073599246810018921L
521	#define M_LOG10El 0.4342944819032518276511289189166051L
522	#define M_LN2l 0.6931471805599453094172321214581766L
523	#define M_LN10l 2.3025850929940456840179914546843642L
524	#define M_PIl 3.1415926535897932384626433832795029L
525	#define M_PI_2l 1.5707963267948966192313216916397514L
526	#define M_PI_4l 0.7853981633974483096156608458198757L
527	#define M_1_PIl 0.3183098861837906715377675267450287L
528	#define M_2_PIl 0.6366197723675813430755350534900574L
529	#define M_2_SQRTPIl 1.1283791670955125738961589031215452L
530	#define M_SQRT2l 1.4142135623730950488016887242096981L
531	#define M_SQRT1_2l 0.7071067811865475244008443621048490L
532	#endif
533
534	#if 1/def __USE_GNU/
535
536	void RT_NOCRT(sincos)(double, double , double );
537	void RT_NOCRT(sincosf)(float, float , float );
538	void RT_NOCRT(sincosl)(long double, long double , long double );
539	float RT_NOCRT(exp10f)(float);
540	double RT_NOCRT(exp10)(double);
541	long double RT_NOCRT(exp10l)(long double);
542	float RT_NOCRT(log2f)(float);
543	double RT_NOCRT(log2)(double);
544	long double RT_NOCRT(log2l)(long double);
545	float RT_NOCRT(tgammaf)(float);
546	long double RT_NOCRT(significandl)(long double);
547	long double RT_NOCRT(j0l)(long double);
548	long double RT_NOCRT(j1l)(long double);
549	long double RT_NOCRT(jnl)(int, long double);
550	long double RT_NOCRT(scalbl)(long double, long double);
551	long double RT_NOCRT(y0l)(long double);
552	long double RT_NOCRT(y1l)(long double);
553	long double RT_NOCRT(ynl)(int, long double);
554	long double RT_NOCRT(lgammal_r)(long double,int *);
555	long double RT_NOCRT(gammal)(long double);
556	#endif
557	__END_DECLS
558
559
560
561	#ifndef RT_WITHOUT_NOCRT_WRAPPERS
562	/* sed -e "/#/d" -e "/RT_NOCRT/!d" -e "s/^.RT_NOCRT($[a-z0-9_]$).$/# define \1 RT_NOCRT(\1)/" /
563	# define __infinity RT_NOCRT(__infinity)
564	# define __nan RT_NOCRT(__nan)
565	# define __fpclassifyf RT_NOCRT(__fpclassifyf)
566	# define __fpclassifyd RT_NOCRT(__fpclassifyd)
567	# define __fpclassifyl RT_NOCRT(__fpclassifyl)
568	# define __isfinitef RT_NOCRT(__isfinitef)
569	# define __isfinite RT_NOCRT(__isfinite)
570	# define __isfinitel RT_NOCRT(__isfinitel)
571	# define __isinff RT_NOCRT(__isinff)
572	# define __isinfl RT_NOCRT(__isinfl)
573	# define __isnanl RT_NOCRT(__isnanl)
574	# define __isnormalf RT_NOCRT(__isnormalf)
575	# define __isnormal RT_NOCRT(__isnormal)
576	# define __isnormall RT_NOCRT(__isnormall)
577	# define __signbitf RT_NOCRT(__signbitf)
578	# define __signbit RT_NOCRT(__signbit)
579	# define __signbitl RT_NOCRT(__signbitl)
580	# define signgam RT_NOCRT(signgam)
581	# define __fpclassifyd RT_NOCRT(__fpclassifyd)
582	# define __fpclassifyf RT_NOCRT(__fpclassifyf)
583	# define __fpclassifyl RT_NOCRT(__fpclassifyl)
584	# define __isfinitef RT_NOCRT(__isfinitef)
585	# define __isfinite RT_NOCRT(__isfinite)
586	# define __isfinitel RT_NOCRT(__isfinitel)
587	# define __isinff RT_NOCRT(__isinff)
588	# define __isinfl RT_NOCRT(__isinfl)
589	# define __isnanl RT_NOCRT(__isnanl)
590	# define __isnormalf RT_NOCRT(__isnormalf)
591	# define __isnormal RT_NOCRT(__isnormal)
592	# define __isnormall RT_NOCRT(__isnormall)
593	# define __signbit RT_NOCRT(__signbit)
594	# define __signbitf RT_NOCRT(__signbitf)
595	# define __signbitl RT_NOCRT(__signbitl)
596	# define acos RT_NOCRT(acos)
597	# define asin RT_NOCRT(asin)
598	# define atan RT_NOCRT(atan)
599	# define atan2 RT_NOCRT(atan2)
600	# define cos RT_NOCRT(cos)
601	# define sin RT_NOCRT(sin)
602	# define tan RT_NOCRT(tan)
603	# define cosh RT_NOCRT(cosh)
604	# define sinh RT_NOCRT(sinh)
605	# define tanh RT_NOCRT(tanh)
606	# define exp RT_NOCRT(exp)
607	# define frexp RT_NOCRT(frexp)
608	# define ldexp RT_NOCRT(ldexp)
609	# define log RT_NOCRT(log)
610	# define log10 RT_NOCRT(log10)
611	# define modf RT_NOCRT(modf)
612	# define pow RT_NOCRT(pow)
613	# define sqrt RT_NOCRT(sqrt)
614	# define ceil RT_NOCRT(ceil)
615	# define fabs RT_NOCRT(fabs)
616	# define floor RT_NOCRT(floor)
617	# define fmod RT_NOCRT(fmod)
618	# define acosh RT_NOCRT(acosh)
619	# define asinh RT_NOCRT(asinh)
620	# define atanh RT_NOCRT(atanh)
621	# define cbrt RT_NOCRT(cbrt)
622	# define erf RT_NOCRT(erf)
623	# define erfc RT_NOCRT(erfc)
624	# define exp2 RT_NOCRT(exp2)
625	# define expm1 RT_NOCRT(expm1)
626	# define fma RT_NOCRT(fma)
627	# define hypot RT_NOCRT(hypot)
628	# define ilogb RT_NOCRT(ilogb)
629	# define lgamma RT_NOCRT(lgamma)
630	# define llrint RT_NOCRT(llrint)
631	# define llround RT_NOCRT(llround)
632	# define log1p RT_NOCRT(log1p)
633	# define logb RT_NOCRT(logb)
634	# define lrint RT_NOCRT(lrint)
635	# define lround RT_NOCRT(lround)
636	# define nextafter RT_NOCRT(nextafter)
637	# define remainder RT_NOCRT(remainder)
638	# define remquo RT_NOCRT(remquo)
639	# define rint RT_NOCRT(rint)
640	# define j0 RT_NOCRT(j0)
641	# define j1 RT_NOCRT(j1)
642	# define jn RT_NOCRT(jn)
643	# define scalb RT_NOCRT(scalb)
644	# define y0 RT_NOCRT(y0)
645	# define y1 RT_NOCRT(y1)
646	# define yn RT_NOCRT(yn)
647	# define gamma RT_NOCRT(gamma)
648	# define copysign RT_NOCRT(copysign)
649	# define fdim RT_NOCRT(fdim)
650	# define fmax RT_NOCRT(fmax)
651	# define fmin RT_NOCRT(fmin)
652	# define nearbyint RT_NOCRT(nearbyint)
653	# define round RT_NOCRT(round)
654	# define scalbln RT_NOCRT(scalbln)
655	# define scalbn RT_NOCRT(scalbn)
656	# define tgamma RT_NOCRT(tgamma)
657	# define trunc RT_NOCRT(trunc)
658	# define drem RT_NOCRT(drem)
659	# define finite RT_NOCRT(finite)
660	# define isnanf RT_NOCRT(isnanf)
661	# define gamma_r RT_NOCRT(gamma_r)
662	# define lgamma_r RT_NOCRT(lgamma_r)
663	# define significand RT_NOCRT(significand)
664	# define acosf RT_NOCRT(acosf)
665	# define asinf RT_NOCRT(asinf)
666	# define atanf RT_NOCRT(atanf)
667	# define atan2f RT_NOCRT(atan2f)
668	# define cosf RT_NOCRT(cosf)
669	# define sinf RT_NOCRT(sinf)
670	# define tanf RT_NOCRT(tanf)
671	# define coshf RT_NOCRT(coshf)
672	# define sinhf RT_NOCRT(sinhf)
673	# define tanhf RT_NOCRT(tanhf)
674	# define exp2f RT_NOCRT(exp2f)
675	# define expf RT_NOCRT(expf)
676	# define expm1f RT_NOCRT(expm1f)
677	# define frexpf RT_NOCRT(frexpf)
678	# define ilogbf RT_NOCRT(ilogbf)
679	# define ldexpf RT_NOCRT(ldexpf)
680	# define log10f RT_NOCRT(log10f)
681	# define log1pf RT_NOCRT(log1pf)
682	# define logf RT_NOCRT(logf)
683	# define modff RT_NOCRT(modff)
684	# define powf RT_NOCRT(powf)
685	# define sqrtf RT_NOCRT(sqrtf)
686	# define ceilf RT_NOCRT(ceilf)
687	# define fabsf RT_NOCRT(fabsf)
688	# define floorf RT_NOCRT(floorf)
689	# define fmodf RT_NOCRT(fmodf)
690	# define roundf RT_NOCRT(roundf)
691	# define erff RT_NOCRT(erff)
692	# define erfcf RT_NOCRT(erfcf)
693	# define hypotf RT_NOCRT(hypotf)
694	# define lgammaf RT_NOCRT(lgammaf)
695	# define acoshf RT_NOCRT(acoshf)
696	# define asinhf RT_NOCRT(asinhf)
697	# define atanhf RT_NOCRT(atanhf)
698	# define cbrtf RT_NOCRT(cbrtf)
699	# define logbf RT_NOCRT(logbf)
700	# define copysignf RT_NOCRT(copysignf)
701	# define llrintf RT_NOCRT(llrintf)
702	# define llroundf RT_NOCRT(llroundf)
703	# define lrintf RT_NOCRT(lrintf)
704	# define lroundf RT_NOCRT(lroundf)
705	# define nearbyintf RT_NOCRT(nearbyintf)
706	# define nextafterf RT_NOCRT(nextafterf)
707	# define remainderf RT_NOCRT(remainderf)
708	# define remquof RT_NOCRT(remquof)
709	# define rintf RT_NOCRT(rintf)
710	# define scalblnf RT_NOCRT(scalblnf)
711	# define scalbnf RT_NOCRT(scalbnf)
712	# define truncf RT_NOCRT(truncf)
713	# define fdimf RT_NOCRT(fdimf)
714	# define fmaf RT_NOCRT(fmaf)
715	# define fmaxf RT_NOCRT(fmaxf)
716	# define fminf RT_NOCRT(fminf)
717	# define dremf RT_NOCRT(dremf)
718	# define finitef RT_NOCRT(finitef)
719	# define gammaf RT_NOCRT(gammaf)
720	# define j0f RT_NOCRT(j0f)
721	# define j1f RT_NOCRT(j1f)
722	# define jnf RT_NOCRT(jnf)
723	# define scalbf RT_NOCRT(scalbf)
724	# define y0f RT_NOCRT(y0f)
725	# define y1f RT_NOCRT(y1f)
726	# define ynf RT_NOCRT(ynf)
727	# define gammaf_r RT_NOCRT(gammaf_r)
728	# define lgammaf_r RT_NOCRT(lgammaf_r)
729	# define significandf RT_NOCRT(significandf)
730	# define acoshl RT_NOCRT(acoshl)
731	# define acosl RT_NOCRT(acosl)
732	# define asinhl RT_NOCRT(asinhl)
733	# define asinl RT_NOCRT(asinl)
734	# define atan2l RT_NOCRT(atan2l)
735	# define atanhl RT_NOCRT(atanhl)
736	# define atanl RT_NOCRT(atanl)
737	# define cbrtl RT_NOCRT(cbrtl)
738	# define ceill RT_NOCRT(ceill)
739	# define copysignl RT_NOCRT(copysignl)
740	# define coshl RT_NOCRT(coshl)
741	# define cosl RT_NOCRT(cosl)
742	# define erfcl RT_NOCRT(erfcl)
743	# define erfl RT_NOCRT(erfl)
744	# define exp2l RT_NOCRT(exp2l)
745	# define expl RT_NOCRT(expl)
746	# define expm1l RT_NOCRT(expm1l)
747	# define fabsl RT_NOCRT(fabsl)
748	# define fdiml RT_NOCRT(fdiml)
749	# define floorl RT_NOCRT(floorl)
750	# define fmal RT_NOCRT(fmal)
751	# define fmaxl RT_NOCRT(fmaxl)
752	# define fminl RT_NOCRT(fminl)
753	# define fmodl RT_NOCRT(fmodl)
754	# define frexpl RT_NOCRT(frexpl)
755	# define hypotl RT_NOCRT(hypotl)
756	# define ilogbl RT_NOCRT(ilogbl)
757	# define ldexpl RT_NOCRT(ldexpl)
758	# define lgammal RT_NOCRT(lgammal)
759	# define llrintl RT_NOCRT(llrintl)
760	# define llroundl RT_NOCRT(llroundl)
761	# define log10l RT_NOCRT(log10l)
762	# define log1pl RT_NOCRT(log1pl)
763	# define log2l RT_NOCRT(log2l)
764	# define logbl RT_NOCRT(logbl)
765	# define logl RT_NOCRT(logl)
766	# define lrintl RT_NOCRT(lrintl)
767	# define lroundl RT_NOCRT(lroundl)
768	# define modfl RT_NOCRT(modfl)
769	# define nanl RT_NOCRT(nanl)
770	# define nearbyintl RT_NOCRT(nearbyintl)
771	# define nextafterl RT_NOCRT(nextafterl)
772	# define nexttoward RT_NOCRT(nexttoward)
773	# define nexttowardf RT_NOCRT(nexttowardf)
774	# define nexttowardl RT_NOCRT(nexttowardl)
775	# define powl RT_NOCRT(powl)
776	# define remainderl RT_NOCRT(remainderl)
777	# define remquol RT_NOCRT(remquol)
778	# define rintl RT_NOCRT(rintl)
779	# define roundl RT_NOCRT(roundl)
780	# define scalblnl RT_NOCRT(scalblnl)
781	# define scalbnl RT_NOCRT(scalbnl)
782	# define sinhl RT_NOCRT(sinhl)
783	# define sinl RT_NOCRT(sinl)
784	# define sqrtl RT_NOCRT(sqrtl)
785	# define tanhl RT_NOCRT(tanhl)
786	# define tanl RT_NOCRT(tanl)
787	# define tgammal RT_NOCRT(tgammal)
788	# define truncl RT_NOCRT(truncl)
789	# define nan RT_NOCRT(nan)
790	# define nanf RT_NOCRT(nanf)
791	# define sincos RT_NOCRT(sincos)
792	# define sincosf RT_NOCRT(sincosf)
793	# define sincosl RT_NOCRT(sincosl)
794	# define exp10f RT_NOCRT(exp10f)
795	# define exp10 RT_NOCRT(exp10)
796	# define exp10l RT_NOCRT(exp10l)
797	# define log2f RT_NOCRT(log2f)
798	# define log2 RT_NOCRT(log2)
799	# define log2l RT_NOCRT(log2l)
800	# define tgammaf RT_NOCRT(tgammaf)
801	# define significandl RT_NOCRT(significandl)
802	# define j0l RT_NOCRT(j0l)
803	# define j1l RT_NOCRT(j1l)
804	# define jnl RT_NOCRT(jnl)
805	# define scalbl RT_NOCRT(scalbl)
806	# define y0l RT_NOCRT(y0l)
807	# define y1l RT_NOCRT(y1l)
808	# define ynl RT_NOCRT(ynl)
809	# define lgammal_r RT_NOCRT(lgammal_r)
810	# define gammal RT_NOCRT(gammal)
811	#endif
812
813	/*
814	* Include inlined implementations.
815	*/
816	#ifdef RT_ARCH_AMD64
817	# include <iprt/nocrt/amd64/math.h>
818	#elif defined(RT_ARCH_X86)
819	# include <iprt/nocrt/x86/math.h>
820	#endif
821
822	#endif
823

Note: See TracBrowser for help on using the repository browser.

source: vbox/trunk/include/iprt/nocrt/math.h@ 6079

Download in other formats: