jbarnette-johnson 1.0.0.200806240111 → 1.0.0.200807291507

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (269) hide show
  1. data/MANIFEST +1 -0
  2. data/Rakefile +3 -10
  3. data/bin/johnson +2 -1
  4. data/ext/spidermonkey/context.c +3 -4
  5. data/ext/spidermonkey/context.h +1 -1
  6. data/ext/spidermonkey/conversions.c +39 -33
  7. data/ext/spidermonkey/debugger.c +5 -5
  8. data/ext/spidermonkey/immutable_node.c.erb +11 -11
  9. data/ext/spidermonkey/jroot.h +4 -4
  10. data/ext/spidermonkey/js_land_proxy.c +9 -8
  11. data/ext/spidermonkey/ruby_land_proxy.c +5 -4
  12. data/ext/spidermonkey/runtime.c +1 -1
  13. data/johnson.gemspec +36 -0
  14. data/lib/hoe.rb +0 -7
  15. data/lib/johnson/cli/options.rb +10 -4
  16. data/lib/johnson/spidermonkey/runtime.rb +2 -2
  17. data/lib/johnson/version.rb +4 -2
  18. data/lib/johnson.rb +1 -0
  19. data/test/johnson/runtime_test.rb +11 -0
  20. data/test/johnson/spidermonkey/ruby_land_proxy_test.rb +6 -0
  21. data/vendor/spidermonkey/.cvsignore +9 -0
  22. data/vendor/spidermonkey/Makefile.in +462 -0
  23. data/vendor/spidermonkey/Makefile.ref +364 -0
  24. data/vendor/spidermonkey/README.html +820 -0
  25. data/vendor/spidermonkey/SpiderMonkey.rsp +12 -0
  26. data/vendor/spidermonkey/Y.js +19 -0
  27. data/vendor/spidermonkey/build.mk +43 -0
  28. data/vendor/spidermonkey/config/AIX4.1.mk +65 -0
  29. data/vendor/spidermonkey/config/AIX4.2.mk +64 -0
  30. data/vendor/spidermonkey/config/AIX4.3.mk +65 -0
  31. data/vendor/spidermonkey/config/Darwin.mk +83 -0
  32. data/vendor/spidermonkey/config/Darwin1.3.mk +81 -0
  33. data/vendor/spidermonkey/config/Darwin1.4.mk +41 -0
  34. data/vendor/spidermonkey/config/Darwin5.2.mk +81 -0
  35. data/vendor/spidermonkey/config/Darwin5.3.mk +81 -0
  36. data/vendor/spidermonkey/config/HP-UXB.10.10.mk +77 -0
  37. data/vendor/spidermonkey/config/HP-UXB.10.20.mk +77 -0
  38. data/vendor/spidermonkey/config/HP-UXB.11.00.mk +80 -0
  39. data/vendor/spidermonkey/config/IRIX.mk +87 -0
  40. data/vendor/spidermonkey/config/IRIX5.3.mk +44 -0
  41. data/vendor/spidermonkey/config/IRIX6.1.mk +44 -0
  42. data/vendor/spidermonkey/config/IRIX6.2.mk +44 -0
  43. data/vendor/spidermonkey/config/IRIX6.3.mk +44 -0
  44. data/vendor/spidermonkey/config/IRIX6.5.mk +44 -0
  45. data/vendor/spidermonkey/config/Linux_All.mk +103 -0
  46. data/vendor/spidermonkey/config/Mac_OS10.0.mk +82 -0
  47. data/vendor/spidermonkey/config/OSF1V4.0.mk +72 -0
  48. data/vendor/spidermonkey/config/OSF1V5.0.mk +69 -0
  49. data/vendor/spidermonkey/config/SunOS4.1.4.mk +101 -0
  50. data/vendor/spidermonkey/config/SunOS5.10.mk +50 -0
  51. data/vendor/spidermonkey/config/SunOS5.3.mk +91 -0
  52. data/vendor/spidermonkey/config/SunOS5.4.mk +92 -0
  53. data/vendor/spidermonkey/config/SunOS5.5.1.mk +44 -0
  54. data/vendor/spidermonkey/config/SunOS5.5.mk +87 -0
  55. data/vendor/spidermonkey/config/SunOS5.6.mk +89 -0
  56. data/vendor/spidermonkey/config/SunOS5.7.mk +44 -0
  57. data/vendor/spidermonkey/config/SunOS5.8.mk +44 -0
  58. data/vendor/spidermonkey/config/SunOS5.9.mk +44 -0
  59. data/vendor/spidermonkey/config/WINNT4.0.mk +117 -0
  60. data/vendor/spidermonkey/config/WINNT5.0.mk +117 -0
  61. data/vendor/spidermonkey/config/WINNT5.1.mk +117 -0
  62. data/vendor/spidermonkey/config/WINNT5.2.mk +117 -0
  63. data/vendor/spidermonkey/config/WINNT6.0.mk +117 -0
  64. data/vendor/spidermonkey/config/dgux.mk +64 -0
  65. data/vendor/spidermonkey/config.mk +192 -0
  66. data/vendor/spidermonkey/editline/Makefile.ref +144 -0
  67. data/vendor/spidermonkey/editline/README +83 -0
  68. data/vendor/spidermonkey/editline/editline.3 +175 -0
  69. data/vendor/spidermonkey/editline/editline.c +1369 -0
  70. data/vendor/spidermonkey/editline/editline.h +135 -0
  71. data/vendor/spidermonkey/editline/sysunix.c +182 -0
  72. data/vendor/spidermonkey/editline/unix.h +82 -0
  73. data/vendor/spidermonkey/fdlibm/.cvsignore +7 -0
  74. data/vendor/spidermonkey/fdlibm/Makefile.in +127 -0
  75. data/vendor/spidermonkey/fdlibm/Makefile.ref +192 -0
  76. data/vendor/spidermonkey/fdlibm/e_acos.c +147 -0
  77. data/vendor/spidermonkey/fdlibm/e_acosh.c +105 -0
  78. data/vendor/spidermonkey/fdlibm/e_asin.c +156 -0
  79. data/vendor/spidermonkey/fdlibm/e_atan2.c +165 -0
  80. data/vendor/spidermonkey/fdlibm/e_atanh.c +110 -0
  81. data/vendor/spidermonkey/fdlibm/e_cosh.c +133 -0
  82. data/vendor/spidermonkey/fdlibm/e_exp.c +202 -0
  83. data/vendor/spidermonkey/fdlibm/e_fmod.c +184 -0
  84. data/vendor/spidermonkey/fdlibm/e_gamma.c +71 -0
  85. data/vendor/spidermonkey/fdlibm/e_gamma_r.c +70 -0
  86. data/vendor/spidermonkey/fdlibm/e_hypot.c +173 -0
  87. data/vendor/spidermonkey/fdlibm/e_j0.c +524 -0
  88. data/vendor/spidermonkey/fdlibm/e_j1.c +523 -0
  89. data/vendor/spidermonkey/fdlibm/e_jn.c +315 -0
  90. data/vendor/spidermonkey/fdlibm/e_lgamma.c +71 -0
  91. data/vendor/spidermonkey/fdlibm/e_lgamma_r.c +347 -0
  92. data/vendor/spidermonkey/fdlibm/e_log.c +184 -0
  93. data/vendor/spidermonkey/fdlibm/e_log10.c +134 -0
  94. data/vendor/spidermonkey/fdlibm/e_pow.c +386 -0
  95. data/vendor/spidermonkey/fdlibm/e_rem_pio2.c +222 -0
  96. data/vendor/spidermonkey/fdlibm/e_remainder.c +120 -0
  97. data/vendor/spidermonkey/fdlibm/e_scalb.c +89 -0
  98. data/vendor/spidermonkey/fdlibm/e_sinh.c +122 -0
  99. data/vendor/spidermonkey/fdlibm/e_sqrt.c +497 -0
  100. data/vendor/spidermonkey/fdlibm/fdlibm.h +273 -0
  101. data/vendor/spidermonkey/fdlibm/fdlibm.mak +1453 -0
  102. data/vendor/spidermonkey/fdlibm/fdlibm.mdp +0 -0
  103. data/vendor/spidermonkey/fdlibm/k_cos.c +135 -0
  104. data/vendor/spidermonkey/fdlibm/k_rem_pio2.c +354 -0
  105. data/vendor/spidermonkey/fdlibm/k_sin.c +114 -0
  106. data/vendor/spidermonkey/fdlibm/k_standard.c +785 -0
  107. data/vendor/spidermonkey/fdlibm/k_tan.c +170 -0
  108. data/vendor/spidermonkey/fdlibm/s_asinh.c +101 -0
  109. data/vendor/spidermonkey/fdlibm/s_atan.c +175 -0
  110. data/vendor/spidermonkey/fdlibm/s_cbrt.c +133 -0
  111. data/vendor/spidermonkey/fdlibm/s_ceil.c +120 -0
  112. data/vendor/spidermonkey/fdlibm/s_copysign.c +72 -0
  113. data/vendor/spidermonkey/fdlibm/s_cos.c +118 -0
  114. data/vendor/spidermonkey/fdlibm/s_erf.c +356 -0
  115. data/vendor/spidermonkey/fdlibm/s_expm1.c +267 -0
  116. data/vendor/spidermonkey/fdlibm/s_fabs.c +70 -0
  117. data/vendor/spidermonkey/fdlibm/s_finite.c +71 -0
  118. data/vendor/spidermonkey/fdlibm/s_floor.c +121 -0
  119. data/vendor/spidermonkey/fdlibm/s_frexp.c +99 -0
  120. data/vendor/spidermonkey/fdlibm/s_ilogb.c +85 -0
  121. data/vendor/spidermonkey/fdlibm/s_isnan.c +74 -0
  122. data/vendor/spidermonkey/fdlibm/s_ldexp.c +66 -0
  123. data/vendor/spidermonkey/fdlibm/s_lib_version.c +73 -0
  124. data/vendor/spidermonkey/fdlibm/s_log1p.c +211 -0
  125. data/vendor/spidermonkey/fdlibm/s_logb.c +79 -0
  126. data/vendor/spidermonkey/fdlibm/s_matherr.c +64 -0
  127. data/vendor/spidermonkey/fdlibm/s_modf.c +132 -0
  128. data/vendor/spidermonkey/fdlibm/s_nextafter.c +124 -0
  129. data/vendor/spidermonkey/fdlibm/s_rint.c +131 -0
  130. data/vendor/spidermonkey/fdlibm/s_scalbn.c +107 -0
  131. data/vendor/spidermonkey/fdlibm/s_signgam.c +40 -0
  132. data/vendor/spidermonkey/fdlibm/s_significand.c +68 -0
  133. data/vendor/spidermonkey/fdlibm/s_sin.c +118 -0
  134. data/vendor/spidermonkey/fdlibm/s_tan.c +112 -0
  135. data/vendor/spidermonkey/fdlibm/s_tanh.c +122 -0
  136. data/vendor/spidermonkey/fdlibm/w_acos.c +78 -0
  137. data/vendor/spidermonkey/fdlibm/w_acosh.c +78 -0
  138. data/vendor/spidermonkey/fdlibm/w_asin.c +80 -0
  139. data/vendor/spidermonkey/fdlibm/w_atan2.c +79 -0
  140. data/vendor/spidermonkey/fdlibm/w_atanh.c +81 -0
  141. data/vendor/spidermonkey/fdlibm/w_cosh.c +77 -0
  142. data/vendor/spidermonkey/fdlibm/w_exp.c +88 -0
  143. data/vendor/spidermonkey/fdlibm/w_fmod.c +78 -0
  144. data/vendor/spidermonkey/fdlibm/w_gamma.c +85 -0
  145. data/vendor/spidermonkey/fdlibm/w_gamma_r.c +81 -0
  146. data/vendor/spidermonkey/fdlibm/w_hypot.c +78 -0
  147. data/vendor/spidermonkey/fdlibm/w_j0.c +105 -0
  148. data/vendor/spidermonkey/fdlibm/w_j1.c +106 -0
  149. data/vendor/spidermonkey/fdlibm/w_jn.c +128 -0
  150. data/vendor/spidermonkey/fdlibm/w_lgamma.c +85 -0
  151. data/vendor/spidermonkey/fdlibm/w_lgamma_r.c +81 -0
  152. data/vendor/spidermonkey/fdlibm/w_log.c +78 -0
  153. data/vendor/spidermonkey/fdlibm/w_log10.c +81 -0
  154. data/vendor/spidermonkey/fdlibm/w_pow.c +99 -0
  155. data/vendor/spidermonkey/fdlibm/w_remainder.c +77 -0
  156. data/vendor/spidermonkey/fdlibm/w_scalb.c +95 -0
  157. data/vendor/spidermonkey/fdlibm/w_sinh.c +77 -0
  158. data/vendor/spidermonkey/fdlibm/w_sqrt.c +77 -0
  159. data/vendor/spidermonkey/javascript-trace.d +73 -0
  160. data/vendor/spidermonkey/js.c +3951 -0
  161. data/vendor/spidermonkey/js.mak +4438 -0
  162. data/vendor/spidermonkey/js.mdp +0 -0
  163. data/vendor/spidermonkey/js.msg +307 -0
  164. data/vendor/spidermonkey/js.pkg +2 -0
  165. data/vendor/spidermonkey/js3240.rc +79 -0
  166. data/vendor/spidermonkey/jsOS240.def +654 -0
  167. data/vendor/spidermonkey/jsapi.c +5836 -0
  168. data/vendor/spidermonkey/jsapi.h +2624 -0
  169. data/vendor/spidermonkey/jsarena.c +450 -0
  170. data/vendor/spidermonkey/jsarena.h +318 -0
  171. data/vendor/spidermonkey/jsarray.c +2988 -0
  172. data/vendor/spidermonkey/jsarray.h +124 -0
  173. data/vendor/spidermonkey/jsatom.c +1045 -0
  174. data/vendor/spidermonkey/jsatom.h +442 -0
  175. data/vendor/spidermonkey/jsbit.h +253 -0
  176. data/vendor/spidermonkey/jsbool.c +176 -0
  177. data/vendor/spidermonkey/jsbool.h +73 -0
  178. data/vendor/spidermonkey/jsclist.h +139 -0
  179. data/vendor/spidermonkey/jscntxt.c +1348 -0
  180. data/vendor/spidermonkey/jscntxt.h +1120 -0
  181. data/vendor/spidermonkey/jscompat.h +57 -0
  182. data/vendor/spidermonkey/jsconfig.h +248 -0
  183. data/vendor/spidermonkey/jsconfig.mk +181 -0
  184. data/vendor/spidermonkey/jscpucfg.c +383 -0
  185. data/vendor/spidermonkey/jscpucfg.h +212 -0
  186. data/vendor/spidermonkey/jsdate.c +2398 -0
  187. data/vendor/spidermonkey/jsdate.h +124 -0
  188. data/vendor/spidermonkey/jsdbgapi.c +1799 -0
  189. data/vendor/spidermonkey/jsdbgapi.h +464 -0
  190. data/vendor/spidermonkey/jsdhash.c +868 -0
  191. data/vendor/spidermonkey/jsdhash.h +592 -0
  192. data/vendor/spidermonkey/jsdtoa.c +3167 -0
  193. data/vendor/spidermonkey/jsdtoa.h +130 -0
  194. data/vendor/spidermonkey/jsdtracef.c +317 -0
  195. data/vendor/spidermonkey/jsdtracef.h +77 -0
  196. data/vendor/spidermonkey/jsemit.c +6909 -0
  197. data/vendor/spidermonkey/jsemit.h +741 -0
  198. data/vendor/spidermonkey/jsexn.c +1371 -0
  199. data/vendor/spidermonkey/jsexn.h +96 -0
  200. data/vendor/spidermonkey/jsfile.c +2736 -0
  201. data/vendor/spidermonkey/jsfile.h +56 -0
  202. data/vendor/spidermonkey/jsfile.msg +90 -0
  203. data/vendor/spidermonkey/jsfun.c +2634 -0
  204. data/vendor/spidermonkey/jsfun.h +254 -0
  205. data/vendor/spidermonkey/jsgc.c +3554 -0
  206. data/vendor/spidermonkey/jsgc.h +403 -0
  207. data/vendor/spidermonkey/jshash.c +476 -0
  208. data/vendor/spidermonkey/jshash.h +151 -0
  209. data/vendor/spidermonkey/jsify.pl +485 -0
  210. data/vendor/spidermonkey/jsinterp.c +6981 -0
  211. data/vendor/spidermonkey/jsinterp.h +521 -0
  212. data/vendor/spidermonkey/jsinvoke.c +43 -0
  213. data/vendor/spidermonkey/jsiter.c +1067 -0
  214. data/vendor/spidermonkey/jsiter.h +122 -0
  215. data/vendor/spidermonkey/jskeyword.tbl +124 -0
  216. data/vendor/spidermonkey/jskwgen.c +460 -0
  217. data/vendor/spidermonkey/jslibmath.h +266 -0
  218. data/vendor/spidermonkey/jslock.c +1309 -0
  219. data/vendor/spidermonkey/jslock.h +313 -0
  220. data/vendor/spidermonkey/jslocko.asm +60 -0
  221. data/vendor/spidermonkey/jslog2.c +94 -0
  222. data/vendor/spidermonkey/jslong.c +264 -0
  223. data/vendor/spidermonkey/jslong.h +412 -0
  224. data/vendor/spidermonkey/jsmath.c +568 -0
  225. data/vendor/spidermonkey/jsmath.h +57 -0
  226. data/vendor/spidermonkey/jsnum.c +1228 -0
  227. data/vendor/spidermonkey/jsnum.h +283 -0
  228. data/vendor/spidermonkey/jsobj.c +5266 -0
  229. data/vendor/spidermonkey/jsobj.h +709 -0
  230. data/vendor/spidermonkey/jsopcode.c +5245 -0
  231. data/vendor/spidermonkey/jsopcode.h +394 -0
  232. data/vendor/spidermonkey/jsopcode.tbl +523 -0
  233. data/vendor/spidermonkey/jsotypes.h +202 -0
  234. data/vendor/spidermonkey/jsparse.c +6680 -0
  235. data/vendor/spidermonkey/jsparse.h +511 -0
  236. data/vendor/spidermonkey/jsprf.c +1262 -0
  237. data/vendor/spidermonkey/jsprf.h +150 -0
  238. data/vendor/spidermonkey/jsproto.tbl +128 -0
  239. data/vendor/spidermonkey/jsprvtd.h +267 -0
  240. data/vendor/spidermonkey/jspubtd.h +744 -0
  241. data/vendor/spidermonkey/jsregexp.c +4352 -0
  242. data/vendor/spidermonkey/jsregexp.h +183 -0
  243. data/vendor/spidermonkey/jsreops.tbl +145 -0
  244. data/vendor/spidermonkey/jsscan.c +2003 -0
  245. data/vendor/spidermonkey/jsscan.h +387 -0
  246. data/vendor/spidermonkey/jsscope.c +1948 -0
  247. data/vendor/spidermonkey/jsscope.h +418 -0
  248. data/vendor/spidermonkey/jsscript.c +1832 -0
  249. data/vendor/spidermonkey/jsscript.h +287 -0
  250. data/vendor/spidermonkey/jsshell.msg +50 -0
  251. data/vendor/spidermonkey/jsstddef.h +83 -0
  252. data/vendor/spidermonkey/jsstr.c +5004 -0
  253. data/vendor/spidermonkey/jsstr.h +641 -0
  254. data/vendor/spidermonkey/jstypes.h +475 -0
  255. data/vendor/spidermonkey/jsutil.c +345 -0
  256. data/vendor/spidermonkey/jsutil.h +157 -0
  257. data/vendor/spidermonkey/jsxdrapi.c +800 -0
  258. data/vendor/spidermonkey/jsxdrapi.h +218 -0
  259. data/vendor/spidermonkey/jsxml.c +8471 -0
  260. data/vendor/spidermonkey/jsxml.h +349 -0
  261. data/vendor/spidermonkey/lock_SunOS.s +119 -0
  262. data/vendor/spidermonkey/perfect.js +39 -0
  263. data/vendor/spidermonkey/plify_jsdhash.sed +36 -0
  264. data/vendor/spidermonkey/prmjtime.c +846 -0
  265. data/vendor/spidermonkey/prmjtime.h +103 -0
  266. data/vendor/spidermonkey/resource.h +15 -0
  267. data/vendor/spidermonkey/rules.mk +197 -0
  268. data/vendor/spidermonkey/win32.order +384 -0
  269. metadata +4 -3
