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.
- data/MANIFEST +1 -0
- data/Rakefile +3 -10
- data/bin/johnson +2 -1
- data/ext/spidermonkey/context.c +3 -4
- data/ext/spidermonkey/context.h +1 -1
- data/ext/spidermonkey/conversions.c +39 -33
- data/ext/spidermonkey/debugger.c +5 -5
- data/ext/spidermonkey/immutable_node.c.erb +11 -11
- data/ext/spidermonkey/jroot.h +4 -4
- data/ext/spidermonkey/js_land_proxy.c +9 -8
- data/ext/spidermonkey/ruby_land_proxy.c +5 -4
- data/ext/spidermonkey/runtime.c +1 -1
- data/johnson.gemspec +36 -0
- data/lib/hoe.rb +0 -7
- data/lib/johnson/cli/options.rb +10 -4
- data/lib/johnson/spidermonkey/runtime.rb +2 -2
- data/lib/johnson/version.rb +4 -2
- data/lib/johnson.rb +1 -0
- data/test/johnson/runtime_test.rb +11 -0
- data/test/johnson/spidermonkey/ruby_land_proxy_test.rb +6 -0
- data/vendor/spidermonkey/.cvsignore +9 -0
- data/vendor/spidermonkey/Makefile.in +462 -0
- data/vendor/spidermonkey/Makefile.ref +364 -0
- data/vendor/spidermonkey/README.html +820 -0
- data/vendor/spidermonkey/SpiderMonkey.rsp +12 -0
- data/vendor/spidermonkey/Y.js +19 -0
- data/vendor/spidermonkey/build.mk +43 -0
- data/vendor/spidermonkey/config/AIX4.1.mk +65 -0
- data/vendor/spidermonkey/config/AIX4.2.mk +64 -0
- data/vendor/spidermonkey/config/AIX4.3.mk +65 -0
- data/vendor/spidermonkey/config/Darwin.mk +83 -0
- data/vendor/spidermonkey/config/Darwin1.3.mk +81 -0
- data/vendor/spidermonkey/config/Darwin1.4.mk +41 -0
- data/vendor/spidermonkey/config/Darwin5.2.mk +81 -0
- data/vendor/spidermonkey/config/Darwin5.3.mk +81 -0
- data/vendor/spidermonkey/config/HP-UXB.10.10.mk +77 -0
- data/vendor/spidermonkey/config/HP-UXB.10.20.mk +77 -0
- data/vendor/spidermonkey/config/HP-UXB.11.00.mk +80 -0
- data/vendor/spidermonkey/config/IRIX.mk +87 -0
- data/vendor/spidermonkey/config/IRIX5.3.mk +44 -0
- data/vendor/spidermonkey/config/IRIX6.1.mk +44 -0
- data/vendor/spidermonkey/config/IRIX6.2.mk +44 -0
- data/vendor/spidermonkey/config/IRIX6.3.mk +44 -0
- data/vendor/spidermonkey/config/IRIX6.5.mk +44 -0
- data/vendor/spidermonkey/config/Linux_All.mk +103 -0
- data/vendor/spidermonkey/config/Mac_OS10.0.mk +82 -0
- data/vendor/spidermonkey/config/OSF1V4.0.mk +72 -0
- data/vendor/spidermonkey/config/OSF1V5.0.mk +69 -0
- data/vendor/spidermonkey/config/SunOS4.1.4.mk +101 -0
- data/vendor/spidermonkey/config/SunOS5.10.mk +50 -0
- data/vendor/spidermonkey/config/SunOS5.3.mk +91 -0
- data/vendor/spidermonkey/config/SunOS5.4.mk +92 -0
- data/vendor/spidermonkey/config/SunOS5.5.1.mk +44 -0
- data/vendor/spidermonkey/config/SunOS5.5.mk +87 -0
- data/vendor/spidermonkey/config/SunOS5.6.mk +89 -0
- data/vendor/spidermonkey/config/SunOS5.7.mk +44 -0
- data/vendor/spidermonkey/config/SunOS5.8.mk +44 -0
- data/vendor/spidermonkey/config/SunOS5.9.mk +44 -0
- data/vendor/spidermonkey/config/WINNT4.0.mk +117 -0
- data/vendor/spidermonkey/config/WINNT5.0.mk +117 -0
- data/vendor/spidermonkey/config/WINNT5.1.mk +117 -0
- data/vendor/spidermonkey/config/WINNT5.2.mk +117 -0
- data/vendor/spidermonkey/config/WINNT6.0.mk +117 -0
- data/vendor/spidermonkey/config/dgux.mk +64 -0
- data/vendor/spidermonkey/config.mk +192 -0
- data/vendor/spidermonkey/editline/Makefile.ref +144 -0
- data/vendor/spidermonkey/editline/README +83 -0
- data/vendor/spidermonkey/editline/editline.3 +175 -0
- data/vendor/spidermonkey/editline/editline.c +1369 -0
- data/vendor/spidermonkey/editline/editline.h +135 -0
- data/vendor/spidermonkey/editline/sysunix.c +182 -0
- data/vendor/spidermonkey/editline/unix.h +82 -0
- data/vendor/spidermonkey/fdlibm/.cvsignore +7 -0
- data/vendor/spidermonkey/fdlibm/Makefile.in +127 -0
- data/vendor/spidermonkey/fdlibm/Makefile.ref +192 -0
- data/vendor/spidermonkey/fdlibm/e_acos.c +147 -0
- data/vendor/spidermonkey/fdlibm/e_acosh.c +105 -0
- data/vendor/spidermonkey/fdlibm/e_asin.c +156 -0
- data/vendor/spidermonkey/fdlibm/e_atan2.c +165 -0
- data/vendor/spidermonkey/fdlibm/e_atanh.c +110 -0
- data/vendor/spidermonkey/fdlibm/e_cosh.c +133 -0
- data/vendor/spidermonkey/fdlibm/e_exp.c +202 -0
- data/vendor/spidermonkey/fdlibm/e_fmod.c +184 -0
- data/vendor/spidermonkey/fdlibm/e_gamma.c +71 -0
- data/vendor/spidermonkey/fdlibm/e_gamma_r.c +70 -0
- data/vendor/spidermonkey/fdlibm/e_hypot.c +173 -0
- data/vendor/spidermonkey/fdlibm/e_j0.c +524 -0
- data/vendor/spidermonkey/fdlibm/e_j1.c +523 -0
- data/vendor/spidermonkey/fdlibm/e_jn.c +315 -0
- data/vendor/spidermonkey/fdlibm/e_lgamma.c +71 -0
- data/vendor/spidermonkey/fdlibm/e_lgamma_r.c +347 -0
- data/vendor/spidermonkey/fdlibm/e_log.c +184 -0
- data/vendor/spidermonkey/fdlibm/e_log10.c +134 -0
- data/vendor/spidermonkey/fdlibm/e_pow.c +386 -0
- data/vendor/spidermonkey/fdlibm/e_rem_pio2.c +222 -0
- data/vendor/spidermonkey/fdlibm/e_remainder.c +120 -0
- data/vendor/spidermonkey/fdlibm/e_scalb.c +89 -0
- data/vendor/spidermonkey/fdlibm/e_sinh.c +122 -0
- data/vendor/spidermonkey/fdlibm/e_sqrt.c +497 -0
- data/vendor/spidermonkey/fdlibm/fdlibm.h +273 -0
- data/vendor/spidermonkey/fdlibm/fdlibm.mak +1453 -0
- data/vendor/spidermonkey/fdlibm/fdlibm.mdp +0 -0
- data/vendor/spidermonkey/fdlibm/k_cos.c +135 -0
- data/vendor/spidermonkey/fdlibm/k_rem_pio2.c +354 -0
- data/vendor/spidermonkey/fdlibm/k_sin.c +114 -0
- data/vendor/spidermonkey/fdlibm/k_standard.c +785 -0
- data/vendor/spidermonkey/fdlibm/k_tan.c +170 -0
- data/vendor/spidermonkey/fdlibm/s_asinh.c +101 -0
- data/vendor/spidermonkey/fdlibm/s_atan.c +175 -0
- data/vendor/spidermonkey/fdlibm/s_cbrt.c +133 -0
- data/vendor/spidermonkey/fdlibm/s_ceil.c +120 -0
- data/vendor/spidermonkey/fdlibm/s_copysign.c +72 -0
- data/vendor/spidermonkey/fdlibm/s_cos.c +118 -0
- data/vendor/spidermonkey/fdlibm/s_erf.c +356 -0
- data/vendor/spidermonkey/fdlibm/s_expm1.c +267 -0
- data/vendor/spidermonkey/fdlibm/s_fabs.c +70 -0
- data/vendor/spidermonkey/fdlibm/s_finite.c +71 -0
- data/vendor/spidermonkey/fdlibm/s_floor.c +121 -0
- data/vendor/spidermonkey/fdlibm/s_frexp.c +99 -0
- data/vendor/spidermonkey/fdlibm/s_ilogb.c +85 -0
- data/vendor/spidermonkey/fdlibm/s_isnan.c +74 -0
- data/vendor/spidermonkey/fdlibm/s_ldexp.c +66 -0
- data/vendor/spidermonkey/fdlibm/s_lib_version.c +73 -0
- data/vendor/spidermonkey/fdlibm/s_log1p.c +211 -0
- data/vendor/spidermonkey/fdlibm/s_logb.c +79 -0
- data/vendor/spidermonkey/fdlibm/s_matherr.c +64 -0
- data/vendor/spidermonkey/fdlibm/s_modf.c +132 -0
- data/vendor/spidermonkey/fdlibm/s_nextafter.c +124 -0
- data/vendor/spidermonkey/fdlibm/s_rint.c +131 -0
- data/vendor/spidermonkey/fdlibm/s_scalbn.c +107 -0
- data/vendor/spidermonkey/fdlibm/s_signgam.c +40 -0
- data/vendor/spidermonkey/fdlibm/s_significand.c +68 -0
- data/vendor/spidermonkey/fdlibm/s_sin.c +118 -0
- data/vendor/spidermonkey/fdlibm/s_tan.c +112 -0
- data/vendor/spidermonkey/fdlibm/s_tanh.c +122 -0
- data/vendor/spidermonkey/fdlibm/w_acos.c +78 -0
- data/vendor/spidermonkey/fdlibm/w_acosh.c +78 -0
- data/vendor/spidermonkey/fdlibm/w_asin.c +80 -0
- data/vendor/spidermonkey/fdlibm/w_atan2.c +79 -0
- data/vendor/spidermonkey/fdlibm/w_atanh.c +81 -0
- data/vendor/spidermonkey/fdlibm/w_cosh.c +77 -0
- data/vendor/spidermonkey/fdlibm/w_exp.c +88 -0
- data/vendor/spidermonkey/fdlibm/w_fmod.c +78 -0
- data/vendor/spidermonkey/fdlibm/w_gamma.c +85 -0
- data/vendor/spidermonkey/fdlibm/w_gamma_r.c +81 -0
- data/vendor/spidermonkey/fdlibm/w_hypot.c +78 -0
- data/vendor/spidermonkey/fdlibm/w_j0.c +105 -0
- data/vendor/spidermonkey/fdlibm/w_j1.c +106 -0
- data/vendor/spidermonkey/fdlibm/w_jn.c +128 -0
- data/vendor/spidermonkey/fdlibm/w_lgamma.c +85 -0
- data/vendor/spidermonkey/fdlibm/w_lgamma_r.c +81 -0
- data/vendor/spidermonkey/fdlibm/w_log.c +78 -0
- data/vendor/spidermonkey/fdlibm/w_log10.c +81 -0
- data/vendor/spidermonkey/fdlibm/w_pow.c +99 -0
- data/vendor/spidermonkey/fdlibm/w_remainder.c +77 -0
- data/vendor/spidermonkey/fdlibm/w_scalb.c +95 -0
- data/vendor/spidermonkey/fdlibm/w_sinh.c +77 -0
- data/vendor/spidermonkey/fdlibm/w_sqrt.c +77 -0
- data/vendor/spidermonkey/javascript-trace.d +73 -0
- data/vendor/spidermonkey/js.c +3951 -0
- data/vendor/spidermonkey/js.mak +4438 -0
- data/vendor/spidermonkey/js.mdp +0 -0
- data/vendor/spidermonkey/js.msg +307 -0
- data/vendor/spidermonkey/js.pkg +2 -0
- data/vendor/spidermonkey/js3240.rc +79 -0
- data/vendor/spidermonkey/jsOS240.def +654 -0
- data/vendor/spidermonkey/jsapi.c +5836 -0
- data/vendor/spidermonkey/jsapi.h +2624 -0
- data/vendor/spidermonkey/jsarena.c +450 -0
- data/vendor/spidermonkey/jsarena.h +318 -0
- data/vendor/spidermonkey/jsarray.c +2988 -0
- data/vendor/spidermonkey/jsarray.h +124 -0
- data/vendor/spidermonkey/jsatom.c +1045 -0
- data/vendor/spidermonkey/jsatom.h +442 -0
- data/vendor/spidermonkey/jsbit.h +253 -0
- data/vendor/spidermonkey/jsbool.c +176 -0
- data/vendor/spidermonkey/jsbool.h +73 -0
- data/vendor/spidermonkey/jsclist.h +139 -0
- data/vendor/spidermonkey/jscntxt.c +1348 -0
- data/vendor/spidermonkey/jscntxt.h +1120 -0
- data/vendor/spidermonkey/jscompat.h +57 -0
- data/vendor/spidermonkey/jsconfig.h +248 -0
- data/vendor/spidermonkey/jsconfig.mk +181 -0
- data/vendor/spidermonkey/jscpucfg.c +383 -0
- data/vendor/spidermonkey/jscpucfg.h +212 -0
- data/vendor/spidermonkey/jsdate.c +2398 -0
- data/vendor/spidermonkey/jsdate.h +124 -0
- data/vendor/spidermonkey/jsdbgapi.c +1799 -0
- data/vendor/spidermonkey/jsdbgapi.h +464 -0
- data/vendor/spidermonkey/jsdhash.c +868 -0
- data/vendor/spidermonkey/jsdhash.h +592 -0
- data/vendor/spidermonkey/jsdtoa.c +3167 -0
- data/vendor/spidermonkey/jsdtoa.h +130 -0
- data/vendor/spidermonkey/jsdtracef.c +317 -0
- data/vendor/spidermonkey/jsdtracef.h +77 -0
- data/vendor/spidermonkey/jsemit.c +6909 -0
- data/vendor/spidermonkey/jsemit.h +741 -0
- data/vendor/spidermonkey/jsexn.c +1371 -0
- data/vendor/spidermonkey/jsexn.h +96 -0
- data/vendor/spidermonkey/jsfile.c +2736 -0
- data/vendor/spidermonkey/jsfile.h +56 -0
- data/vendor/spidermonkey/jsfile.msg +90 -0
- data/vendor/spidermonkey/jsfun.c +2634 -0
- data/vendor/spidermonkey/jsfun.h +254 -0
- data/vendor/spidermonkey/jsgc.c +3554 -0
- data/vendor/spidermonkey/jsgc.h +403 -0
- data/vendor/spidermonkey/jshash.c +476 -0
- data/vendor/spidermonkey/jshash.h +151 -0
- data/vendor/spidermonkey/jsify.pl +485 -0
- data/vendor/spidermonkey/jsinterp.c +6981 -0
- data/vendor/spidermonkey/jsinterp.h +521 -0
- data/vendor/spidermonkey/jsinvoke.c +43 -0
- data/vendor/spidermonkey/jsiter.c +1067 -0
- data/vendor/spidermonkey/jsiter.h +122 -0
- data/vendor/spidermonkey/jskeyword.tbl +124 -0
- data/vendor/spidermonkey/jskwgen.c +460 -0
- data/vendor/spidermonkey/jslibmath.h +266 -0
- data/vendor/spidermonkey/jslock.c +1309 -0
- data/vendor/spidermonkey/jslock.h +313 -0
- data/vendor/spidermonkey/jslocko.asm +60 -0
- data/vendor/spidermonkey/jslog2.c +94 -0
- data/vendor/spidermonkey/jslong.c +264 -0
- data/vendor/spidermonkey/jslong.h +412 -0
- data/vendor/spidermonkey/jsmath.c +568 -0
- data/vendor/spidermonkey/jsmath.h +57 -0
- data/vendor/spidermonkey/jsnum.c +1228 -0
- data/vendor/spidermonkey/jsnum.h +283 -0
- data/vendor/spidermonkey/jsobj.c +5266 -0
- data/vendor/spidermonkey/jsobj.h +709 -0
- data/vendor/spidermonkey/jsopcode.c +5245 -0
- data/vendor/spidermonkey/jsopcode.h +394 -0
- data/vendor/spidermonkey/jsopcode.tbl +523 -0
- data/vendor/spidermonkey/jsotypes.h +202 -0
- data/vendor/spidermonkey/jsparse.c +6680 -0
- data/vendor/spidermonkey/jsparse.h +511 -0
- data/vendor/spidermonkey/jsprf.c +1262 -0
- data/vendor/spidermonkey/jsprf.h +150 -0
- data/vendor/spidermonkey/jsproto.tbl +128 -0
- data/vendor/spidermonkey/jsprvtd.h +267 -0
- data/vendor/spidermonkey/jspubtd.h +744 -0
- data/vendor/spidermonkey/jsregexp.c +4352 -0
- data/vendor/spidermonkey/jsregexp.h +183 -0
- data/vendor/spidermonkey/jsreops.tbl +145 -0
- data/vendor/spidermonkey/jsscan.c +2003 -0
- data/vendor/spidermonkey/jsscan.h +387 -0
- data/vendor/spidermonkey/jsscope.c +1948 -0
- data/vendor/spidermonkey/jsscope.h +418 -0
- data/vendor/spidermonkey/jsscript.c +1832 -0
- data/vendor/spidermonkey/jsscript.h +287 -0
- data/vendor/spidermonkey/jsshell.msg +50 -0
- data/vendor/spidermonkey/jsstddef.h +83 -0
- data/vendor/spidermonkey/jsstr.c +5004 -0
- data/vendor/spidermonkey/jsstr.h +641 -0
- data/vendor/spidermonkey/jstypes.h +475 -0
- data/vendor/spidermonkey/jsutil.c +345 -0
- data/vendor/spidermonkey/jsutil.h +157 -0
- data/vendor/spidermonkey/jsxdrapi.c +800 -0
- data/vendor/spidermonkey/jsxdrapi.h +218 -0
- data/vendor/spidermonkey/jsxml.c +8471 -0
- data/vendor/spidermonkey/jsxml.h +349 -0
- data/vendor/spidermonkey/lock_SunOS.s +119 -0
- data/vendor/spidermonkey/perfect.js +39 -0
- data/vendor/spidermonkey/plify_jsdhash.sed +36 -0
- data/vendor/spidermonkey/prmjtime.c +846 -0
- data/vendor/spidermonkey/prmjtime.h +103 -0
- data/vendor/spidermonkey/resource.h +15 -0
- data/vendor/spidermonkey/rules.mk +197 -0
- data/vendor/spidermonkey/win32.order +384 -0
- 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
|
+
}
|