/* Return the quantum 1 * 10^exponent(x) */
_Bool __samequantumd64 (_Decimal64 x, _Decimal64 y)
{
int xnan = __isnand64 (x);
int ynan = __isnand64 (y);
if (xnan | ynan)
return xnan & ynan;
return __quantumd64 (x) == __quantumd64 (y);
}
CMPINT
__BACKEND_(unorddd2) (_Decimal64 x, _Decimal64 y)
{
return __isnand64 (x) || isnand64 (y);
}
int
__BACKEND_(unordsd2) (_Decimal32 x, _Decimal32 y)
{
return __isnand64 ((_Decimal64)x) || __isnand64 ((_Decimal64)y);
}
static DEC_TYPE
IEEE_FUNCTION_NAME (DEC_TYPE val)
{
_Decimal64 result, top, fraction, a, x, tp;
_Decimal64 t1, t2, t3, t4, t5, t6, t7, t8;
_Decimal64 t9, t10, t11, t12, t13, t14, t15, t16;
int exp;
long tmp;
if (__isnand64(val))
return val+val;
if (val == DFP_CONSTANT(0.0))
{
DFP_EXCEPT (FE_DIVBYZERO);
return -DFP_HUGE_VAL;
}
if (val < DFP_CONSTANT(0.0))
{
DFP_EXCEPT (FE_INVALID);
return DFP_NAN;
}
if (__isinfd64(val))
return val;
result = __frexpd64 (val, &exp);
tmp = result * 100.0DD;
top = tmp;
top = top * 0.01DD;
fraction = result - top;
if (fraction != 0.00DD)
{
a = top;
x = result;
t1 = (x - a) / a;
tp = t1 * t1; /* tp == ((x-a)/a)^2 */
t2 = tp / 2.0DD; /* t2 == (((x-a)/a)^2)/2 */
tp = tp * t1; /* tp == ((x-a)/a)^3 */
t3 = tp / 3.0DD; /* t3 == (((x-a)/a)^3)/3 */
tp = tp * t1; /* tp == ((x-a)/a)^4 */
t4 = tp / 4.0DD; /* t4 == (((x-a)/a)^4)/4 */
tp = tp * t1; /* tp == ((x-a)/a)^5 */
t5 = tp / 5.0DD; /* t5 == (((x-a)/a)^5)/5 */
tp = tp * t1; /* tp == ((x-a)/a)^6 */
t6 = tp / 6.0DD; /* t6 == (((x-a)/a)^6)/6 */
tp = tp * t1; /* tp == ((x-a)/a)^7 */
t7 = tp / 7.0DD; /* t7 == (((x-a)/a)^7)/7 */
tp = tp * t1; /* tp == ((x-a)/a)^8 */
t8 = tp / 8.0DD; /* t8 == (((x-a)/a)^8)/8 */
if ( t8 > 1.0E-16DD)
{
tp = tp * t1; /* tp == ((x-a)/a)^9 */
t9 = tp / 9.0DD; /* t9 == (((x-a)/a)^9)/9 */
tp = tp * t1; /* tp == ((x-a)/a)^10 */
t10 = tp / 10.0DD; /* t10 == (((x-a)/a)^10)/10 */
tp = tp * t1; /* tp == ((x-a)/a)^11 */
t11 = tp / 11.0DD; /* t12 == (((x-a)/a)^11)/11 */
tp = tp * t1; /* tp == ((x-a)/a)^12 */
t12 = tp / 12.0DD; /* t12 == (((x-a)/a)^11)/12 */
if (t12 > 1.0E-16DD)
{
tp = tp * t1; /* tp == ((x-a)/a)^13 */
t13 = tp / 13.0DD; /* t13 == (((x-a)/a)^13)/13 */
tp = tp * t1; /* tp == ((x-a)/a)^14 */
t14 = tp / 14.0DD; /* t14 == (((x-a)/a)^14)/14 */
tp = tp * t1; /* tp == ((x-a)/a)^15 */
t15 = tp / 15.0DD; /* t15 == (((x-a)/a)^15)/15 */
tp = tp * t1; /* tp == ((x-a)/a)^16 */
t16 = tp / 16.0DD; /* t16 == (((x-a)/a)^16)/16 */
tp = t15 - t16;
tp = tp + (t13 - t14);
tp = tp + (t11 - t12);
}
else
{
tp = t11 - t12;
}
tp = tp + (t9 - t10);
tp = tp + (t7 - t8);
}
else
tp = t7 - t8;
/* now sum the terms from smallest to largest to avoid loss of percision */
tp = tp + (t5 - t6);
tp = tp + (t3 - t4);
tp = tp + (t1 - t2);
if (exp!=0)
result = (__LN_10 * exp) + __LN2digits[tmp] + tp;
else
result = __LN2digits[tmp] + tp;
}
else
{
if (exp!=0)
result = (__LN_10 * exp) + __LN2digits[tmp];
else
result = __LN2digits[tmp];
}
return result;
}
