johnson 2.0.0.pre1 → 2.0.0.pre2

Sign up to get free protection for your applications and to get access to all the features.
Files changed (311) hide show
  1. data/CHANGELOG.rdoc +12 -0
  2. data/Manifest.txt +4 -285
  3. data/Rakefile +13 -20
  4. data/ext/tracemonkey/global.cc +4 -1
  5. data/ext/tracemonkey/js.cc +30 -33
  6. data/ext/tracemonkey/runtime.cc +96 -6
  7. data/ext/tracemonkey/split_global.cc +0 -5
  8. data/ext/tracemonkey/tracemonkey.h +2 -2
  9. data/lib/johnson.rb +2 -2
  10. data/lib/johnson/runtime.rb +25 -15
  11. data/lib/johnson/tracemonkey/runtime.rb +6 -3
  12. data/vendor/tracemonkey/config/system-headers +1 -3
  13. data/vendor/tracemonkey/jscntxt.h +5 -2
  14. data/vendor/tracemonkey/jsdbgapi.cpp +9 -1
  15. data/vendor/tracemonkey/jsdbgapi.h +4 -0
  16. data/vendor/tracemonkey/tests/ecma/NativeObjects/browser.js +0 -0
  17. data/vendor/tracemonkey/tests/ecma/NativeObjects/jstests.list +0 -0
  18. data/vendor/tracemonkey/tests/ecma_3_1/Object/jstests.list +1 -1
  19. data/vendor/tracemonkey/tests/js1_3/misc/browser.js +0 -0
  20. data/vendor/tracemonkey/tests/js1_3/misc/jstests.list +0 -0
  21. data/vendor/tracemonkey/tests/js1_5/Regress/jstests.list +4 -4
  22. data/vendor/tracemonkey/tests/js1_5/Scope/jstests.list +1 -1
  23. data/vendor/tracemonkey/tests/js1_5/decompilation/jstests.list +2 -2
  24. data/vendor/tracemonkey/tests/js1_7/decompilation/jstests.list +1 -1
  25. data/vendor/tracemonkey/tests/shell.js +2 -1
  26. metadata +51 -309
  27. data/ext/spidermonkey/context.c +0 -116
  28. data/ext/spidermonkey/context.h +0 -19
  29. data/ext/spidermonkey/conversions.c +0 -361
  30. data/ext/spidermonkey/conversions.h +0 -31
  31. data/ext/spidermonkey/debugger.c +0 -234
  32. data/ext/spidermonkey/debugger.h +0 -10
  33. data/ext/spidermonkey/extconf.rb +0 -32
  34. data/ext/spidermonkey/extensions.c +0 -37
  35. data/ext/spidermonkey/extensions.h +0 -12
  36. data/ext/spidermonkey/global.c +0 -40
  37. data/ext/spidermonkey/global.h +0 -11
  38. data/ext/spidermonkey/idhash.c +0 -16
  39. data/ext/spidermonkey/idhash.h +0 -8
  40. data/ext/spidermonkey/immutable_node.c +0 -1153
  41. data/ext/spidermonkey/immutable_node.c.erb +0 -523
  42. data/ext/spidermonkey/immutable_node.h +0 -22
  43. data/ext/spidermonkey/jroot.h +0 -197
  44. data/ext/spidermonkey/js_land_proxy.c +0 -620
  45. data/ext/spidermonkey/js_land_proxy.h +0 -20
  46. data/ext/spidermonkey/ruby_land_proxy.c +0 -618
  47. data/ext/spidermonkey/ruby_land_proxy.h +0 -38
  48. data/ext/spidermonkey/runtime.c +0 -396
  49. data/ext/spidermonkey/runtime.h +0 -27
  50. data/ext/spidermonkey/spidermonkey.c +0 -22
  51. data/ext/spidermonkey/spidermonkey.h +0 -29
  52. data/lib/johnson/spidermonkey.rb +0 -12
  53. data/lib/johnson/spidermonkey/context.rb +0 -10
  54. data/lib/johnson/spidermonkey/debugger.rb +0 -67
  55. data/lib/johnson/spidermonkey/immutable_node.rb +0 -282
  56. data/lib/johnson/spidermonkey/js_land_proxy.rb +0 -64
  57. data/lib/johnson/spidermonkey/mutable_tree_visitor.rb +0 -242
  58. data/lib/johnson/spidermonkey/ruby_land_proxy.rb +0 -17
  59. data/lib/johnson/spidermonkey/runtime.rb +0 -74
  60. data/test/johnson/spidermonkey/context_test.rb +0 -21
  61. data/test/johnson/spidermonkey/immutable_node_test.rb +0 -34
  62. data/test/johnson/spidermonkey/js_land_proxy_test.rb +0 -273
  63. data/test/johnson/spidermonkey/ruby_land_proxy_test.rb +0 -274
  64. data/test/johnson/spidermonkey/runtime_test.rb +0 -41
  65. data/vendor/spidermonkey/.cvsignore +0 -9
  66. data/vendor/spidermonkey/Makefile.in +0 -449
  67. data/vendor/spidermonkey/Makefile.ref +0 -365
  68. data/vendor/spidermonkey/README.html +0 -820
  69. data/vendor/spidermonkey/SpiderMonkey.rsp +0 -12
  70. data/vendor/spidermonkey/Y.js +0 -19
  71. data/vendor/spidermonkey/build.mk +0 -43
  72. data/vendor/spidermonkey/config.mk +0 -192
  73. data/vendor/spidermonkey/config/AIX4.1.mk +0 -65
  74. data/vendor/spidermonkey/config/AIX4.2.mk +0 -64
  75. data/vendor/spidermonkey/config/AIX4.3.mk +0 -65
  76. data/vendor/spidermonkey/config/Darwin.mk +0 -83
  77. data/vendor/spidermonkey/config/Darwin1.3.mk +0 -81
  78. data/vendor/spidermonkey/config/Darwin1.4.mk +0 -41
  79. data/vendor/spidermonkey/config/Darwin5.2.mk +0 -81
  80. data/vendor/spidermonkey/config/Darwin5.3.mk +0 -81
  81. data/vendor/spidermonkey/config/HP-UXB.10.10.mk +0 -77
  82. data/vendor/spidermonkey/config/HP-UXB.10.20.mk +0 -77
  83. data/vendor/spidermonkey/config/HP-UXB.11.00.mk +0 -80
  84. data/vendor/spidermonkey/config/IRIX.mk +0 -87
  85. data/vendor/spidermonkey/config/IRIX5.3.mk +0 -44
  86. data/vendor/spidermonkey/config/IRIX6.1.mk +0 -44
  87. data/vendor/spidermonkey/config/IRIX6.2.mk +0 -44
  88. data/vendor/spidermonkey/config/IRIX6.3.mk +0 -44
  89. data/vendor/spidermonkey/config/IRIX6.5.mk +0 -44
  90. data/vendor/spidermonkey/config/Linux_All.mk +0 -103
  91. data/vendor/spidermonkey/config/Mac_OS10.0.mk +0 -82
  92. data/vendor/spidermonkey/config/OSF1V4.0.mk +0 -72
  93. data/vendor/spidermonkey/config/OSF1V5.0.mk +0 -69
  94. data/vendor/spidermonkey/config/SunOS4.1.4.mk +0 -101
  95. data/vendor/spidermonkey/config/SunOS5.10.mk +0 -50
  96. data/vendor/spidermonkey/config/SunOS5.3.mk +0 -91
  97. data/vendor/spidermonkey/config/SunOS5.4.mk +0 -92
  98. data/vendor/spidermonkey/config/SunOS5.5.1.mk +0 -44
  99. data/vendor/spidermonkey/config/SunOS5.5.mk +0 -87
  100. data/vendor/spidermonkey/config/SunOS5.6.mk +0 -89
  101. data/vendor/spidermonkey/config/SunOS5.7.mk +0 -44
  102. data/vendor/spidermonkey/config/SunOS5.8.mk +0 -44
  103. data/vendor/spidermonkey/config/SunOS5.9.mk +0 -44
  104. data/vendor/spidermonkey/config/WINNT4.0.mk +0 -117
  105. data/vendor/spidermonkey/config/WINNT5.0.mk +0 -117
  106. data/vendor/spidermonkey/config/WINNT5.1.mk +0 -117
  107. data/vendor/spidermonkey/config/WINNT5.2.mk +0 -117
  108. data/vendor/spidermonkey/config/WINNT6.0.mk +0 -117
  109. data/vendor/spidermonkey/config/dgux.mk +0 -64
  110. data/vendor/spidermonkey/editline/Makefile.ref +0 -144
  111. data/vendor/spidermonkey/editline/README +0 -83
  112. data/vendor/spidermonkey/editline/editline.3 +0 -175
  113. data/vendor/spidermonkey/editline/editline.c +0 -1369
  114. data/vendor/spidermonkey/editline/editline.h +0 -135
  115. data/vendor/spidermonkey/editline/sysunix.c +0 -182
  116. data/vendor/spidermonkey/editline/unix.h +0 -82
  117. data/vendor/spidermonkey/fdlibm/.cvsignore +0 -7
  118. data/vendor/spidermonkey/fdlibm/Makefile.in +0 -127
  119. data/vendor/spidermonkey/fdlibm/Makefile.ref +0 -192
  120. data/vendor/spidermonkey/fdlibm/e_acos.c +0 -147
  121. data/vendor/spidermonkey/fdlibm/e_acosh.c +0 -105
  122. data/vendor/spidermonkey/fdlibm/e_asin.c +0 -156
  123. data/vendor/spidermonkey/fdlibm/e_atan2.c +0 -165
  124. data/vendor/spidermonkey/fdlibm/e_atanh.c +0 -110
  125. data/vendor/spidermonkey/fdlibm/e_cosh.c +0 -133
  126. data/vendor/spidermonkey/fdlibm/e_exp.c +0 -202
  127. data/vendor/spidermonkey/fdlibm/e_fmod.c +0 -184
  128. data/vendor/spidermonkey/fdlibm/e_gamma.c +0 -71
  129. data/vendor/spidermonkey/fdlibm/e_gamma_r.c +0 -70
  130. data/vendor/spidermonkey/fdlibm/e_hypot.c +0 -173
  131. data/vendor/spidermonkey/fdlibm/e_j0.c +0 -524
  132. data/vendor/spidermonkey/fdlibm/e_j1.c +0 -523
  133. data/vendor/spidermonkey/fdlibm/e_jn.c +0 -315
  134. data/vendor/spidermonkey/fdlibm/e_lgamma.c +0 -71
  135. data/vendor/spidermonkey/fdlibm/e_lgamma_r.c +0 -347
  136. data/vendor/spidermonkey/fdlibm/e_log.c +0 -184
  137. data/vendor/spidermonkey/fdlibm/e_log10.c +0 -134
  138. data/vendor/spidermonkey/fdlibm/e_pow.c +0 -386
  139. data/vendor/spidermonkey/fdlibm/e_rem_pio2.c +0 -222
  140. data/vendor/spidermonkey/fdlibm/e_remainder.c +0 -120
  141. data/vendor/spidermonkey/fdlibm/e_scalb.c +0 -89
  142. data/vendor/spidermonkey/fdlibm/e_sinh.c +0 -122
  143. data/vendor/spidermonkey/fdlibm/e_sqrt.c +0 -497
  144. data/vendor/spidermonkey/fdlibm/fdlibm.h +0 -273
  145. data/vendor/spidermonkey/fdlibm/fdlibm.mak +0 -1453
  146. data/vendor/spidermonkey/fdlibm/fdlibm.mdp +0 -0
  147. data/vendor/spidermonkey/fdlibm/k_cos.c +0 -135
  148. data/vendor/spidermonkey/fdlibm/k_rem_pio2.c +0 -354
  149. data/vendor/spidermonkey/fdlibm/k_sin.c +0 -114
  150. data/vendor/spidermonkey/fdlibm/k_standard.c +0 -785
  151. data/vendor/spidermonkey/fdlibm/k_tan.c +0 -170
  152. data/vendor/spidermonkey/fdlibm/s_asinh.c +0 -101
  153. data/vendor/spidermonkey/fdlibm/s_atan.c +0 -175
  154. data/vendor/spidermonkey/fdlibm/s_cbrt.c +0 -133
  155. data/vendor/spidermonkey/fdlibm/s_ceil.c +0 -120
  156. data/vendor/spidermonkey/fdlibm/s_copysign.c +0 -72
  157. data/vendor/spidermonkey/fdlibm/s_cos.c +0 -118
  158. data/vendor/spidermonkey/fdlibm/s_erf.c +0 -356
  159. data/vendor/spidermonkey/fdlibm/s_expm1.c +0 -267
  160. data/vendor/spidermonkey/fdlibm/s_fabs.c +0 -70
  161. data/vendor/spidermonkey/fdlibm/s_finite.c +0 -71
  162. data/vendor/spidermonkey/fdlibm/s_floor.c +0 -121
  163. data/vendor/spidermonkey/fdlibm/s_frexp.c +0 -99
  164. data/vendor/spidermonkey/fdlibm/s_ilogb.c +0 -85
  165. data/vendor/spidermonkey/fdlibm/s_isnan.c +0 -74
  166. data/vendor/spidermonkey/fdlibm/s_ldexp.c +0 -66
  167. data/vendor/spidermonkey/fdlibm/s_lib_version.c +0 -73
  168. data/vendor/spidermonkey/fdlibm/s_log1p.c +0 -211
  169. data/vendor/spidermonkey/fdlibm/s_logb.c +0 -79
  170. data/vendor/spidermonkey/fdlibm/s_matherr.c +0 -64
  171. data/vendor/spidermonkey/fdlibm/s_modf.c +0 -132
  172. data/vendor/spidermonkey/fdlibm/s_nextafter.c +0 -124
  173. data/vendor/spidermonkey/fdlibm/s_rint.c +0 -131
  174. data/vendor/spidermonkey/fdlibm/s_scalbn.c +0 -107
  175. data/vendor/spidermonkey/fdlibm/s_signgam.c +0 -40
  176. data/vendor/spidermonkey/fdlibm/s_significand.c +0 -68
  177. data/vendor/spidermonkey/fdlibm/s_sin.c +0 -118
  178. data/vendor/spidermonkey/fdlibm/s_tan.c +0 -112
  179. data/vendor/spidermonkey/fdlibm/s_tanh.c +0 -122
  180. data/vendor/spidermonkey/fdlibm/w_acos.c +0 -78
  181. data/vendor/spidermonkey/fdlibm/w_acosh.c +0 -78
  182. data/vendor/spidermonkey/fdlibm/w_asin.c +0 -80
  183. data/vendor/spidermonkey/fdlibm/w_atan2.c +0 -79
  184. data/vendor/spidermonkey/fdlibm/w_atanh.c +0 -81
  185. data/vendor/spidermonkey/fdlibm/w_cosh.c +0 -77
  186. data/vendor/spidermonkey/fdlibm/w_exp.c +0 -88
  187. data/vendor/spidermonkey/fdlibm/w_fmod.c +0 -78
  188. data/vendor/spidermonkey/fdlibm/w_gamma.c +0 -85
  189. data/vendor/spidermonkey/fdlibm/w_gamma_r.c +0 -81
  190. data/vendor/spidermonkey/fdlibm/w_hypot.c +0 -78
  191. data/vendor/spidermonkey/fdlibm/w_j0.c +0 -105
  192. data/vendor/spidermonkey/fdlibm/w_j1.c +0 -106
  193. data/vendor/spidermonkey/fdlibm/w_jn.c +0 -128
  194. data/vendor/spidermonkey/fdlibm/w_lgamma.c +0 -85
  195. data/vendor/spidermonkey/fdlibm/w_lgamma_r.c +0 -81
  196. data/vendor/spidermonkey/fdlibm/w_log.c +0 -78
  197. data/vendor/spidermonkey/fdlibm/w_log10.c +0 -81
  198. data/vendor/spidermonkey/fdlibm/w_pow.c +0 -99
  199. data/vendor/spidermonkey/fdlibm/w_remainder.c +0 -77
  200. data/vendor/spidermonkey/fdlibm/w_scalb.c +0 -95
  201. data/vendor/spidermonkey/fdlibm/w_sinh.c +0 -77
  202. data/vendor/spidermonkey/fdlibm/w_sqrt.c +0 -77
  203. data/vendor/spidermonkey/javascript-trace.d +0 -73
  204. data/vendor/spidermonkey/js.c +0 -3951
  205. data/vendor/spidermonkey/js.mdp +0 -0
  206. data/vendor/spidermonkey/js.msg +0 -308
  207. data/vendor/spidermonkey/js.pkg +0 -2
  208. data/vendor/spidermonkey/js3240.rc +0 -79
  209. data/vendor/spidermonkey/jsOS240.def +0 -654
  210. data/vendor/spidermonkey/jsapi.c +0 -5836
  211. data/vendor/spidermonkey/jsapi.h +0 -2624
  212. data/vendor/spidermonkey/jsarena.c +0 -450
  213. data/vendor/spidermonkey/jsarena.h +0 -318
  214. data/vendor/spidermonkey/jsarray.c +0 -2996
  215. data/vendor/spidermonkey/jsarray.h +0 -127
  216. data/vendor/spidermonkey/jsatom.c +0 -1045
  217. data/vendor/spidermonkey/jsatom.h +0 -442
  218. data/vendor/spidermonkey/jsbit.h +0 -253
  219. data/vendor/spidermonkey/jsbool.c +0 -176
  220. data/vendor/spidermonkey/jsbool.h +0 -73
  221. data/vendor/spidermonkey/jsclist.h +0 -139
  222. data/vendor/spidermonkey/jscntxt.c +0 -1348
  223. data/vendor/spidermonkey/jscntxt.h +0 -1120
  224. data/vendor/spidermonkey/jscompat.h +0 -57
  225. data/vendor/spidermonkey/jsconfig.h +0 -248
  226. data/vendor/spidermonkey/jsconfig.mk +0 -181
  227. data/vendor/spidermonkey/jscpucfg.c +0 -396
  228. data/vendor/spidermonkey/jscpucfg.h +0 -212
  229. data/vendor/spidermonkey/jsdate.c +0 -2390
  230. data/vendor/spidermonkey/jsdate.h +0 -124
  231. data/vendor/spidermonkey/jsdbgapi.c +0 -1802
  232. data/vendor/spidermonkey/jsdbgapi.h +0 -464
  233. data/vendor/spidermonkey/jsdhash.c +0 -868
  234. data/vendor/spidermonkey/jsdhash.h +0 -592
  235. data/vendor/spidermonkey/jsdtoa.c +0 -3167
  236. data/vendor/spidermonkey/jsdtoa.h +0 -130
  237. data/vendor/spidermonkey/jsdtracef.c +0 -317
  238. data/vendor/spidermonkey/jsdtracef.h +0 -77
  239. data/vendor/spidermonkey/jsemit.c +0 -6909
  240. data/vendor/spidermonkey/jsemit.h +0 -741
  241. data/vendor/spidermonkey/jsexn.c +0 -1371
  242. data/vendor/spidermonkey/jsexn.h +0 -96
  243. data/vendor/spidermonkey/jsfile.c +0 -2736
  244. data/vendor/spidermonkey/jsfile.h +0 -56
  245. data/vendor/spidermonkey/jsfile.msg +0 -90
  246. data/vendor/spidermonkey/jsfun.c +0 -2634
  247. data/vendor/spidermonkey/jsfun.h +0 -254
  248. data/vendor/spidermonkey/jsgc.c +0 -3562
  249. data/vendor/spidermonkey/jsgc.h +0 -403
  250. data/vendor/spidermonkey/jshash.c +0 -476
  251. data/vendor/spidermonkey/jshash.h +0 -151
  252. data/vendor/spidermonkey/jsify.pl +0 -485
  253. data/vendor/spidermonkey/jsinterp.c +0 -7007
  254. data/vendor/spidermonkey/jsinterp.h +0 -525
  255. data/vendor/spidermonkey/jsinvoke.c +0 -43
  256. data/vendor/spidermonkey/jsiter.c +0 -1067
  257. data/vendor/spidermonkey/jsiter.h +0 -122
  258. data/vendor/spidermonkey/jskeyword.tbl +0 -124
  259. data/vendor/spidermonkey/jskwgen.c +0 -460
  260. data/vendor/spidermonkey/jslibmath.h +0 -266
  261. data/vendor/spidermonkey/jslock.c +0 -1309
  262. data/vendor/spidermonkey/jslock.h +0 -313
  263. data/vendor/spidermonkey/jslocko.asm +0 -60
  264. data/vendor/spidermonkey/jslog2.c +0 -94
  265. data/vendor/spidermonkey/jslong.c +0 -264
  266. data/vendor/spidermonkey/jslong.h +0 -412
  267. data/vendor/spidermonkey/jsmath.c +0 -567
  268. data/vendor/spidermonkey/jsmath.h +0 -57
  269. data/vendor/spidermonkey/jsnum.c +0 -1239
  270. data/vendor/spidermonkey/jsnum.h +0 -283
  271. data/vendor/spidermonkey/jsobj.c +0 -5282
  272. data/vendor/spidermonkey/jsobj.h +0 -709
  273. data/vendor/spidermonkey/jsopcode.c +0 -5245
  274. data/vendor/spidermonkey/jsopcode.h +0 -394
  275. data/vendor/spidermonkey/jsopcode.tbl +0 -523
  276. data/vendor/spidermonkey/jsotypes.h +0 -202
  277. data/vendor/spidermonkey/jsparse.c +0 -6704
  278. data/vendor/spidermonkey/jsparse.h +0 -511
  279. data/vendor/spidermonkey/jsprf.c +0 -1264
  280. data/vendor/spidermonkey/jsprf.h +0 -150
  281. data/vendor/spidermonkey/jsproto.tbl +0 -128
  282. data/vendor/spidermonkey/jsprvtd.h +0 -267
  283. data/vendor/spidermonkey/jspubtd.h +0 -744
  284. data/vendor/spidermonkey/jsregexp.c +0 -4364
  285. data/vendor/spidermonkey/jsregexp.h +0 -183
  286. data/vendor/spidermonkey/jsreops.tbl +0 -145
  287. data/vendor/spidermonkey/jsscan.c +0 -2012
  288. data/vendor/spidermonkey/jsscan.h +0 -387
  289. data/vendor/spidermonkey/jsscope.c +0 -1957
  290. data/vendor/spidermonkey/jsscope.h +0 -418
  291. data/vendor/spidermonkey/jsscript.c +0 -1832
  292. data/vendor/spidermonkey/jsscript.h +0 -287
  293. data/vendor/spidermonkey/jsshell.msg +0 -50
  294. data/vendor/spidermonkey/jsstddef.h +0 -83
  295. data/vendor/spidermonkey/jsstr.c +0 -5005
  296. data/vendor/spidermonkey/jsstr.h +0 -641
  297. data/vendor/spidermonkey/jstypes.h +0 -475
  298. data/vendor/spidermonkey/jsutil.c +0 -345
  299. data/vendor/spidermonkey/jsutil.h +0 -157
  300. data/vendor/spidermonkey/jsxdrapi.c +0 -800
  301. data/vendor/spidermonkey/jsxdrapi.h +0 -218
  302. data/vendor/spidermonkey/jsxml.c +0 -8476
  303. data/vendor/spidermonkey/jsxml.h +0 -349
  304. data/vendor/spidermonkey/lock_SunOS.s +0 -119
  305. data/vendor/spidermonkey/perfect.js +0 -39
  306. data/vendor/spidermonkey/plify_jsdhash.sed +0 -36
  307. data/vendor/spidermonkey/prmjtime.c +0 -846
  308. data/vendor/spidermonkey/prmjtime.h +0 -103
  309. data/vendor/spidermonkey/resource.h +0 -15
  310. data/vendor/spidermonkey/rules.mk +0 -197
  311. data/vendor/spidermonkey/win32.order +0 -384
