static long
internal_equality (long al, long ah, long bl, long bh, long neq)
{
FP_DECL_EX;
FP_DECL_Q(A); FP_DECL_Q(B);
long r;
FP_UNPACK_RAW_Q(A, a);
FP_UNPACK_RAW_Q(B, b);
if ((A_e == _FP_EXPMAX_Q && !_FP_FRAC_ZEROP_2(A))
|| (B_e == _FP_EXPMAX_Q && !_FP_FRAC_ZEROP_2(B)))
{
/* EQ and NE signal invalid operation only if either operand is SNaN. */
if (FP_ISSIGNAN_Q(A) || FP_ISSIGNAN_Q(B))
{
FP_SET_EXCEPTION(FP_EX_INVALID);
FP_HANDLE_EXCEPTIONS;
}
return -1;
}
r = (A_e == B_e
&& _FP_FRAC_EQ_2 (A, B)
&& (A_s == B_s || (!A_e && _FP_FRAC_ZEROP_2(A))));
r ^= neq;
return r;
}
int __ieee754_ilogbl (long double x)
{
FP_DECL_EX;
FP_DECL_Q(X);
/*
FP_UNPACK_Q(X, x);
switch (X_c)
{
case FP_CLS_ZERO:
return FP_ILOGB0;
case FP_CLS_NAN:
case FP_CLS_INF:
return FP_ILOGBNAN;
default:
return X_e;
}
*/
FP_UNPACK_RAW_Q(X, x);
switch (X_e)
{
default:
return X_e - _FP_EXPBIAS_Q;
case 0:
#if (2 * _FP_W_TYPE_SIZE) < _FP_FRACBITS_Q
if (_FP_FRAC_ZEROP_4(X))
return FP_ILOGB0;
else
{
_FP_I_TYPE shift;
_FP_FRAC_CLZ_4(shift, X);
shift -= _FP_FRACXBITS_Q;
return X_e - _FP_EXPBIAS_Q - 1 + shift;
}
#else
if (_FP_FRAC_ZEROP_2(X))
return FP_ILOGB0;
else
{
_FP_I_TYPE shift;
_FP_FRAC_CLZ_2(shift, X);
shift -= _FP_FRACXBITS_Q;
return X_e - _FP_EXPBIAS_Q - 1 + shift;
}
#endif
case _FP_EXPBIAS_Q:
return FP_ILOGBNAN;
}
}