diff options
author | Szabolcs Nagy <szabolcs.nagy@arm.com> | 2018-06-29 09:56:54 +0100 |
---|---|---|
committer | Szabolcs Nagy <szabolcs.nagy@arm.com> | 2018-06-29 13:25:58 +0100 |
commit | 5e8389113b47a48622d1de8beb9a4adc2b332694 (patch) | |
tree | 88e1303e741d50737c310a13e3fa49a2d2c92d37 /math/exp.c | |
parent | 4aa92161d340899eaf3cd1b599a72109ba818c22 (diff) | |
download | arm-optimized-routines-5e8389113b47a48622d1de8beb9a4adc2b332694.tar.gz |
Fix GNU style issues
Whitespace changes only.
Diffstat (limited to 'math/exp.c')
-rw-r--r-- | math/exp.c | 12 |
1 files changed, 6 insertions, 6 deletions
@@ -124,27 +124,27 @@ exp_inline (double x, double xtail, int hastail) ki = asuint64 (kd); kd -= Shift; #endif - r = x + kd*NegLn2hiN + kd*NegLn2loN; + r = x + kd * NegLn2hiN + kd * NegLn2loN; /* The code assumes 2^-200 < |xtail| < 2^-8/N. */ if (hastail) r += xtail; /* 2^(k/N) ~= scale * (1 + tail). */ - idx = 2*(ki % N); + idx = 2 * (ki % N); top = ki << (52 - EXP_TABLE_BITS); tail = asdouble (T[idx]); /* This is only a valid scale when -1023*N < k < 1024*N. */ sbits = T[idx + 1] + top; /* exp(x) = 2^(k/N) * exp(r) ~= scale + scale * (tail + exp(r) - 1). */ /* Evaluation is optimized assuming superscalar pipelined execution. */ - r2 = r*r; + r2 = r * r; /* Without fma the worst case error is 0.25/N ulp larger. */ /* Worst case error is less than 0.5+1.11/N+(abs poly error * 2^53) ulp. */ #if EXP_POLY_ORDER == 4 - tmp = tail + r + r2*C2 + r*r2*(C3 + r*C4); + tmp = tail + r + r2 * C2 + r * r2 * (C3 + r * C4); #elif EXP_POLY_ORDER == 5 - tmp = tail + r + r2*(C2 + r*C3) + r2*r2*(C4 + r*C5); + tmp = tail + r + r2 * (C2 + r * C3) + r2 * r2 * (C4 + r * C5); #elif EXP_POLY_ORDER == 6 - tmp = tail + r + r2*(0.5 + r*C3) + r2*r2*(C4 + r*C5 + r2*C6); + tmp = tail + r + r2 * (0.5 + r * C3) + r2 * r2 * (C4 + r * C5 + r2 * C6); #endif if (unlikely (abstop == 0)) return specialcase (tmp, sbits, ki); |