math.h 11 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430
  1. #ifndef _MATH_H
  2. #define _MATH_H
  3. #ifdef __cplusplus
  4. extern "C" {
  5. #endif
  6. #include <features.h>
  7. #define __NEED_float_t
  8. #define __NEED_double_t
  9. #include <bits/alltypes.h>
  10. #if 100*__GNUC__+__GNUC_MINOR__ >= 303
  11. #define NAN __builtin_nanf("")
  12. #define INFINITY __builtin_inff()
  13. #else
  14. #define NAN (0.0f/0.0f)
  15. #define INFINITY 1e5000f
  16. #endif
  17. #define HUGE_VALF INFINITY
  18. #define HUGE_VAL ((double)INFINITY)
  19. #define HUGE_VALL ((long double)INFINITY)
  20. #define MATH_ERRNO 1
  21. #define MATH_ERREXCEPT 2
  22. #define math_errhandling 2
  23. #define FP_ILOGBNAN (-1-(int)(((unsigned)-1)>>1))
  24. #define FP_ILOGB0 FP_ILOGBNAN
  25. #define FP_NAN 0
  26. #define FP_INFINITE 1
  27. #define FP_ZERO 2
  28. #define FP_SUBNORMAL 3
  29. #define FP_NORMAL 4
  30. int __fpclassify(double);
  31. int __fpclassifyf(float);
  32. int __fpclassifyl(long double);
  33. static __inline unsigned __FLOAT_BITS(float __f)
  34. {
  35. union {float __f; unsigned __i;} __u;
  36. __u.__f = __f;
  37. return __u.__i;
  38. }
  39. static __inline unsigned long long __DOUBLE_BITS(double __f)
  40. {
  41. union {double __f; unsigned long long __i;} __u;
  42. __u.__f = __f;
  43. return __u.__i;
  44. }
  45. #define fpclassify(x) ( \
  46. sizeof(x) == sizeof(float) ? __fpclassifyf(x) : \
  47. sizeof(x) == sizeof(double) ? __fpclassify(x) : \
  48. __fpclassifyl(x) )
  49. #define isinf(x) ( \
  50. sizeof(x) == sizeof(float) ? (__FLOAT_BITS(x) & 0x7fffffff) == 0x7f800000 : \
  51. sizeof(x) == sizeof(double) ? (__DOUBLE_BITS(x) & -1ULL>>1) == 0x7ffULL<<52 : \
  52. __fpclassifyl(x) == FP_INFINITE)
  53. #define isnan(x) ( \
  54. sizeof(x) == sizeof(float) ? (__FLOAT_BITS(x) & 0x7fffffff) > 0x7f800000 : \
  55. sizeof(x) == sizeof(double) ? (__DOUBLE_BITS(x) & -1ULL>>1) > 0x7ffULL<<52 : \
  56. __fpclassifyl(x) == FP_NAN)
  57. #define isnormal(x) ( \
  58. sizeof(x) == sizeof(float) ? ((__FLOAT_BITS(x)+0x00800000) & 0x7fffffff) >= 0x01000000 : \
  59. sizeof(x) == sizeof(double) ? ((__DOUBLE_BITS(x)+(1ULL<<52)) & -1ULL>>1) >= 1ULL<<53 : \
  60. __fpclassifyl(x) == FP_NORMAL)
  61. #define isfinite(x) ( \
  62. sizeof(x) == sizeof(float) ? (__FLOAT_BITS(x) & 0x7fffffff) < 0x7f800000 : \
  63. sizeof(x) == sizeof(double) ? (__DOUBLE_BITS(x) & -1ULL>>1) < 0x7ffULL<<52 : \
  64. __fpclassifyl(x) > FP_INFINITE)
  65. int __signbit(double);
  66. int __signbitf(float);
  67. int __signbitl(long double);
  68. #define signbit(x) ( \
  69. sizeof(x) == sizeof(float) ? (int)(__FLOAT_BITS(x)>>31) : \
  70. sizeof(x) == sizeof(double) ? (int)(__DOUBLE_BITS(x)>>63) : \
  71. __signbitl(x) )
  72. #define isunordered(x,y) (isnan((x)) ? ((void)(y),1) : isnan((y)))
  73. #define __ISREL_DEF(rel, op, type) \
  74. static __inline int __is##rel(type __x, type __y) \
  75. { return !isunordered(__x,__y) && __x op __y; }
  76. __ISREL_DEF(lessf, <, float_t)
  77. __ISREL_DEF(less, <, double_t)
  78. __ISREL_DEF(lessl, <, long double)
  79. __ISREL_DEF(lessequalf, <=, float_t)
  80. __ISREL_DEF(lessequal, <=, double_t)
  81. __ISREL_DEF(lessequall, <=, long double)
  82. __ISREL_DEF(lessgreaterf, !=, float_t)
  83. __ISREL_DEF(lessgreater, !=, double_t)
  84. __ISREL_DEF(lessgreaterl, !=, long double)
  85. __ISREL_DEF(greaterf, >, float_t)
  86. __ISREL_DEF(greater, >, double_t)
  87. __ISREL_DEF(greaterl, >, long double)
  88. __ISREL_DEF(greaterequalf, >=, float_t)
  89. __ISREL_DEF(greaterequal, >=, double_t)
  90. __ISREL_DEF(greaterequall, >=, long double)
  91. #define __tg_pred_2(x, y, p) ( \
  92. sizeof((x)+(y)) == sizeof(float) ? p##f(x, y) : \
  93. sizeof((x)+(y)) == sizeof(double) ? p(x, y) : \
  94. p##l(x, y) )
  95. #define isless(x, y) __tg_pred_2(x, y, __isless)
  96. #define islessequal(x, y) __tg_pred_2(x, y, __islessequal)
  97. #define islessgreater(x, y) __tg_pred_2(x, y, __islessgreater)
  98. #define isgreater(x, y) __tg_pred_2(x, y, __isgreater)
  99. #define isgreaterequal(x, y) __tg_pred_2(x, y, __isgreaterequal)
  100. double acos(double);
  101. float acosf(float);
  102. long double acosl(long double);
  103. double acosh(double);
  104. float acoshf(float);
  105. long double acoshl(long double);
  106. double asin(double);
  107. float asinf(float);
  108. long double asinl(long double);
  109. double asinh(double);
  110. float asinhf(float);
  111. long double asinhl(long double);
  112. double atan(double);
  113. float atanf(float);
  114. long double atanl(long double);
  115. double atan2(double, double);
  116. float atan2f(float, float);
  117. long double atan2l(long double, long double);
  118. double atanh(double);
  119. float atanhf(float);
  120. long double atanhl(long double);
  121. double cbrt(double);
  122. float cbrtf(float);
  123. long double cbrtl(long double);
  124. double ceil(double);
  125. float ceilf(float);
  126. long double ceill(long double);
  127. double copysign(double, double);
  128. float copysignf(float, float);
  129. long double copysignl(long double, long double);
  130. double cos(double);
  131. float cosf(float);
  132. long double cosl(long double);
  133. double cosh(double);
  134. float coshf(float);
  135. long double coshl(long double);
  136. double erf(double);
  137. float erff(float);
  138. long double erfl(long double);
  139. double erfc(double);
  140. float erfcf(float);
  141. long double erfcl(long double);
  142. double exp(double);
  143. float expf(float);
  144. long double expl(long double);
  145. double exp2(double);
  146. float exp2f(float);
  147. long double exp2l(long double);
  148. double expm1(double);
  149. float expm1f(float);
  150. long double expm1l(long double);
  151. double fabs(double);
  152. float fabsf(float);
  153. long double fabsl(long double);
  154. double fdim(double, double);
  155. float fdimf(float, float);
  156. long double fdiml(long double, long double);
  157. double floor(double);
  158. float floorf(float);
  159. long double floorl(long double);
  160. double fma(double, double, double);
  161. float fmaf(float, float, float);
  162. long double fmal(long double, long double, long double);
  163. double fmax(double, double);
  164. float fmaxf(float, float);
  165. long double fmaxl(long double, long double);
  166. double fmin(double, double);
  167. float fminf(float, float);
  168. long double fminl(long double, long double);
  169. double fmod(double, double);
  170. float fmodf(float, float);
  171. long double fmodl(long double, long double);
  172. double frexp(double, int *);
  173. float frexpf(float, int *);
  174. long double frexpl(long double, int *);
  175. double hypot(double, double);
  176. float hypotf(float, float);
  177. long double hypotl(long double, long double);
  178. int ilogb(double);
  179. int ilogbf(float);
  180. int ilogbl(long double);
  181. double ldexp(double, int);
  182. float ldexpf(float, int);
  183. long double ldexpl(long double, int);
  184. double lgamma(double);
  185. float lgammaf(float);
  186. long double lgammal(long double);
  187. long long llrint(double);
  188. long long llrintf(float);
  189. long long llrintl(long double);
  190. long long llround(double);
  191. long long llroundf(float);
  192. long long llroundl(long double);
  193. double log(double);
  194. float logf(float);
  195. long double logl(long double);
  196. double log10(double);
  197. float log10f(float);
  198. long double log10l(long double);
  199. double log1p(double);
  200. float log1pf(float);
  201. long double log1pl(long double);
  202. double log2(double);
  203. float log2f(float);
  204. long double log2l(long double);
  205. double logb(double);
  206. float logbf(float);
  207. long double logbl(long double);
  208. long lrint(double);
  209. long lrintf(float);
  210. long lrintl(long double);
  211. long lround(double);
  212. long lroundf(float);
  213. long lroundl(long double);
  214. double modf(double, double *);
  215. float modff(float, float *);
  216. long double modfl(long double, long double *);
  217. double nan(const char *);
  218. float nanf(const char *);
  219. long double nanl(const char *);
  220. double nearbyint(double);
  221. float nearbyintf(float);
  222. long double nearbyintl(long double);
  223. double nextafter(double, double);
  224. float nextafterf(float, float);
  225. long double nextafterl(long double, long double);
  226. double nexttoward(double, long double);
  227. float nexttowardf(float, long double);
  228. long double nexttowardl(long double, long double);
  229. double pow(double, double);
  230. float powf(float, float);
  231. long double powl(long double, long double);
  232. double remainder(double, double);
  233. float remainderf(float, float);
  234. long double remainderl(long double, long double);
  235. double remquo(double, double, int *);
  236. float remquof(float, float, int *);
  237. long double remquol(long double, long double, int *);
  238. double rint(double);
  239. float rintf(float);
  240. long double rintl(long double);
  241. double round(double);
  242. float roundf(float);
  243. long double roundl(long double);
  244. double scalbln(double, long);
  245. float scalblnf(float, long);
  246. long double scalblnl(long double, long);
  247. double scalbn(double, int);
  248. float scalbnf(float, int);
  249. long double scalbnl(long double, int);
  250. double sin(double);
  251. float sinf(float);
  252. long double sinl(long double);
  253. double sinh(double);
  254. float sinhf(float);
  255. long double sinhl(long double);
  256. double sqrt(double);
  257. float sqrtf(float);
  258. long double sqrtl(long double);
  259. double tan(double);
  260. float tanf(float);
  261. long double tanl(long double);
  262. double tanh(double);
  263. float tanhf(float);
  264. long double tanhl(long double);
  265. double tgamma(double);
  266. float tgammaf(float);
  267. long double tgammal(long double);
  268. double trunc(double);
  269. float truncf(float);
  270. long double truncl(long double);
  271. #if defined(_XOPEN_SOURCE) || defined(_BSD_SOURCE)
  272. #undef MAXFLOAT
  273. #define MAXFLOAT 3.40282346638528859812e+38F
  274. #endif
  275. #if defined(_XOPEN_SOURCE) || defined(_GNU_SOURCE) || defined(_BSD_SOURCE)
  276. #define M_E 2.7182818284590452354 /* e */
  277. #define M_LOG2E 1.4426950408889634074 /* log_2 e */
  278. #define M_LOG10E 0.43429448190325182765 /* log_10 e */
  279. #define M_LN2 0.69314718055994530942 /* log_e 2 */
  280. #define M_LN10 2.30258509299404568402 /* log_e 10 */
  281. #define M_PI 3.14159265358979323846 /* pi */
  282. #define M_PI_2 1.57079632679489661923 /* pi/2 */
  283. #define M_PI_4 0.78539816339744830962 /* pi/4 */
  284. #define M_1_PI 0.31830988618379067154 /* 1/pi */
  285. #define M_2_PI 0.63661977236758134308 /* 2/pi */
  286. #define M_2_SQRTPI 1.12837916709551257390 /* 2/sqrt(pi) */
  287. #define M_SQRT2 1.41421356237309504880 /* sqrt(2) */
  288. #define M_SQRT1_2 0.70710678118654752440 /* 1/sqrt(2) */
  289. extern int signgam;
  290. double j0(double);
  291. double j1(double);
  292. double jn(int, double);
  293. double y0(double);
  294. double y1(double);
  295. double yn(int, double);
  296. #endif
  297. #if defined(_GNU_SOURCE) || defined(_BSD_SOURCE)
  298. #define HUGE 3.40282346638528859812e+38F
  299. double drem(double, double);
  300. float dremf(float, float);
  301. int finite(double);
  302. int finitef(float);
  303. double scalb(double, double);
  304. float scalbf(float, float);
  305. double significand(double);
  306. float significandf(float);
  307. double lgamma_r(double, int*);
  308. float lgammaf_r(float, int*);
  309. float j0f(float);
  310. float j1f(float);
  311. float jnf(int, float);
  312. float y0f(float);
  313. float y1f(float);
  314. float ynf(int, float);
  315. #endif
  316. #ifdef _GNU_SOURCE
  317. long double lgammal_r(long double, int*);
  318. void sincos(double, double*, double*);
  319. void sincosf(float, float*, float*);
  320. void sincosl(long double, long double*, long double*);
  321. double exp10(double);
  322. float exp10f(float);
  323. long double exp10l(long double);
  324. double pow10(double);
  325. float pow10f(float);
  326. long double pow10l(long double);
  327. #endif
  328. #ifdef __cplusplus
  329. }
  330. #endif
  331. #endif