Chromium Code Reviews
chromiumcodereview-hr@appspot.gserviceaccount.com (chromiumcodereview-hr) | Please choose your nickname with Settings | Help | Chromium Project | Gerrit Changes | Sign out
(561)

Side by Side Diff: fusl/src/math/expf.c

Issue 1714623002: [fusl] clang-format fusl (Closed) Base URL: git@github.com:domokit/mojo.git@master
Patch Set: headers too Created 4 years, 10 months ago
Use n/p to move between diff chunks; N/P to move between comments. Draft comments are only viewable by you.
Jump to:
View unified diff | Download patch
OLDNEW
1 /* origin: FreeBSD /usr/src/lib/msun/src/e_expf.c */ 1 /* origin: FreeBSD /usr/src/lib/msun/src/e_expf.c */
2 /* 2 /*
3 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com. 3 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
4 */ 4 */
5 /* 5 /*
6 * ==================================================== 6 * ====================================================
7 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 7 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8 * 8 *
9 * Developed at SunPro, a Sun Microsystems, Inc. business. 9 * Developed at SunPro, a Sun Microsystems, Inc. business.
10 * Permission to use, copy, modify, and distribute this 10 * Permission to use, copy, modify, and distribute this
11 * software is freely granted, provided that this notice 11 * software is freely granted, provided that this notice
12 * is preserved. 12 * is preserved.
13 * ==================================================== 13 * ====================================================
14 */ 14 */
15 15
16 #include "libm.h" 16 #include "libm.h"
17 17
18 static const float 18 static const float half[2] = {0.5, -0.5},
19 half[2] = {0.5,-0.5}, 19 ln2hi = 6.9314575195e-1f, /* 0x3f317200 */
20 ln2hi = 6.9314575195e-1f, /* 0x3f317200 */ 20 ln2lo = 1.4286067653e-6f, /* 0x35bfbe8e */
21 ln2lo = 1.4286067653e-6f, /* 0x35bfbe8e */ 21 invln2 = 1.4426950216e+0f, /* 0x3fb8aa3b */
22 invln2 = 1.4426950216e+0f, /* 0x3fb8aa3b */ 22 /*
23 /* 23 * Domain [-0.34568, 0.34568], rang e ~[-4.278e-9, 4.447e-9]:
24 * Domain [-0.34568, 0.34568], range ~[-4.278e-9, 4.447e-9]: 24 * |x*(exp(x)+1)/(exp(x)-1) - p(x)| < 2**-27.74
25 * |x*(exp(x)+1)/(exp(x)-1) - p(x)| < 2**-27.74 25 */
26 */ 26 P1 = 1.6666625440e-1f, /* 0xaaaa8f.0p-26 */
27 P1 = 1.6666625440e-1f, /* 0xaaaa8f.0p-26 */ 27 P2 = -2.7667332906e-3f; /* -0xb55215.0p-32 */
28 P2 = -2.7667332906e-3f; /* -0xb55215.0p-32 */
29 28
30 float expf(float x) 29 float expf(float x) {
31 { 30 float_t hi, lo, c, xx, y;
32 » float_t hi, lo, c, xx, y; 31 int k, sign;
33 » int k, sign; 32 uint32_t hx;
34 » uint32_t hx;
35 33
36 » GET_FLOAT_WORD(hx, x); 34 GET_FLOAT_WORD(hx, x);
37 » sign = hx >> 31; /* sign bit of x */ 35 sign = hx >> 31; /* sign bit of x */
38 » hx &= 0x7fffffff; /* high word of |x| */ 36 hx &= 0x7fffffff; /* high word of |x| */
39 37
40 » /* special cases */ 38 /* special cases */
41 » if (hx >= 0x42aeac50) { /* if |x| >= -87.33655f or NaN */ 39 if (hx >= 0x42aeac50) { /* if |x| >= -87.33655f or NaN */
42 » » if (hx >= 0x42b17218 && !sign) { /* x >= 88.722839f */ 40 if (hx >= 0x42b17218 && !sign) { /* x >= 88.722839f */
43 » » » /* overflow */ 41 /* overflow */
44 » » » x *= 0x1p127f; 42 x *= 0x1p127f;
45 » » » return x; 43 return x;
46 » » } 44 }
47 » » if (sign) { 45 if (sign) {
48 » » » /* underflow */ 46 /* underflow */
49 » » » FORCE_EVAL(-0x1p-149f/x); 47 FORCE_EVAL(-0x1p-149f / x);
50 » » » if (hx >= 0x42cff1b5) /* x <= -103.972084f */ 48 if (hx >= 0x42cff1b5) /* x <= -103.972084f */
51 » » » » return 0; 49 return 0;
52 » » } 50 }
53 » } 51 }
54 52
55 » /* argument reduction */ 53 /* argument reduction */
56 » if (hx > 0x3eb17218) { /* if |x| > 0.5 ln2 */ 54 if (hx > 0x3eb17218) { /* if |x| > 0.5 ln2 */
57 » » if (hx > 0x3f851592) /* if |x| > 1.5 ln2 */ 55 if (hx > 0x3f851592) /* if |x| > 1.5 ln2 */
58 » » » k = invln2*x + half[sign]; 56 k = invln2 * x + half[sign];
59 » » else 57 else
60 » » » k = 1 - sign - sign; 58 k = 1 - sign - sign;
61 » » hi = x - k*ln2hi; /* k*ln2hi is exact here */ 59 hi = x - k * ln2hi; /* k*ln2hi is exact here */
62 » » lo = k*ln2lo; 60 lo = k * ln2lo;
63 » » x = hi - lo; 61 x = hi - lo;
64 » } else if (hx > 0x39000000) { /* |x| > 2**-14 */ 62 } else if (hx > 0x39000000) { /* |x| > 2**-14 */
65 » » k = 0; 63 k = 0;
66 » » hi = x; 64 hi = x;
67 » » lo = 0; 65 lo = 0;
68 » } else { 66 } else {
69 » » /* raise inexact */ 67 /* raise inexact */
70 » » FORCE_EVAL(0x1p127f + x); 68 FORCE_EVAL(0x1p127f + x);
71 » » return 1 + x; 69 return 1 + x;
72 » } 70 }
73 71
74 » /* x is now in primary range */ 72 /* x is now in primary range */
75 » xx = x*x; 73 xx = x * x;
76 » c = x - xx*(P1+xx*P2); 74 c = x - xx * (P1 + xx * P2);
77 » y = 1 + (x*c/(2-c) - lo + hi); 75 y = 1 + (x * c / (2 - c) - lo + hi);
78 » if (k == 0) 76 if (k == 0)
79 » » return y; 77 return y;
80 » return scalbnf(y, k); 78 return scalbnf(y, k);
81 } 79 }
OLDNEW

Powered by Google App Engine
This is Rietveld 408576698