@@ -0,0 +1,170 @@
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
+ /* @(#)k_tan.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
+ /* __kernel_tan( x, y, k )
53
+ * kernel tan function on [-pi/4, pi/4], pi/4 ~ 0.7854
54
+ * Input x is assumed to be bounded by ~pi/4 in magnitude.
55
+ * Input y is the tail of x.
56
+ * Input k indicates whether tan (if k=1) or
57
+ * -1/tan (if k= -1) is returned.
58
+ *
59
+ * Algorithm
60
+ * 1. Since tan(-x) = -tan(x), we need only to consider positive x.
61
+ * 2. if x < 2^-28 (hx<0x3e300000 0), return x with inexact if x!=0.
62
+ * 3. tan(x) is approximated by a odd polynomial of degree 27 on
63
+ * [0,0.67434]
64
+ * 3 27
65
+ * tan(x) ~ x + T1*x + ... + T13*x
66
+ * where
67
+ *
68
+ * |tan(x) 2 4 26 | -59.2
69
+ * |----- - (1+T1*x +T2*x +.... +T13*x )| <= 2
70
+ * | x |
71
+ *
72
+ * Note: tan(x+y) = tan(x) + tan'(x)*y
73
+ * ~ tan(x) + (1+x*x)*y
74
+ * Therefore, for better accuracy in computing tan(x+y), let
75
+ * 3 2 2 2 2
76
+ * r = x *(T2+x *(T3+x *(...+x *(T12+x *T13))))
77
+ * then
78
+ * 3 2
79
+ * tan(x+y) = x + (T1*x + (x *(r+y)+y))
80
+ *
81
+ * 4. For x in [0.67434,pi/4], let y = pi/4 - x, then
82
+ * tan(x) = tan(pi/4-y) = (1-tan(y))/(1+tan(y))
83
+ * = 1 - 2*(tan(y) - (tan(y)^2)/(1+tan(y)))
84
+ */
85
+
86
+ #include "fdlibm.h"
87
+ #ifdef __STDC__
88
+ static const double
89
+ #else
90
+ static double
91
+ #endif
92
+ one = 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */
93
+ pio4 = 7.85398163397448278999e-01, /* 0x3FE921FB, 0x54442D18 */
94
+ pio4lo= 3.06161699786838301793e-17, /* 0x3C81A626, 0x33145C07 */
95
+ T[] = {
96
+ 3.33333333333334091986e-01, /* 0x3FD55555, 0x55555563 */
97
+ 1.33333333333201242699e-01, /* 0x3FC11111, 0x1110FE7A */
98
+ 5.39682539762260521377e-02, /* 0x3FABA1BA, 0x1BB341FE */
99
+ 2.18694882948595424599e-02, /* 0x3F9664F4, 0x8406D637 */
100
+ 8.86323982359930005737e-03, /* 0x3F8226E3, 0xE96E8493 */
101
+ 3.59207910759131235356e-03, /* 0x3F6D6D22, 0xC9560328 */
102
+ 1.45620945432529025516e-03, /* 0x3F57DBC8, 0xFEE08315 */
103
+ 5.88041240820264096874e-04, /* 0x3F4344D8, 0xF2F26501 */
104
+ 2.46463134818469906812e-04, /* 0x3F3026F7, 0x1A8D1068 */
105
+ 7.81794442939557092300e-05, /* 0x3F147E88, 0xA03792A6 */
106
+ 7.14072491382608190305e-05, /* 0x3F12B80F, 0x32F0A7E9 */
107
+ -1.85586374855275456654e-05, /* 0xBEF375CB, 0xDB605373 */
108
+ 2.59073051863633712884e-05, /* 0x3EFB2A70, 0x74BF7AD4 */
109
+ };
110
+
111
+ #ifdef __STDC__
112
+ double __kernel_tan(double x, double y, int iy)
113
+ #else
114
+ double __kernel_tan(x, y, iy)
115
+ double x,y; int iy;
116
+ #endif
117
+ {
118
+ fd_twoints u;
119
+ double z,r,v,w,s;
120
+ int ix,hx;
121
+ u.d = x;
122
+ hx = __HI(u); /* high word of x */
123
+ ix = hx&0x7fffffff; /* high word of |x| */
124
+ if(ix<0x3e300000) /* x < 2**-28 */
125
+ {if((int)x==0) { /* generate inexact */
126
+ u.d =x;
127
+ if(((ix|__LO(u))|(iy+1))==0) return one/fd_fabs(x);
128
+ else return (iy==1)? x: -one/x;
129
+ }
130
+ }
131
+ if(ix>=0x3FE59428) { /* |x|>=0.6744 */
132
+ if(hx<0) {x = -x; y = -y;}
133
+ z = pio4-x;
134
+ w = pio4lo-y;
135
+ x = z+w; y = 0.0;
136
+ }
137
+ z = x*x;
138
+ w = z*z;
139
+ /* Break x^5*(T[1]+x^2*T[2]+...) into
140
+ * x^5(T[1]+x^4*T[3]+...+x^20*T[11]) +
141
+ * x^5(x^2*(T[2]+x^4*T[4]+...+x^22*[T12]))
142
+ */
143
+ r = T[1]+w*(T[3]+w*(T[5]+w*(T[7]+w*(T[9]+w*T[11]))));
144
+ v = z*(T[2]+w*(T[4]+w*(T[6]+w*(T[8]+w*(T[10]+w*T[12])))));
145
+ s = z*x;
146
+ r = y + z*(s*(r+v)+y);
147
+ r += T[0]*s;
148
+ w = x+r;
149
+ if(ix>=0x3FE59428) {
150
+ v = (double)iy;
151
+ return (double)(1-((hx>>30)&2))*(v-2.0*(x-(w*w/(w+v)-r)));
152
+ }
153
+ if(iy==1) return w;
154
+ else { /* if allow error up to 2 ulp,
155
+ simply return -1.0/(x+r) here */
156
+ /* compute -1.0/(x+r) accurately */
157
+ double a,t;
158
+ z = w;
159
+ u.d = z;
160
+ __LO(u) = 0;
161
+ z = u.d;
162
+ v = r-(z - x); /* z+v = r+x */
163
+ t = a = -1.0/w; /* a = -1.0/w */
164
+ u.d = t;
165
+ __LO(u) = 0;
166
+ t = u.d;
167
+ s = 1.0+t*z;
168
+ return t+a*(s+t*v);
169
+ }
170
+ }
@@ -0,0 +1,101 @@
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
+ /* @(#)s_asinh.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
+ /* asinh(x)
53
+ * Method :
54
+ * Based on
55
+ * asinh(x) = sign(x) * log [ |x| + sqrt(x*x+1) ]
56
+ * we have
57
+ * asinh(x) := x if 1+x*x=1,
58
+ * := sign(x)*(log(x)+ln2)) for large |x|, else
59
+ * := sign(x)*log(2|x|+1/(|x|+sqrt(x*x+1))) if|x|>2, else
60
+ * := sign(x)*log1p(|x| + x^2/(1 + sqrt(1+x^2)))
61
+ */
62
+
63
+ #include "fdlibm.h"
64
+
65
+ #ifdef __STDC__
66
+ static const double
67
+ #else
68
+ static double
69
+ #endif
70
+ one = 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */
71
+ ln2 = 6.93147180559945286227e-01, /* 0x3FE62E42, 0xFEFA39EF */
72
+ really_big= 1.00000000000000000000e+300;
73
+
74
+ #ifdef __STDC__
75
+ double fd_asinh(double x)
76
+ #else
77
+ double fd_asinh(x)
78
+ double x;
79
+ #endif
80
+ {
81
+ fd_twoints u;
82
+ double t,w;
83
+ int hx,ix;
84
+ u.d = x;
85
+ hx = __HI(u);
86
+ ix = hx&0x7fffffff;
87
+ if(ix>=0x7ff00000) return x+x; /* x is inf or NaN */
88
+ if(ix< 0x3e300000) { /* |x|<2**-28 */
89
+ if(really_big+x>one) return x; /* return x inexact except 0 */
90
+ }
91
+ if(ix>0x41b00000) { /* |x| > 2**28 */
92
+ w = __ieee754_log(fd_fabs(x))+ln2;
93
+ } else if (ix>0x40000000) { /* 2**28 > |x| > 2.0 */
94
+ t = fd_fabs(x);
95
+ w = __ieee754_log(2.0*t+one/(fd_sqrt(x*x+one)+t));
96
+ } else { /* 2.0 > |x| > 2**-28 */
97
+ t = x*x;
98
+ w =fd_log1p(fd_fabs(x)+t/(one+fd_sqrt(one+t)));
99
+ }
100
+ if(hx>0) return w; else return -w;
101
+ }
@@ -0,0 +1,175 @@
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
+ /* @(#)s_atan.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
+ /* atan(x)
54
+ * Method
55
+ * 1. Reduce x to positive by atan(x) = -atan(-x).
56
+ * 2. According to the integer k=4t+0.25 chopped, t=x, the argument
57
+ * is further reduced to one of the following intervals and the
58
+ * arctangent of t is evaluated by the corresponding formula:
59
+ *
60
+ * [0,7/16] atan(x) = t-t^3*(a1+t^2*(a2+...(a10+t^2*a11)...)
61
+ * [7/16,11/16] atan(x) = atan(1/2) + atan( (t-0.5)/(1+t/2) )
62
+ * [11/16.19/16] atan(x) = atan( 1 ) + atan( (t-1)/(1+t) )
63
+ * [19/16,39/16] atan(x) = atan(3/2) + atan( (t-1.5)/(1+1.5t) )
64
+ * [39/16,INF] atan(x) = atan(INF) + atan( -1/t )
65
+ *
66
+ * Constants:
67
+ * The hexadecimal values are the intended ones for the following
68
+ * constants. The decimal values may be used, provided that the
69
+ * compiler will convert from decimal to binary accurately enough
70
+ * to produce the hexadecimal values shown.
71
+ */
72
+
73
+ #include "fdlibm.h"
74
+
75
+ #ifdef __STDC__
76
+ static const double atanhi[] = {
77
+ #else
78
+ static double atanhi[] = {
79
+ #endif
80
+ 4.63647609000806093515e-01, /* atan(0.5)hi 0x3FDDAC67, 0x0561BB4F */
81
+ 7.85398163397448278999e-01, /* atan(1.0)hi 0x3FE921FB, 0x54442D18 */
82
+ 9.82793723247329054082e-01, /* atan(1.5)hi 0x3FEF730B, 0xD281F69B */
83
+ 1.57079632679489655800e+00, /* atan(inf)hi 0x3FF921FB, 0x54442D18 */
84
+ };
85
+
86
+ #ifdef __STDC__
87
+ static const double atanlo[] = {
88
+ #else
89
+ static double atanlo[] = {
90
+ #endif
91
+ 2.26987774529616870924e-17, /* atan(0.5)lo 0x3C7A2B7F, 0x222F65E2 */
92
+ 3.06161699786838301793e-17, /* atan(1.0)lo 0x3C81A626, 0x33145C07 */
93
+ 1.39033110312309984516e-17, /* atan(1.5)lo 0x3C700788, 0x7AF0CBBD */
94
+ 6.12323399573676603587e-17, /* atan(inf)lo 0x3C91A626, 0x33145C07 */
95
+ };
96
+
97
+ #ifdef __STDC__
98
+ static const double aT[] = {
99
+ #else
100
+ static double aT[] = {
101
+ #endif
102
+ 3.33333333333329318027e-01, /* 0x3FD55555, 0x5555550D */
103
+ -1.99999999998764832476e-01, /* 0xBFC99999, 0x9998EBC4 */
104
+ 1.42857142725034663711e-01, /* 0x3FC24924, 0x920083FF */
105
+ -1.11111104054623557880e-01, /* 0xBFBC71C6, 0xFE231671 */
106
+ 9.09088713343650656196e-02, /* 0x3FB745CD, 0xC54C206E */
107
+ -7.69187620504482999495e-02, /* 0xBFB3B0F2, 0xAF749A6D */
108
+ 6.66107313738753120669e-02, /* 0x3FB10D66, 0xA0D03D51 */
109
+ -5.83357013379057348645e-02, /* 0xBFADDE2D, 0x52DEFD9A */
110
+ 4.97687799461593236017e-02, /* 0x3FA97B4B, 0x24760DEB */
111
+ -3.65315727442169155270e-02, /* 0xBFA2B444, 0x2C6A6C2F */
112
+ 1.62858201153657823623e-02, /* 0x3F90AD3A, 0xE322DA11 */
113
+ };
114
+
115
+ #ifdef __STDC__
116
+ static const double
117
+ #else
118
+ static double
119
+ #endif
120
+ one = 1.0,
121
+ really_big = 1.0e300;
122
+
123
+ #ifdef __STDC__
124
+ double fd_atan(double x)
125
+ #else
126
+ double fd_atan(x)
127
+ double x;
128
+ #endif
129
+ {
130
+ fd_twoints u;
131
+ double w,s1,s2,z;
132
+ int ix,hx,id;
133
+
134
+ u.d = x;
135
+ hx = __HI(u);
136
+ ix = hx&0x7fffffff;
137
+ if(ix>=0x44100000) { /* if |x| >= 2^66 */
138
+ u.d = x;
139
+ if(ix>0x7ff00000||
140
+ (ix==0x7ff00000&&(__LO(u)!=0)))
141
+ return x+x; /* NaN */
142
+ if(hx>0) return atanhi[3]+atanlo[3];
143
+ else return -atanhi[3]-atanlo[3];
144
+ } if (ix < 0x3fdc0000) { /* |x| < 0.4375 */
145
+ if (ix < 0x3e200000) { /* |x| < 2^-29 */
146
+ if(really_big+x>one) return x; /* raise inexact */
147
+ }
148
+ id = -1;
149
+ } else {
150
+ x = fd_fabs(x);
151
+ if (ix < 0x3ff30000) { /* |x| < 1.1875 */
152
+ if (ix < 0x3fe60000) { /* 7/16 <=|x|<11/16 */
153
+ id = 0; x = (2.0*x-one)/(2.0+x);
154
+ } else { /* 11/16<=|x|< 19/16 */
155
+ id = 1; x = (x-one)/(x+one);
156
+ }
157
+ } else {
158
+ if (ix < 0x40038000) { /* |x| < 2.4375 */
159
+ id = 2; x = (x-1.5)/(one+1.5*x);
160
+ } else { /* 2.4375 <= |x| < 2^66 */
161
+ id = 3; x = -1.0/x;
162
+ }
163
+ }}
164
+ /* end of argument reduction */
165
+ z = x*x;
166
+ w = z*z;
167
+ /* break sum from i=0 to 10 aT[i]z**(i+1) into odd and even poly */
168
+ s1 = z*(aT[0]+w*(aT[2]+w*(aT[4]+w*(aT[6]+w*(aT[8]+w*aT[10])))));
169
+ s2 = w*(aT[1]+w*(aT[3]+w*(aT[5]+w*(aT[7]+w*aT[9]))));
170
+ if (id<0) return x - x*(s1+s2);
171
+ else {
172
+ z = atanhi[id] - ((x*(s1+s2) - atanlo[id]) - x);
173
+ return (hx<0)? -z:z;
174
+ }
175
+ }
@@ -0,0 +1,133 @@
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
+ /* @(#)s_cbrt.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
+ #include "fdlibm.h"
54
+
55
+ /* cbrt(x)
56
+ * Return cube root of x
57
+ */
58
+ #ifdef __STDC__
59
+ static const unsigned
60
+ #else
61
+ static unsigned
62
+ #endif
63
+ B1 = 715094163, /* B1 = (682-0.03306235651)*2**20 */
64
+ B2 = 696219795; /* B2 = (664-0.03306235651)*2**20 */
65
+
66
+ #ifdef __STDC__
67
+ static const double
68
+ #else
69
+ static double
70
+ #endif
71
+ C = 5.42857142857142815906e-01, /* 19/35 = 0x3FE15F15, 0xF15F15F1 */
72
+ D = -7.05306122448979611050e-01, /* -864/1225 = 0xBFE691DE, 0x2532C834 */
73
+ E = 1.41428571428571436819e+00, /* 99/70 = 0x3FF6A0EA, 0x0EA0EA0F */
74
+ F = 1.60714285714285720630e+00, /* 45/28 = 0x3FF9B6DB, 0x6DB6DB6E */
75
+ G = 3.57142857142857150787e-01; /* 5/14 = 0x3FD6DB6D, 0xB6DB6DB7 */
76
+
77
+ #ifdef __STDC__
78
+ double fd_cbrt(double x)
79
+ #else
80
+ double fd_cbrt(x)
81
+ double x;
82
+ #endif
83
+ {
84
+ fd_twoints u;
85
+ int hx;
86
+ double r,s,t=0.0,w;
87
+ unsigned sign;
88
+
89
+ u.d = x;
90
+ hx = __HI(u); /* high word of x */
91
+ sign=hx&0x80000000; /* sign= sign(x) */
92
+ hx ^=sign;
93
+ if(hx>=0x7ff00000) return(x+x); /* cbrt(NaN,INF) is itself */
94
+ if((hx|__LO(u))==0) {
95
+ x = u.d;
96
+ return(x); /* cbrt(0) is itself */
97
+ }
98
+ u.d = x;
99
+ __HI(u) = hx; /* x <- |x| */
100
+ x = u.d;
101
+ /* rough cbrt to 5 bits */
102
+ if(hx<0x00100000) /* subnormal number */
103
+ {u.d = t; __HI(u)=0x43500000; t=u.d; /* set t= 2**54 */
104
+ t*=x; __HI(u)=__HI(u)/3+B2;
105
+ }
106
+ else {
107
+ u.d = t; __HI(u)=hx/3+B1; t = u.d;
108
+ }
109
+
110
+
111
+ /* new cbrt to 23 bits, may be implemented in single precision */
112
+ r=t*t/x;
113
+ s=C+r*t;
114
+ t*=G+F/(s+E+D/s);
115
+
116
+ /* chopped to 20 bits and make it larger than cbrt(x) */
117
+ u.d = t;
118
+ __LO(u)=0; __HI(u)+=0x00000001;
119
+ t = u.d;
120
+
121
+ /* one step newton iteration to 53 bits with error less than 0.667 ulps */
122
+ s=t*t; /* t*t is exact */
123
+ r=x/s;
124
+ w=t+t;
125
+ r=(r-t)/(w+r); /* r-s is exact */
126
+ t=t+t*r;
127
+
128
+ /* retore the sign bit */
129
+ u.d = t;
130
+ __HI(u) |= sign;
131
+ t = u.d;
132
+ return(t);
133
+ }