@@ -1,192 +0,0 @@
1
- # -*- Mode: makefile -*-
2
- #
3
- # ***** BEGIN LICENSE BLOCK *****
4
- # Version: MPL 1.1/GPL 2.0/LGPL 2.1
5
- #
6
- # The contents of this file are subject to the Mozilla Public License Version
7
- # 1.1 (the "License"); you may not use this file except in compliance with
8
- # the License. You may obtain a copy of the License at
9
- # http://www.mozilla.org/MPL/
10
- #
11
- # Software distributed under the License is distributed on an "AS IS" basis,
12
- # WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
13
- # for the specific language governing rights and limitations under the
14
- # License.
15
- #
16
- # The Original Code is Mozilla Communicator client code, released
17
- # March 31, 1998.
18
- #
19
- # The Initial Developer of the Original Code is
20
- # Sun Microsystems, Inc.
21
- # Portions created by the Initial Developer are Copyright (C) 1998
22
- # the Initial Developer. All Rights Reserved.
23
- #
24
- # Contributor(s):
25
- #
26
- # Alternatively, the contents of this file may be used under the terms of
27
- # either of the GNU General Public License Version 2 or later (the "GPL"),
28
- # or the GNU Lesser General Public License Version 2.1 or later (the "LGPL"),
29
- # in which case the provisions of the GPL or the LGPL are applicable instead
30
- # of those above. If you wish to allow use of your version of this file only
31
- # under the terms of either the GPL or the LGPL, and not to allow others to
32
- # use your version of this file under the terms of the MPL, indicate your
33
- # decision by deleting the provisions above and replace them with the notice
34
- # and other provisions required by the GPL or the LGPL. If you do not delete
35
- # the provisions above, a recipient may use your version of this file under
36
- # the terms of any one of the MPL, the GPL or the LGPL.
37
- #
38
- # ***** END LICENSE BLOCK *****
39
-
40
- #
41
- # @(#)Makefile 1.4 95/01/18
42
- #
43
- # ====================================================
44
- # Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
45
- #
46
- # Developed at SunSoft, a Sun Microsystems, Inc. business.
47
- # Permission to use, copy, modify, and distribute this
48
- # software is freely granted, provided that this notice
49
- # is preserved.
50
- # ====================================================
51
- #
52
- #
53
-
54
- #
55
- # There are two options in making libm at fdlibm compile time:
56
- # _IEEE_LIBM --- IEEE libm; smaller, and somewhat faster
57
- # _MULTI_LIBM --- Support multi-standard at runtime by
58
- # imposing wrapper functions defined in
59
- # fdlibm.h:
60
- # _IEEE_MODE -- IEEE
61
- # _XOPEN_MODE -- X/OPEN
62
- # _POSIX_MODE -- POSIX/ANSI
63
- # _SVID3_MODE -- SVID
64
- #
65
- # Here is how to set up CFLAGS to create the desired libm at
66
- # compile time:
67
- #
68
- # CFLAGS = -D_IEEE_LIBM ... IEEE libm (recommended)
69
- # CFLAGS = -D_SVID3_MODE ... Multi-standard supported
70
- # libm with SVID as the
71
- # default standard
72
- # CFLAGS = -D_XOPEN_MODE ... Multi-standard supported
73
- # libm with XOPEN as the
74
- # default standard
75
- # CFLAGS = -D_POSIX_MODE ... Multi-standard supported
76
- # libm with POSIX as the
77
- # default standard
78
- # CFLAGS = ... Multi-standard supported
79
- # libm with IEEE as the
80
- # default standard
81
- #
82
- # NOTE: if scalb's second arguement is an int, then one must
83
- # define _SCALB_INT in CFLAGS. The default prototype of scalb
84
- # is double scalb(double, double)
85
- #
86
-
87
- DEPTH = ..
88
-
89
- include $(DEPTH)/config.mk
90
-
91
- #
92
- # Default IEEE libm
93
- #
94
- CFLAGS += -DXP_UNIX $(OPTIMIZER) $(OS_CFLAGS) $(DEFINES) $(INCLUDES) \
95
- -DJSFILE $(XCFLAGS) -D_IEEE_LIBM
96
-
97
- # Need for jstypes.h and friends
98
- INCLUDES += -I..
99
- INCLUDES += -I../$(OBJDIR)
100
-
101
- #CC = cc
102
-
103
- INCFILES = fdlibm.h
104
- .INIT: $(INCFILES)
105
- .KEEP_STATE:
106
- FDLIBM_CFILES = \
107
- k_standard.c k_rem_pio2.c \
108
- k_cos.c k_sin.c k_tan.c \
109
- e_acos.c e_acosh.c e_asin.c e_atan2.c \
110
- e_atanh.c e_cosh.c e_exp.c e_fmod.c \
111
- e_gamma.c e_gamma_r.c e_hypot.c e_j0.c \
112
- e_j1.c e_jn.c e_lgamma.c e_lgamma_r.c \
113
- e_log.c e_log10.c e_pow.c e_rem_pio2.c e_remainder.c \
114
- e_scalb.c e_sinh.c e_sqrt.c \
115
- w_acos.c w_acosh.c w_asin.c w_atan2.c \
116
- w_atanh.c w_cosh.c w_exp.c w_fmod.c \
117
- w_gamma.c w_gamma_r.c w_hypot.c w_j0.c \
118
- w_j1.c w_jn.c w_lgamma.c w_lgamma_r.c \
119
- w_log.c w_log10.c w_pow.c w_remainder.c \
120
- w_scalb.c w_sinh.c w_sqrt.c \
121
- s_asinh.c s_atan.c s_cbrt.c s_ceil.c s_copysign.c \
122
- s_cos.c s_erf.c s_expm1.c s_fabs.c s_finite.c s_floor.c \
123
- s_frexp.c s_ilogb.c s_isnan.c s_ldexp.c s_lib_version.c \
124
- s_log1p.c s_logb.c s_matherr.c s_modf.c s_nextafter.c \
125
- s_rint.c s_scalbn.c s_signgam.c s_significand.c s_sin.c \
126
- s_tan.c s_tanh.c
127
-
128
- ifdef USE_MSVC
129
- FDLIBM_OBJS = $(addprefix $(OBJDIR)/, $(FDLIBM_CFILES:.c=.obj))
130
- else
131
- FDLIBM_OBJS = $(addprefix $(OBJDIR)/, $(FDLIBM_CFILES:.c=.o))
132
- endif
133
-
134
- ifdef USE_MSVC
135
- LIBRARY = $(OBJDIR)/fdlibm.lib
136
- else
137
- LIBRARY = $(OBJDIR)/libfdm.a
138
- endif
139
-
140
- define MAKE_OBJDIR
141
- if test ! -d $(@D); then rm -rf $(@D); mkdir -p $(@D); fi
142
- endef
143
-
144
- all: $(LIBRARY)
145
-
146
- export:
147
-
148
- $(OBJDIR)/%: %.c
149
- @$(MAKE_OBJDIR)
150
- $(CC) -o $@ $(CFLAGS) $*.c $(LDFLAGS)
151
-
152
- $(OBJDIR)/%.o: %.c
153
- @$(MAKE_OBJDIR)
154
- $(CC) -o $@ -c $(CFLAGS) $*.c
155
-
156
- $(OBJDIR)/%.o: %.s
157
- @$(MAKE_OBJDIR)
158
- $(AS) -o $@ $(ASFLAGS) $*.s
159
-
160
- # windows only
161
- $(OBJDIR)/%.obj: %.c
162
- @$(MAKE_OBJDIR)
163
- $(CC) -Fo$(OBJDIR)/ -c $(CFLAGS) $*.c
164
-
165
- ifeq ($(OS_ARCH),OS2)
166
- $(LIBRARY): $(FDLIBM_OBJS)
167
- $(AR) $@ $? $(AR_OS2_SUFFIX)
168
- $(RANLIB) $@
169
- else
170
- ifdef USE_MSVC
171
- $(LIBRARY): $(FDLIBM_OBJS)
172
- lib.exe /out:"$@" $?
173
- else
174
- $(LIBRARY): $(FDLIBM_OBJS)
175
- $(AR) rv $@ $?
176
- $(RANLIB) $@
177
- endif
178
- endif
179
-
180
- libfdm.a : $(FDLIBM_OBJS)
181
- $(AR) cru $(OBJDIR)/libfdm.a $(FDLIBM_OBJS)
182
- $(RANLIB) $(OBJDIR)/libfdm.a
183
-
184
- clean:
185
- rm -rf $(FDLIBM_OBJS)
186
-
187
- clobber:
188
- rm -rf $(FDLIBM_OBJS) $(LIBRARY) $(DEPENDENCIES)
189
-
190
- SUFFIXES: .i
191
- %.i: %.c
192
- $(CC) -C -E $(CFLAGS) $< > $*.i
@@ -1,147 +0,0 @@
1
- /* -*- Mode: C; tab-width: 8; indent-tabs-mode: nil; c-basic-offset: 4 -*-
2
- *
3
- * ***** BEGIN LICENSE BLOCK *****
4
- * Version: MPL 1.1/GPL 2.0/LGPL 2.1
5
- *
6
- * The contents of this file are subject to the Mozilla Public License Version
7
- * 1.1 (the "License"); you may not use this file except in compliance with
8
- * the License. You may obtain a copy of the License at
9
- * http://www.mozilla.org/MPL/
10
- *
11
- * Software distributed under the License is distributed on an "AS IS" basis,
12
- * WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
13
- * for the specific language governing rights and limitations under the
14
- * License.
15
- *
16
- * The Original Code is Mozilla Communicator client code, released
17
- * March 31, 1998.
18
- *
19
- * The Initial Developer of the Original Code is
20
- * Sun Microsystems, Inc.
21
- * Portions created by the Initial Developer are Copyright (C) 1998
22
- * the Initial Developer. All Rights Reserved.
23
- *
24
- * Contributor(s):
25
- *
26
- * Alternatively, the contents of this file may be used under the terms of
27
- * either of the GNU General Public License Version 2 or later (the "GPL"),
28
- * or the GNU Lesser General Public License Version 2.1 or later (the "LGPL"),
29
- * in which case the provisions of the GPL or the LGPL are applicable instead
30
- * of those above. If you wish to allow use of your version of this file only
31
- * under the terms of either the GPL or the LGPL, and not to allow others to
32
- * use your version of this file under the terms of the MPL, indicate your
33
- * decision by deleting the provisions above and replace them with the notice
34
- * and other provisions required by the GPL or the LGPL. If you do not delete
35
- * the provisions above, a recipient may use your version of this file under
36
- * the terms of any one of the MPL, the GPL or the LGPL.
37
- *
38
- * ***** END LICENSE BLOCK ***** */
39
-
40
- /* @(#)e_acos.c 1.3 95/01/18 */
41
- /*
42
- * ====================================================
43
- * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
44
- *
45
- * Developed at SunSoft, a Sun Microsystems, Inc. business.
46
- * Permission to use, copy, modify, and distribute this
47
- * software is freely granted, provided that this notice
48
- * is preserved.
49
- * ====================================================
50
- */
51
-
52
- /* __ieee754_acos(x)
53
- * Method :
54
- * acos(x) = pi/2 - asin(x)
55
- * acos(-x) = pi/2 + asin(x)
56
- * For |x|<=0.5
57
- * acos(x) = pi/2 - (x + x*x^2*R(x^2)) (see asin.c)
58
- * For x>0.5
59
- * acos(x) = pi/2 - (pi/2 - 2asin(sqrt((1-x)/2)))
60
- * = 2asin(sqrt((1-x)/2))
61
- * = 2s + 2s*z*R(z) ...z=(1-x)/2, s=sqrt(z)
62
- * = 2f + (2c + 2s*z*R(z))
63
- * where f=hi part of s, and c = (z-f*f)/(s+f) is the correction term
64
- * for f so that f+c ~ sqrt(z).
65
- * For x<-0.5
66
- * acos(x) = pi - 2asin(sqrt((1-|x|)/2))
67
- * = pi - 0.5*(s+s*z*R(z)), where z=(1-|x|)/2,s=sqrt(z)
68
- *
69
- * Special cases:
70
- * if x is NaN, return x itself;
71
- * if |x|>1, return NaN with invalid signal.
72
- *
73
- * Function needed: sqrt
74
- */
75
-
76
- #include "fdlibm.h"
77
-
78
- #ifdef __STDC__
79
- static const double
80
- #else
81
- static double
82
- #endif
83
- one= 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */
84
- pi = 3.14159265358979311600e+00, /* 0x400921FB, 0x54442D18 */
85
- pio2_hi = 1.57079632679489655800e+00, /* 0x3FF921FB, 0x54442D18 */
86
- pio2_lo = 6.12323399573676603587e-17, /* 0x3C91A626, 0x33145C07 */
87
- pS0 = 1.66666666666666657415e-01, /* 0x3FC55555, 0x55555555 */
88
- pS1 = -3.25565818622400915405e-01, /* 0xBFD4D612, 0x03EB6F7D */
89
- pS2 = 2.01212532134862925881e-01, /* 0x3FC9C155, 0x0E884455 */
90
- pS3 = -4.00555345006794114027e-02, /* 0xBFA48228, 0xB5688F3B */
91
- pS4 = 7.91534994289814532176e-04, /* 0x3F49EFE0, 0x7501B288 */
92
- pS5 = 3.47933107596021167570e-05, /* 0x3F023DE1, 0x0DFDF709 */
93
- qS1 = -2.40339491173441421878e+00, /* 0xC0033A27, 0x1C8A2D4B */
94
- qS2 = 2.02094576023350569471e+00, /* 0x40002AE5, 0x9C598AC8 */
95
- qS3 = -6.88283971605453293030e-01, /* 0xBFE6066C, 0x1B8D0159 */
96
- qS4 = 7.70381505559019352791e-02; /* 0x3FB3B8C5, 0xB12E9282 */
97
-
98
- #ifdef __STDC__
99
- double __ieee754_acos(double x)
100
- #else
101
- double __ieee754_acos(x)
102
- double x;
103
- #endif
104
- {
105
- fd_twoints u;
106
- double df;
107
- double z,p,q,r,w,s,c;
108
- int hx,ix;
109
- u.d = x;
110
- hx = __HI(u);
111
- ix = hx&0x7fffffff;
112
- if(ix>=0x3ff00000) { /* |x| >= 1 */
113
- if(((ix-0x3ff00000)|__LO(u))==0) { /* |x|==1 */
114
- if(hx>0) return 0.0; /* acos(1) = 0 */
115
- else return pi+2.0*pio2_lo; /* acos(-1)= pi */
116
- }
117
- return (x-x)/(x-x); /* acos(|x|>1) is NaN */
118
- }
119
- if(ix<0x3fe00000) { /* |x| < 0.5 */
120
- if(ix<=0x3c600000) return pio2_hi+pio2_lo;/*if|x|<2**-57*/
121
- z = x*x;
122
- p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5)))));
123
- q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4)));
124
- r = p/q;
125
- return pio2_hi - (x - (pio2_lo-x*r));
126
- } else if (hx<0) { /* x < -0.5 */
127
- z = (one+x)*0.5;
128
- p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5)))));
129
- q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4)));
130
- s = fd_sqrt(z);
131
- r = p/q;
132
- w = r*s-pio2_lo;
133
- return pi - 2.0*(s+w);
134
- } else { /* x > 0.5 */
135
- z = (one-x)*0.5;
136
- s = fd_sqrt(z);
137
- u.d = s;
138
- __LO(u) = 0;
139
- df = u.d;
140
- c = (z-df*df)/(s+df);
141
- p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5)))));
142
- q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4)));
143
- r = p/q;
144
- w = r*s+c;
145
- return 2.0*(df+w);
146
- }
147
- }
@@ -1,105 +0,0 @@
1
- /* -*- Mode: C; tab-width: 8; indent-tabs-mode: nil; c-basic-offset: 4 -*-
2
- *
3
- * ***** BEGIN LICENSE BLOCK *****
4
- * Version: MPL 1.1/GPL 2.0/LGPL 2.1
5
- *
6
- * The contents of this file are subject to the Mozilla Public License Version
7
- * 1.1 (the "License"); you may not use this file except in compliance with
8
- * the License. You may obtain a copy of the License at
9
- * http://www.mozilla.org/MPL/
10
- *
11
- * Software distributed under the License is distributed on an "AS IS" basis,
12
- * WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
13
- * for the specific language governing rights and limitations under the
14
- * License.
15
- *
16
- * The Original Code is Mozilla Communicator client code, released
17
- * March 31, 1998.
18
- *
19
- * The Initial Developer of the Original Code is
20
- * Sun Microsystems, Inc.
21
- * Portions created by the Initial Developer are Copyright (C) 1998
22
- * the Initial Developer. All Rights Reserved.
23
- *
24
- * Contributor(s):
25
- *
26
- * Alternatively, the contents of this file may be used under the terms of
27
- * either of the GNU General Public License Version 2 or later (the "GPL"),
28
- * or the GNU Lesser General Public License Version 2.1 or later (the "LGPL"),
29
- * in which case the provisions of the GPL or the LGPL are applicable instead
30
- * of those above. If you wish to allow use of your version of this file only
31
- * under the terms of either the GPL or the LGPL, and not to allow others to
32
- * use your version of this file under the terms of the MPL, indicate your
33
- * decision by deleting the provisions above and replace them with the notice
34
- * and other provisions required by the GPL or the LGPL. If you do not delete
35
- * the provisions above, a recipient may use your version of this file under
36
- * the terms of any one of the MPL, the GPL or the LGPL.
37
- *
38
- * ***** END LICENSE BLOCK ***** */
39
-
40
- /* @(#)e_acosh.c 1.3 95/01/18 */
41
- /*
42
- * ====================================================
43
- * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
44
- *
45
- * Developed at SunSoft, a Sun Microsystems, Inc. business.
46
- * Permission to use, copy, modify, and distribute this
47
- * software is freely granted, provided that this notice
48
- * is preserved.
49
- * ====================================================
50
- *
51
- */
52
-
53
- /* __ieee754_acosh(x)
54
- * Method :
55
- * Based on
56
- * acosh(x) = log [ x + sqrt(x*x-1) ]
57
- * we have
58
- * acosh(x) := log(x)+ln2, if x is large; else
59
- * acosh(x) := log(2x-1/(sqrt(x*x-1)+x)) if x>2; else
60
- * acosh(x) := log1p(t+sqrt(2.0*t+t*t)); where t=x-1.
61
- *
62
- * Special cases:
63
- * acosh(x) is NaN with signal if x<1.
64
- * acosh(NaN) is NaN without signal.
65
- */
66
-
67
- #include "fdlibm.h"
68
-
69
- #ifdef __STDC__
70
- static const double
71
- #else
72
- static double
73
- #endif
74
- one = 1.0,
75
- ln2 = 6.93147180559945286227e-01; /* 0x3FE62E42, 0xFEFA39EF */
76
-
77
- #ifdef __STDC__
78
- double __ieee754_acosh(double x)
79
- #else
80
- double __ieee754_acosh(x)
81
- double x;
82
- #endif
83
- {
84
- fd_twoints u;
85
- double t;
86
- int hx;
87
- u.d = x;
88
- hx = __HI(u);
89
- if(hx<0x3ff00000) { /* x < 1 */
90
- return (x-x)/(x-x);
91
- } else if(hx >=0x41b00000) { /* x > 2**28 */
92
- if(hx >=0x7ff00000) { /* x is inf of NaN */
93
- return x+x;
94
- } else
95
- return __ieee754_log(x)+ln2; /* acosh(huge)=log(2x) */
96
- } else if(((hx-0x3ff00000)|__LO(u))==0) {
97
- return 0.0; /* acosh(1) = 0 */
98
- } else if (hx > 0x40000000) { /* 2**28 > x > 2 */
99
- t=x*x;
100
- return __ieee754_log(2.0*x-one/(x+fd_sqrt(t-one)));
101
- } else { /* 1<x<2 */
102
- t = x-one;
103
- return fd_log1p(t+fd_sqrt(2.0*t+t*t));
104
- }
105
- }
@@ -1,156 +0,0 @@
1
- /* -*- Mode: C; tab-width: 8; indent-tabs-mode: nil; c-basic-offset: 4 -*-
2
- *
3
- * ***** BEGIN LICENSE BLOCK *****
4
- * Version: MPL 1.1/GPL 2.0/LGPL 2.1
5
- *
6
- * The contents of this file are subject to the Mozilla Public License Version
7
- * 1.1 (the "License"); you may not use this file except in compliance with
8
- * the License. You may obtain a copy of the License at
9
- * http://www.mozilla.org/MPL/
10
- *
11
- * Software distributed under the License is distributed on an "AS IS" basis,
12
- * WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
13
- * for the specific language governing rights and limitations under the
14
- * License.
15
- *
16
- * The Original Code is Mozilla Communicator client code, released
17
- * March 31, 1998.
18
- *
19
- * The Initial Developer of the Original Code is
20
- * Sun Microsystems, Inc.
21
- * Portions created by the Initial Developer are Copyright (C) 1998
22
- * the Initial Developer. All Rights Reserved.
23
- *
24
- * Contributor(s):
25
- *
26
- * Alternatively, the contents of this file may be used under the terms of
27
- * either of the GNU General Public License Version 2 or later (the "GPL"),
28
- * or the GNU Lesser General Public License Version 2.1 or later (the "LGPL"),
29
- * in which case the provisions of the GPL or the LGPL are applicable instead
30
- * of those above. If you wish to allow use of your version of this file only
31
- * under the terms of either the GPL or the LGPL, and not to allow others to
32
- * use your version of this file under the terms of the MPL, indicate your
33
- * decision by deleting the provisions above and replace them with the notice
34
- * and other provisions required by the GPL or the LGPL. If you do not delete
35
- * the provisions above, a recipient may use your version of this file under
36
- * the terms of any one of the MPL, the GPL or the LGPL.
37
- *
38
- * ***** END LICENSE BLOCK ***** */
39
-
40
- /* @(#)e_asin.c 1.3 95/01/18 */
41
- /*
42
- * ====================================================
43
- * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
44
- *
45
- * Developed at SunSoft, a Sun Microsystems, Inc. business.
46
- * Permission to use, copy, modify, and distribute this
47
- * software is freely granted, provided that this notice
48
- * is preserved.
49
- * ====================================================
50
- */
51
-
52
- /* __ieee754_asin(x)
53
- * Method :
54
- * Since asin(x) = x + x^3/6 + x^5*3/40 + x^7*15/336 + ...
55
- * we approximate asin(x) on [0,0.5] by
56
- * asin(x) = x + x*x^2*R(x^2)
57
- * where
58
- * R(x^2) is a rational approximation of (asin(x)-x)/x^3
59
- * and its remez error is bounded by
60
- * |(asin(x)-x)/x^3 - R(x^2)| < 2^(-58.75)
61
- *
62
- * For x in [0.5,1]
63
- * asin(x) = pi/2-2*asin(sqrt((1-x)/2))
64
- * Let y = (1-x), z = y/2, s := sqrt(z), and pio2_hi+pio2_lo=pi/2;
65
- * then for x>0.98
66
- * asin(x) = pi/2 - 2*(s+s*z*R(z))
67
- * = pio2_hi - (2*(s+s*z*R(z)) - pio2_lo)
68
- * For x<=0.98, let pio4_hi = pio2_hi/2, then
69
- * f = hi part of s;
70
- * c = sqrt(z) - f = (z-f*f)/(s+f) ...f+c=sqrt(z)
71
- * and
72
- * asin(x) = pi/2 - 2*(s+s*z*R(z))
73
- * = pio4_hi+(pio4-2s)-(2s*z*R(z)-pio2_lo)
74
- * = pio4_hi+(pio4-2f)-(2s*z*R(z)-(pio2_lo+2c))
75
- *
76
- * Special cases:
77
- * if x is NaN, return x itself;
78
- * if |x|>1, return NaN with invalid signal.
79
- *
80
- */
81
-
82
-
83
- #include "fdlibm.h"
84
-
85
- #ifdef __STDC__
86
- static const double
87
- #else
88
- static double
89
- #endif
90
- one = 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */
91
- really_big = 1.000e+300,
92
- pio2_hi = 1.57079632679489655800e+00, /* 0x3FF921FB, 0x54442D18 */
93
- pio2_lo = 6.12323399573676603587e-17, /* 0x3C91A626, 0x33145C07 */
94
- pio4_hi = 7.85398163397448278999e-01, /* 0x3FE921FB, 0x54442D18 */
95
- /* coefficient for R(x^2) */
96
- pS0 = 1.66666666666666657415e-01, /* 0x3FC55555, 0x55555555 */
97
- pS1 = -3.25565818622400915405e-01, /* 0xBFD4D612, 0x03EB6F7D */
98
- pS2 = 2.01212532134862925881e-01, /* 0x3FC9C155, 0x0E884455 */
99
- pS3 = -4.00555345006794114027e-02, /* 0xBFA48228, 0xB5688F3B */
100
- pS4 = 7.91534994289814532176e-04, /* 0x3F49EFE0, 0x7501B288 */
101
- pS5 = 3.47933107596021167570e-05, /* 0x3F023DE1, 0x0DFDF709 */
102
- qS1 = -2.40339491173441421878e+00, /* 0xC0033A27, 0x1C8A2D4B */
103
- qS2 = 2.02094576023350569471e+00, /* 0x40002AE5, 0x9C598AC8 */
104
- qS3 = -6.88283971605453293030e-01, /* 0xBFE6066C, 0x1B8D0159 */
105
- qS4 = 7.70381505559019352791e-02; /* 0x3FB3B8C5, 0xB12E9282 */
106
-
107
- #ifdef __STDC__
108
- double __ieee754_asin(double x)
109
- #else
110
- double __ieee754_asin(x)
111
- double x;
112
- #endif
113
- {
114
- fd_twoints u;
115
- double w,t,p,q,c,r,s;
116
- int hx,ix;
117
- u.d = x;
118
- hx = __HI(u);
119
- x = u.d;
120
- ix = hx&0x7fffffff;
121
- if(ix>= 0x3ff00000) { /* |x|>= 1 */
122
- if(((ix-0x3ff00000)|__LO(u))==0)
123
- /* asin(1)=+-pi/2 with inexact */
124
- return x*pio2_hi+x*pio2_lo;
125
- return (x-x)/(x-x); /* asin(|x|>1) is NaN */
126
- } else if (ix<0x3fe00000) { /* |x|<0.5 */
127
- if(ix<0x3e400000) { /* if |x| < 2**-27 */
128
- if(really_big+x>one) return x;/* return x with inexact if x!=0*/
129
- } else
130
- t = x*x;
131
- p = t*(pS0+t*(pS1+t*(pS2+t*(pS3+t*(pS4+t*pS5)))));
132
- q = one+t*(qS1+t*(qS2+t*(qS3+t*qS4)));
133
- w = p/q;
134
- return x+x*w;
135
- }
136
- /* 1> |x|>= 0.5 */
137
- w = one-fd_fabs(x);
138
- t = w*0.5;
139
- p = t*(pS0+t*(pS1+t*(pS2+t*(pS3+t*(pS4+t*pS5)))));
140
- q = one+t*(qS1+t*(qS2+t*(qS3+t*qS4)));
141
- s = fd_sqrt(t);
142
- if(ix>=0x3FEF3333) { /* if |x| > 0.975 */
143
- w = p/q;
144
- t = pio2_hi-(2.0*(s+s*w)-pio2_lo);
145
- } else {
146
- u.d = s;
147
- __LO(u) = 0;
148
- w = u.d;
149
- c = (t-w*w)/(s+w);
150
- r = p/q;
151
- p = 2.0*s*r-(pio2_lo-2.0*c);
152
- q = pio4_hi-2.0*w;
153
- t = pio4_hi-(p-q);
154
- }
155
- if(hx>0) return t; else return -t;
156
- }