ldouble_wrappers.c 7.7 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428
  1. /* vi: set sw=4 ts=4: */
  2. /*
  3. * Wrapper functions implementing all the long double math functions
  4. * defined by SuSv3 by actually calling the double version of
  5. * each function and then casting the result back to a long double
  6. * to return to the user.
  7. *
  8. * Copyright (C) 2005 by Erik Andersen <andersen@uclibc.org>
  9. *
  10. * Licensed under the LGPL v2.1, see the file COPYING.LIB in this tarball.
  11. */
  12. #include <features.h>
  13. /* Prevent math.h from defining colliding inlines */
  14. #undef __USE_EXTERN_INLINES
  15. #include "math.h"
  16. #include <complex.h>
  17. #define WRAPPER1(func) \
  18. long double func##l(long double x) \
  19. { \
  20. return (long double) func((double) x); \
  21. }
  22. #define int_WRAPPER1(func) \
  23. int func##l(long double x) \
  24. { \
  25. return func((double) x); \
  26. }
  27. #define long_WRAPPER1(func) \
  28. long func##l(long double x) \
  29. { \
  30. return func((double) x); \
  31. }
  32. #define long_long_WRAPPER1(func) \
  33. long long func##l(long double x) \
  34. { \
  35. return func((double) x); \
  36. }
  37. /* Implement the following, as defined by SuSv3 */
  38. #if 0
  39. long double acoshl(long double);
  40. long double acosl(long double);
  41. long double asinhl(long double);
  42. long double asinl(long double);
  43. long double atan2l(long double, long double);
  44. long double atanhl(long double);
  45. long double atanl(long double);
  46. long double cargl(long double complex);
  47. long double cbrtl(long double);
  48. long double ceill(long double);
  49. long double copysignl(long double, long double);
  50. long double coshl(long double);
  51. long double cosl(long double);
  52. long double erfcl(long double);
  53. long double erfl(long double);
  54. long double exp2l(long double);
  55. long double expl(long double);
  56. long double expm1l(long double);
  57. long double fabsl(long double);
  58. long double fdiml(long double, long double);
  59. long double floorl(long double);
  60. long double fmal(long double, long double, long double);
  61. long double fmaxl(long double, long double);
  62. long double fminl(long double, long double);
  63. long double fmodl(long double, long double);
  64. long double frexpl(long double value, int *);
  65. long double hypotl(long double, long double);
  66. int ilogbl(long double);
  67. long double ldexpl(long double, int);
  68. long double lgammal(long double);
  69. long long llrintl(long double);
  70. long long llroundl(long double);
  71. long double log10l(long double);
  72. long double log1pl(long double);
  73. long double log2l(long double);
  74. long double logbl(long double);
  75. long double logl(long double);
  76. long lrintl(long double);
  77. long lroundl(long double);
  78. long double modfl(long double, long double *);
  79. long double nearbyintl(long double);
  80. long double nextafterl(long double, long double);
  81. long double nexttowardl(long double, long double);
  82. long double powl(long double, long double);
  83. long double remainderl(long double, long double);
  84. long double remquol(long double, long double, int *);
  85. long double rintl(long double);
  86. long double roundl(long double);
  87. long double scalblnl(long double, long);
  88. long double scalbnl(long double, int);
  89. long double sinhl(long double);
  90. long double sinl(long double);
  91. long double sqrtl(long double);
  92. long double tanhl(long double);
  93. long double tanl(long double);
  94. long double tgammal(long double);
  95. long double truncl(long double);
  96. #endif
  97. #ifdef L_acoshl
  98. WRAPPER1(acosh)
  99. #endif
  100. #ifdef L_acosl
  101. WRAPPER1(acos)
  102. #endif
  103. #ifdef L_asinhl
  104. WRAPPER1(asinh)
  105. #endif
  106. #ifdef L_asinl
  107. WRAPPER1(asin)
  108. #endif
  109. #ifdef L_atan2l
  110. long double atan2l (long double x, long double y)
  111. {
  112. return (long double) atan2( (double)x, (double)y );
  113. }
  114. #endif
  115. #ifdef L_atanhl
  116. WRAPPER1(atanh)
  117. #endif
  118. #ifdef L_atanl
  119. WRAPPER1(atan)
  120. #endif
  121. #ifdef L_cargl
  122. long double cargl (long double complex x)
  123. {
  124. return (long double) carg( (double complex)x );
  125. }
  126. #endif
  127. #ifdef L_cbrtl
  128. WRAPPER1(cbrt)
  129. #endif
  130. #ifdef L_ceill
  131. WRAPPER1(ceil)
  132. #endif
  133. #ifdef L_copysignl
  134. long double copysignl (long double x, long double y)
  135. {
  136. return (long double) copysign( (double)x, (double)y );
  137. }
  138. #endif
  139. #ifdef L_coshl
  140. WRAPPER1(cosh)
  141. #endif
  142. #ifdef L_cosl
  143. WRAPPER1(cos)
  144. #endif
  145. #ifdef L_erfcl
  146. WRAPPER1(erfc)
  147. #endif
  148. #ifdef L_erfl
  149. WRAPPER1(erf)
  150. #endif
  151. #ifdef L_exp2l
  152. WRAPPER1(exp2)
  153. #endif
  154. #ifdef L_expl
  155. WRAPPER1(exp)
  156. #endif
  157. #ifdef L_expm1l
  158. WRAPPER1(expm1)
  159. #endif
  160. #ifdef L_fabsl
  161. WRAPPER1(fabs)
  162. #endif
  163. #ifdef L_fdiml
  164. long double fdiml (long double x, long double y)
  165. {
  166. return (long double) fdim( (double)x, (double)y );
  167. }
  168. #endif
  169. #ifdef L_floorl
  170. WRAPPER1(floor)
  171. #endif
  172. #ifdef L_fmal
  173. long double fmal (long double x, long double y, long double z)
  174. {
  175. return (long double) fma( (double)x, (double)y, (double)z );
  176. }
  177. #endif
  178. #ifdef L_fmaxl
  179. long double fmaxl (long double x, long double y)
  180. {
  181. return (long double) fmax( (double)x, (double)y );
  182. }
  183. #endif
  184. #ifdef L_fminl
  185. long double fminl (long double x, long double y)
  186. {
  187. return (long double) fmin( (double)x, (double)y );
  188. }
  189. #endif
  190. #ifdef L_fmodl
  191. long double fmodl (long double x, long double y)
  192. {
  193. return (long double) fmod( (double)x, (double)y );
  194. }
  195. #endif
  196. #ifdef L_frexpl
  197. long double frexpl (long double x, int *ex)
  198. {
  199. return (long double) frexp( (double)x, ex );
  200. }
  201. #endif
  202. #ifdef L_gammal
  203. WRAPPER1(gamma)
  204. #endif
  205. #ifdef L_hypotl
  206. long double hypotl (long double x, long double y)
  207. {
  208. return (long double) hypot( (double)x, (double)y );
  209. }
  210. #endif
  211. #ifdef L_ilogbl
  212. int_WRAPPER1(ilogb)
  213. #endif
  214. #ifdef L_ldexpl
  215. long double ldexpl (long double x, int ex)
  216. {
  217. return (long double) ldexp( (double)x, ex );
  218. }
  219. #endif
  220. #ifdef L_lgammal
  221. WRAPPER1(lgamma)
  222. #endif
  223. #ifdef L_llrintl
  224. long_long_WRAPPER1(llrint)
  225. #endif
  226. #ifdef L_llroundl
  227. long_long_WRAPPER1(llround)
  228. #endif
  229. #ifdef L_log10l
  230. WRAPPER1(log10)
  231. #endif
  232. #ifdef L_log1pl
  233. WRAPPER1(log1p)
  234. #endif
  235. #ifdef L_log2l
  236. WRAPPER1(log2)
  237. #endif
  238. #ifdef L_logbl
  239. WRAPPER1(logb)
  240. #endif
  241. #ifdef L_logl
  242. WRAPPER1(log)
  243. #endif
  244. #ifdef L_lrintl
  245. long_WRAPPER1(lrint)
  246. #endif
  247. #ifdef L_lroundl
  248. long_WRAPPER1(lround)
  249. #endif
  250. #ifdef L_modfl
  251. long double modfl (long double x, long double *iptr)
  252. {
  253. double y, result;
  254. result = modf ( x, &y );
  255. *iptr = (long double)y;
  256. return (long double) result;
  257. }
  258. #endif
  259. #ifdef L_nearbyintl
  260. WRAPPER1(nearbyint)
  261. #endif
  262. #ifdef L_nextafterl
  263. long double nextafterl (long double x, long double y)
  264. {
  265. return (long double) nextafter( (double)x, (double)y );
  266. }
  267. #endif
  268. /* Disabled in Makefile.in */
  269. #if 0 /* def L_nexttowardl */
  270. long double nexttowardl (long double x, long double y)
  271. {
  272. return (long double) nexttoward( (double)x, (double)y );
  273. }
  274. libm_hidden_def(nexttowardl)
  275. #endif
  276. #ifdef L_powl
  277. long double powl (long double x, long double y)
  278. {
  279. return (long double) pow( (double)x, (double)y );
  280. }
  281. #endif
  282. #ifdef L_remainderl
  283. long double remainderl (long double x, long double y)
  284. {
  285. return (long double) remainder( (double)x, (double)y );
  286. }
  287. #endif
  288. #ifdef L_remquol
  289. long double remquol (long double x, long double y, int *quo)
  290. {
  291. return (long double) remquo( (double)x, (double)y, quo );
  292. }
  293. #endif
  294. #ifdef L_rintl
  295. WRAPPER1(rint)
  296. #endif
  297. #ifdef L_roundl
  298. WRAPPER1(round)
  299. #endif
  300. #ifdef L_scalblnl
  301. long double scalblnl (long double x, long ex)
  302. {
  303. return (long double) scalbln( (double)x, ex );
  304. }
  305. #endif
  306. #ifdef L_scalbnl
  307. long double scalbnl (long double x, int ex)
  308. {
  309. return (long double) scalbn( (double)x, ex );
  310. }
  311. #endif
  312. /* scalb is an obsolete function */
  313. #ifdef L_sinhl
  314. WRAPPER1(sinh)
  315. #endif
  316. #ifdef L_sinl
  317. WRAPPER1(sin)
  318. #endif
  319. #ifdef L_sqrtl
  320. WRAPPER1(sqrt)
  321. #endif
  322. #ifdef L_tanhl
  323. WRAPPER1(tanh)
  324. #endif
  325. #ifdef L_tanl
  326. WRAPPER1(tan)
  327. #endif
  328. #ifdef L_tgammal
  329. WRAPPER1(tgamma)
  330. #endif
  331. #ifdef L_truncl
  332. WRAPPER1(trunc)
  333. #endif
  334. #ifdef L_significandl
  335. WRAPPER1(significand)
  336. #endif
  337. #if defined __DO_C99_MATH__ && !defined __NO_LONG_DOUBLE_MATH
  338. # ifdef L___fpclassifyl
  339. int_WRAPPER1(__fpclassify)
  340. libm_hidden_def(__fpclassifyl)
  341. # endif
  342. # ifdef L___finitel
  343. int_WRAPPER1(__finite)
  344. libm_hidden_def(__finitel)
  345. # endif
  346. # ifdef L___signbitl
  347. int_WRAPPER1(__signbit)
  348. libm_hidden_def(__signbitl)
  349. # endif
  350. # ifdef L___isnanl
  351. int_WRAPPER1(__isnan)
  352. libm_hidden_def(__isnanl)
  353. # endif
  354. # ifdef L___isinfl
  355. int_WRAPPER1(__isinf)
  356. libm_hidden_def(__isinfl)
  357. # endif
  358. #endif