z = absZ;
if ( zSign ) z = - z;
if ( ( absZ>>32 ) || ( z && ( ( z < 0 ) ^ zSign ) ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return zSign ? (int32_t) 0x80000000 : 0x7FFFFFFF;
}
if ( roundBits ) STATUS(float_exception_flags) |= float_flag_inexact;
if ( zSign ) z = - z;
if ( z && ( ( z < 0 ) ^ zSign ) ) {
overflow:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return
zSign ? (int64_t) LIT64( 0x8000000000000000 )
: LIT64( 0x7FFFFFFFFFFFFFFF );
if (increment) {
++absZ0;
if (absZ0 == 0) {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return LIT64(0xFFFFFFFFFFFFFFFF);
}
absZ0 &= ~(((uint64_t)(absZ1<<1) == 0) & roundNearestEven);
}
if (zSign && absZ0) {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
{
if (STATUS(flush_inputs_to_zero)) {
if (extractFloat32Exp(a) == 0 && extractFloat32Frac(a) != 0) {
- float_raise(float_flag_input_denormal STATUS_VAR);
+ float_raise(float_flag_input_denormal, status);
return make_float32(float32_val(a) & 0x80000000);
}
}
|| ( ( zExp == 0xFD )
&& ( (int32_t) ( zSig + roundIncrement ) < 0 ) )
) {
- float_raise( float_flag_overflow | float_flag_inexact STATUS_VAR);
+ float_raise(float_flag_overflow | float_flag_inexact, status);
return packFloat32( zSign, 0xFF, - ( roundIncrement == 0 ));
}
if ( zExp < 0 ) {
if (STATUS(flush_to_zero)) {
- float_raise(float_flag_output_denormal STATUS_VAR);
+ float_raise(float_flag_output_denormal, status);
return packFloat32(zSign, 0, 0);
}
isTiny =
shift32RightJamming( zSig, - zExp, &zSig );
zExp = 0;
roundBits = zSig & 0x7F;
- if ( isTiny && roundBits ) float_raise( float_flag_underflow STATUS_VAR);
+ if (isTiny && roundBits) {
+ float_raise(float_flag_underflow, status);
+ }
}
}
if ( roundBits ) STATUS(float_exception_flags) |= float_flag_inexact;
int8 shiftCount;
shiftCount = countLeadingZeros32( zSig ) - 1;
- return roundAndPackFloat32( zSign, zExp - shiftCount, zSig<<shiftCount STATUS_VAR);
+ return roundAndPackFloat32(zSign, zExp - shiftCount, zSig<<shiftCount,
+ status);
}
{
if (STATUS(flush_inputs_to_zero)) {
if (extractFloat64Exp(a) == 0 && extractFloat64Frac(a) != 0) {
- float_raise(float_flag_input_denormal STATUS_VAR);
+ float_raise(float_flag_input_denormal, status);
return make_float64(float64_val(a) & (1ULL << 63));
}
}
|| ( ( zExp == 0x7FD )
&& ( (int64_t) ( zSig + roundIncrement ) < 0 ) )
) {
- float_raise( float_flag_overflow | float_flag_inexact STATUS_VAR);
+ float_raise(float_flag_overflow | float_flag_inexact, status);
return packFloat64( zSign, 0x7FF, - ( roundIncrement == 0 ));
}
if ( zExp < 0 ) {
if (STATUS(flush_to_zero)) {
- float_raise(float_flag_output_denormal STATUS_VAR);
+ float_raise(float_flag_output_denormal, status);
return packFloat64(zSign, 0, 0);
}
isTiny =
shift64RightJamming( zSig, - zExp, &zSig );
zExp = 0;
roundBits = zSig & 0x3FF;
- if ( isTiny && roundBits ) float_raise( float_flag_underflow STATUS_VAR);
+ if (isTiny && roundBits) {
+ float_raise(float_flag_underflow, status);
+ }
}
}
if ( roundBits ) STATUS(float_exception_flags) |= float_flag_inexact;
int8 shiftCount;
shiftCount = countLeadingZeros64( zSig ) - 1;
- return roundAndPackFloat64( zSign, zExp - shiftCount, zSig<<shiftCount STATUS_VAR);
+ return roundAndPackFloat64(zSign, zExp - shiftCount, zSig<<shiftCount,
+ status);
}
}
if ( zExp <= 0 ) {
if (STATUS(flush_to_zero)) {
- float_raise(float_flag_output_denormal STATUS_VAR);
+ float_raise(float_flag_output_denormal, status);
return packFloatx80(zSign, 0, 0);
}
isTiny =
shift64RightJamming( zSig0, 1 - zExp, &zSig0 );
zExp = 0;
roundBits = zSig0 & roundMask;
- if ( isTiny && roundBits ) float_raise( float_flag_underflow STATUS_VAR);
+ if (isTiny && roundBits) {
+ float_raise(float_flag_underflow, status);
+ }
if ( roundBits ) STATUS(float_exception_flags) |= float_flag_inexact;
zSig0 += roundIncrement;
if ( (int64_t) zSig0 < 0 ) zExp = 1;
) {
roundMask = 0;
overflow:
- float_raise( float_flag_overflow | float_flag_inexact STATUS_VAR);
+ float_raise(float_flag_overflow | float_flag_inexact, status);
if ( ( roundingMode == float_round_to_zero )
|| ( zSign && ( roundingMode == float_round_up ) )
|| ( ! zSign && ( roundingMode == float_round_down ) )
|| ( zSig0 < LIT64( 0xFFFFFFFFFFFFFFFF ) );
shift64ExtraRightJamming( zSig0, zSig1, 1 - zExp, &zSig0, &zSig1 );
zExp = 0;
- if ( isTiny && zSig1 ) float_raise( float_flag_underflow STATUS_VAR);
+ if (isTiny && zSig1) {
+ float_raise(float_flag_underflow, status);
+ }
if ( zSig1 ) STATUS(float_exception_flags) |= float_flag_inexact;
switch (roundingMode) {
case float_round_nearest_even:
shiftCount = countLeadingZeros64( zSig0 );
shortShift128Left( zSig0, zSig1, shiftCount, &zSig0, &zSig1 );
zExp -= shiftCount;
- return
- roundAndPackFloatx80( roundingPrecision, zSign, zExp, zSig0, zSig1 STATUS_VAR);
+ return roundAndPackFloatx80(roundingPrecision, zSign, zExp,
+ zSig0, zSig1, status);
}
&& increment
)
) {
- float_raise( float_flag_overflow | float_flag_inexact STATUS_VAR);
+ float_raise(float_flag_overflow | float_flag_inexact, status);
if ( ( roundingMode == float_round_to_zero )
|| ( zSign && ( roundingMode == float_round_up ) )
|| ( ! zSign && ( roundingMode == float_round_down ) )
}
if ( zExp < 0 ) {
if (STATUS(flush_to_zero)) {
- float_raise(float_flag_output_denormal STATUS_VAR);
+ float_raise(float_flag_output_denormal, status);
return packFloat128(zSign, 0, 0, 0);
}
isTiny =
shift128ExtraRightJamming(
zSig0, zSig1, zSig2, - zExp, &zSig0, &zSig1, &zSig2 );
zExp = 0;
- if ( isTiny && zSig2 ) float_raise( float_flag_underflow STATUS_VAR);
+ if (isTiny && zSig2) {
+ float_raise(float_flag_underflow, status);
+ }
switch (roundingMode) {
case float_round_nearest_even:
case float_round_ties_away:
zSig0, zSig1, 0, - shiftCount, &zSig0, &zSig1, &zSig2 );
}
zExp -= shiftCount;
- return roundAndPackFloat128( zSign, zExp, zSig0, zSig1, zSig2 STATUS_VAR);
+ return roundAndPackFloat128(zSign, zExp, zSig0, zSig1, zSig2, status);
}
if ( a == 0 ) return float32_zero;
if ( a == (int32_t) 0x80000000 ) return packFloat32( 1, 0x9E, 0 );
zSign = ( a < 0 );
- return normalizeRoundAndPackFloat32( zSign, 0x9C, zSign ? - a : a STATUS_VAR );
-
+ return normalizeRoundAndPackFloat32(zSign, 0x9C, zSign ? -a : a, status);
}
/*----------------------------------------------------------------------------
else {
absA <<= shiftCount;
}
- return roundAndPackFloat32( zSign, 0x9C - shiftCount, absA STATUS_VAR );
+ return roundAndPackFloat32(zSign, 0x9C - shiftCount, absA, status);
}
}
return packFloat64( 1, 0x43E, 0 );
}
zSign = ( a < 0 );
- return normalizeRoundAndPackFloat64( zSign, 0x43C, zSign ? - a : a STATUS_VAR );
-
+ return normalizeRoundAndPackFloat64(zSign, 0x43C, zSign ? -a : a, status);
}
/*----------------------------------------------------------------------------
a <<= shiftcount;
}
- return roundAndPackFloat32(0, 0x9c - shiftcount, a STATUS_VAR);
+ return roundAndPackFloat32(0, 0x9c - shiftcount, a, status);
}
/*----------------------------------------------------------------------------
} else {
a <<= shiftcount;
}
- return roundAndPackFloat64(0, exp - shiftcount, a STATUS_VAR);
+ return roundAndPackFloat64(0, exp - shiftcount, a, status);
}
/*----------------------------------------------------------------------------
if (a == 0) {
return float128_zero;
}
- return normalizeRoundAndPackFloat128(0, 0x406E, a, 0 STATUS_VAR);
+ return normalizeRoundAndPackFloat128(0, 0x406E, a, 0, status);
}
/*----------------------------------------------------------------------------
uint32_t aSig;
uint64_t aSig64;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
aSign = extractFloat32Sign( a );
aSig64 = aSig;
aSig64 <<= 32;
if ( 0 < shiftCount ) shift64RightJamming( aSig64, shiftCount, &aSig64 );
- return roundAndPackInt32( aSign, aSig64 STATUS_VAR );
+ return roundAndPackInt32(aSign, aSig64, status);
}
int_fast16_t aExp, shiftCount;
uint32_t aSig;
int32_t z;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
shiftCount = aExp - 0x9E;
if ( 0 <= shiftCount ) {
if ( float32_val(a) != 0xCF000000 ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
if ( ! aSign || ( ( aExp == 0xFF ) && aSig ) ) return 0x7FFFFFFF;
}
return (int32_t) 0x80000000;
shiftCount = aExp - 0x8E;
if ( 0 <= shiftCount ) {
if ( float32_val(a) != 0xC7000000 ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
if ( ! aSign || ( ( aExp == 0xFF ) && aSig ) ) {
return 0x7FFF;
}
int_fast16_t aExp, shiftCount;
uint32_t aSig;
uint64_t aSig64, aSigExtra;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
aSign = extractFloat32Sign( a );
shiftCount = 0xBE - aExp;
if ( shiftCount < 0 ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
if ( ! aSign || ( ( aExp == 0xFF ) && aSig ) ) {
return LIT64( 0x7FFFFFFFFFFFFFFF );
}
aSig64 = aSig;
aSig64 <<= 40;
shift64ExtraRightJamming( aSig64, 0, shiftCount, &aSig64, &aSigExtra );
- return roundAndPackInt64( aSign, aSig64, aSigExtra STATUS_VAR );
+ return roundAndPackInt64(aSign, aSig64, aSigExtra, status);
}
int_fast16_t aExp, shiftCount;
uint32_t aSig;
uint64_t aSig64, aSigExtra;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac(a);
aExp = extractFloat32Exp(a);
aSign = extractFloat32Sign(a);
if ((aSign) && (aExp > 126)) {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
if (float32_is_any_nan(a)) {
return LIT64(0xFFFFFFFFFFFFFFFF);
} else {
aSig |= 0x00800000;
}
if (shiftCount < 0) {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return LIT64(0xFFFFFFFFFFFFFFFF);
}
aSig64 = aSig;
aSig64 <<= 40;
shift64ExtraRightJamming(aSig64, 0, shiftCount, &aSig64, &aSigExtra);
- return roundAndPackUint64(aSign, aSig64, aSigExtra STATUS_VAR);
+ return roundAndPackUint64(aSign, aSig64, aSigExtra, status);
}
/*----------------------------------------------------------------------------
uint64 float32_to_uint64_round_to_zero(float32 a, float_status *status)
{
signed char current_rounding_mode = STATUS(float_rounding_mode);
- set_float_rounding_mode(float_round_to_zero STATUS_VAR);
- int64_t v = float32_to_uint64(a STATUS_VAR);
- set_float_rounding_mode(current_rounding_mode STATUS_VAR);
+ set_float_rounding_mode(float_round_to_zero, status);
+ int64_t v = float32_to_uint64(a, status);
+ set_float_rounding_mode(current_rounding_mode, status);
return v;
}
uint32_t aSig;
uint64_t aSig64;
int64 z;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
shiftCount = aExp - 0xBE;
if ( 0 <= shiftCount ) {
if ( float32_val(a) != 0xDF000000 ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
if ( ! aSign || ( ( aExp == 0xFF ) && aSig ) ) {
return LIT64( 0x7FFFFFFFFFFFFFFF );
}
flag aSign;
int_fast16_t aExp;
uint32_t aSig;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
aSign = extractFloat32Sign( a );
if ( aExp == 0xFF ) {
- if ( aSig ) return commonNaNToFloat64( float32ToCommonNaN( a STATUS_VAR ) STATUS_VAR );
+ if (aSig) {
+ return commonNaNToFloat64(float32ToCommonNaN(a, status), status);
+ }
return packFloat64( aSign, 0x7FF, 0 );
}
if ( aExp == 0 ) {
int_fast16_t aExp;
uint32_t aSig;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
aSign = extractFloat32Sign( a );
if ( aExp == 0xFF ) {
- if ( aSig ) return commonNaNToFloatx80( float32ToCommonNaN( a STATUS_VAR ) STATUS_VAR );
+ if (aSig) {
+ return commonNaNToFloatx80(float32ToCommonNaN(a, status), status);
+ }
return packFloatx80( aSign, 0x7FFF, LIT64( 0x8000000000000000 ) );
}
if ( aExp == 0 ) {
int_fast16_t aExp;
uint32_t aSig;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
aSign = extractFloat32Sign( a );
if ( aExp == 0xFF ) {
- if ( aSig ) return commonNaNToFloat128( float32ToCommonNaN( a STATUS_VAR ) STATUS_VAR );
+ if (aSig) {
+ return commonNaNToFloat128(float32ToCommonNaN(a, status), status);
+ }
return packFloat128( aSign, 0x7FFF, 0, 0 );
}
if ( aExp == 0 ) {
int_fast16_t aExp;
uint32_t lastBitMask, roundBitsMask;
uint32_t z;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aExp = extractFloat32Exp( a );
if ( 0x96 <= aExp ) {
if ( ( aExp == 0xFF ) && extractFloat32Frac( a ) ) {
- return propagateFloat32NaN( a, a STATUS_VAR );
+ return propagateFloat32NaN(a, a, status);
}
return a;
}
bSig <<= 6;
if ( 0 < expDiff ) {
if ( aExp == 0xFF ) {
- if ( aSig ) return propagateFloat32NaN( a, b STATUS_VAR );
+ if (aSig) {
+ return propagateFloat32NaN(a, b, status);
+ }
return a;
}
if ( bExp == 0 ) {
}
else if ( expDiff < 0 ) {
if ( bExp == 0xFF ) {
- if ( bSig ) return propagateFloat32NaN( a, b STATUS_VAR );
+ if (bSig) {
+ return propagateFloat32NaN(a, b, status);
+ }
return packFloat32( zSign, 0xFF, 0 );
}
if ( aExp == 0 ) {
}
else {
if ( aExp == 0xFF ) {
- if ( aSig | bSig ) return propagateFloat32NaN( a, b STATUS_VAR );
+ if (aSig | bSig) {
+ return propagateFloat32NaN(a, b, status);
+ }
return a;
}
if ( aExp == 0 ) {
if (STATUS(flush_to_zero)) {
if (aSig | bSig) {
- float_raise(float_flag_output_denormal STATUS_VAR);
+ float_raise(float_flag_output_denormal, status);
}
return packFloat32(zSign, 0, 0);
}
++zExp;
}
roundAndPack:
- return roundAndPackFloat32( zSign, zExp, zSig STATUS_VAR );
+ return roundAndPackFloat32(zSign, zExp, zSig, status);
}
if ( 0 < expDiff ) goto aExpBigger;
if ( expDiff < 0 ) goto bExpBigger;
if ( aExp == 0xFF ) {
- if ( aSig | bSig ) return propagateFloat32NaN( a, b STATUS_VAR );
- float_raise( float_flag_invalid STATUS_VAR);
+ if (aSig | bSig) {
+ return propagateFloat32NaN(a, b, status);
+ }
+ float_raise(float_flag_invalid, status);
return float32_default_nan;
}
if ( aExp == 0 ) {
return packFloat32( STATUS(float_rounding_mode) == float_round_down, 0, 0 );
bExpBigger:
if ( bExp == 0xFF ) {
- if ( bSig ) return propagateFloat32NaN( a, b STATUS_VAR );
+ if (bSig) {
+ return propagateFloat32NaN(a, b, status);
+ }
return packFloat32( zSign ^ 1, 0xFF, 0 );
}
if ( aExp == 0 ) {
goto normalizeRoundAndPack;
aExpBigger:
if ( aExp == 0xFF ) {
- if ( aSig ) return propagateFloat32NaN( a, b STATUS_VAR );
+ if (aSig) {
+ return propagateFloat32NaN(a, b, status);
+ }
return a;
}
if ( bExp == 0 ) {
zExp = aExp;
normalizeRoundAndPack:
--zExp;
- return normalizeRoundAndPackFloat32( zSign, zExp, zSig STATUS_VAR );
+ return normalizeRoundAndPackFloat32(zSign, zExp, zSig, status);
}
float32 float32_add(float32 a, float32 b, float_status *status)
{
flag aSign, bSign;
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
aSign = extractFloat32Sign( a );
bSign = extractFloat32Sign( b );
if ( aSign == bSign ) {
- return addFloat32Sigs( a, b, aSign STATUS_VAR);
+ return addFloat32Sigs(a, b, aSign, status);
}
else {
- return subFloat32Sigs( a, b, aSign STATUS_VAR );
+ return subFloat32Sigs(a, b, aSign, status);
}
}
float32 float32_sub(float32 a, float32 b, float_status *status)
{
flag aSign, bSign;
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
aSign = extractFloat32Sign( a );
bSign = extractFloat32Sign( b );
if ( aSign == bSign ) {
- return subFloat32Sigs( a, b, aSign STATUS_VAR );
+ return subFloat32Sigs(a, b, aSign, status);
}
else {
- return addFloat32Sigs( a, b, aSign STATUS_VAR );
+ return addFloat32Sigs(a, b, aSign, status);
}
}
uint64_t zSig64;
uint32_t zSig;
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
zSign = aSign ^ bSign;
if ( aExp == 0xFF ) {
if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {
- return propagateFloat32NaN( a, b STATUS_VAR );
+ return propagateFloat32NaN(a, b, status);
}
if ( ( bExp | bSig ) == 0 ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float32_default_nan;
}
return packFloat32( zSign, 0xFF, 0 );
}
if ( bExp == 0xFF ) {
- if ( bSig ) return propagateFloat32NaN( a, b STATUS_VAR );
+ if (bSig) {
+ return propagateFloat32NaN(a, b, status);
+ }
if ( ( aExp | aSig ) == 0 ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float32_default_nan;
}
return packFloat32( zSign, 0xFF, 0 );
zSig <<= 1;
--zExp;
}
- return roundAndPackFloat32( zSign, zExp, zSig STATUS_VAR );
+ return roundAndPackFloat32(zSign, zExp, zSig, status);
}
flag aSign, bSign, zSign;
int_fast16_t aExp, bExp, zExp;
uint32_t aSig, bSig, zSig;
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
bSign = extractFloat32Sign( b );
zSign = aSign ^ bSign;
if ( aExp == 0xFF ) {
- if ( aSig ) return propagateFloat32NaN( a, b STATUS_VAR );
+ if (aSig) {
+ return propagateFloat32NaN(a, b, status);
+ }
if ( bExp == 0xFF ) {
- if ( bSig ) return propagateFloat32NaN( a, b STATUS_VAR );
- float_raise( float_flag_invalid STATUS_VAR);
+ if (bSig) {
+ return propagateFloat32NaN(a, b, status);
+ }
+ float_raise(float_flag_invalid, status);
return float32_default_nan;
}
return packFloat32( zSign, 0xFF, 0 );
}
if ( bExp == 0xFF ) {
- if ( bSig ) return propagateFloat32NaN( a, b STATUS_VAR );
+ if (bSig) {
+ return propagateFloat32NaN(a, b, status);
+ }
return packFloat32( zSign, 0, 0 );
}
if ( bExp == 0 ) {
if ( bSig == 0 ) {
if ( ( aExp | aSig ) == 0 ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float32_default_nan;
}
- float_raise( float_flag_divbyzero STATUS_VAR);
+ float_raise(float_flag_divbyzero, status);
return packFloat32( zSign, 0xFF, 0 );
}
normalizeFloat32Subnormal( bSig, &bExp, &bSig );
if ( ( zSig & 0x3F ) == 0 ) {
zSig |= ( (uint64_t) bSig * zSig != ( (uint64_t) aSig )<<32 );
}
- return roundAndPackFloat32( zSign, zExp, zSig STATUS_VAR );
+ return roundAndPackFloat32(zSign, zExp, zSig, status);
}
uint64_t aSig64, bSig64, q64;
uint32_t alternateASig;
int32_t sigMean;
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
bExp = extractFloat32Exp( b );
if ( aExp == 0xFF ) {
if ( aSig || ( ( bExp == 0xFF ) && bSig ) ) {
- return propagateFloat32NaN( a, b STATUS_VAR );
+ return propagateFloat32NaN(a, b, status);
}
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float32_default_nan;
}
if ( bExp == 0xFF ) {
- if ( bSig ) return propagateFloat32NaN( a, b STATUS_VAR );
+ if (bSig) {
+ return propagateFloat32NaN(a, b, status);
+ }
return a;
}
if ( bExp == 0 ) {
if ( bSig == 0 ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float32_default_nan;
}
normalizeFloat32Subnormal( bSig, &bExp, &bSig );
}
zSign = ( (int32_t) aSig < 0 );
if ( zSign ) aSig = - aSig;
- return normalizeRoundAndPackFloat32( aSign ^ zSign, bExp, aSig STATUS_VAR );
-
+ return normalizeRoundAndPackFloat32(aSign ^ zSign, bExp, aSig, status);
}
/*----------------------------------------------------------------------------
int shiftcount;
flag signflip, infzero;
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
- c = float32_squash_input_denormal(c STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
+ c = float32_squash_input_denormal(c, status);
aSig = extractFloat32Frac(a);
aExp = extractFloat32Exp(a);
aSign = extractFloat32Sign(a);
if (((aExp == 0xff) && aSig) ||
((bExp == 0xff) && bSig) ||
((cExp == 0xff) && cSig)) {
- return propagateFloat32MulAddNaN(a, b, c, infzero STATUS_VAR);
+ return propagateFloat32MulAddNaN(a, b, c, infzero, status);
}
if (infzero) {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float32_default_nan;
}
if (cExp == 0xff) {
if (pInf && (pSign ^ cSign)) {
/* addition of opposite-signed infinities => InvalidOperation */
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float32_default_nan;
}
/* Otherwise generate an infinity of the same sign */
}
/* Exact zero plus a denorm */
if (STATUS(flush_to_zero)) {
- float_raise(float_flag_output_denormal STATUS_VAR);
+ float_raise(float_flag_output_denormal, status);
return packFloat32(cSign ^ signflip, 0, 0);
}
}
*/
cExp -= 2;
cSig = (cSig | 0x00800000) << 7;
- return roundAndPackFloat32(cSign ^ signflip, cExp, cSig STATUS_VAR);
+ return roundAndPackFloat32(cSign ^ signflip, cExp, cSig, status);
}
return packFloat32(cSign ^ signflip, cExp, cSig);
}
pExp--;
}
return roundAndPackFloat32(zSign, pExp - 1,
- pSig STATUS_VAR);
+ pSig, status);
}
normalizeFloat32Subnormal(cSig, &cExp, &cSig);
}
}
shift64RightJamming(zSig64, 32, &zSig64);
- return roundAndPackFloat32(zSign, zExp, zSig64 STATUS_VAR);
+ return roundAndPackFloat32(zSign, zExp, zSig64, status);
}
int_fast16_t aExp, zExp;
uint32_t aSig, zSig;
uint64_t rem, term;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
aSign = extractFloat32Sign( a );
if ( aExp == 0xFF ) {
- if ( aSig ) return propagateFloat32NaN( a, float32_zero STATUS_VAR );
+ if (aSig) {
+ return propagateFloat32NaN(a, float32_zero, status);
+ }
if ( ! aSign ) return a;
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float32_default_nan;
}
if ( aSign ) {
if ( ( aExp | aSig ) == 0 ) return a;
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float32_default_nan;
}
if ( aExp == 0 ) {
}
shift32RightJamming( zSig, 1, &zSig );
roundAndPack:
- return roundAndPackFloat32( 0, zExp, zSig STATUS_VAR );
+ return roundAndPackFloat32(0, zExp, zSig, status);
}
uint32_t aSig;
float64 r, x, xn;
int i;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
aSign = extractFloat32Sign( a );
if ( aExp == 0xFF) {
- if ( aSig ) return propagateFloat32NaN( a, float32_zero STATUS_VAR );
+ if (aSig) {
+ return propagateFloat32NaN(a, float32_zero, status);
+ }
return (aSign) ? float32_zero : a;
}
if (aExp == 0) {
if (aSig == 0) return float32_one;
}
- float_raise( float_flag_inexact STATUS_VAR);
+ float_raise(float_flag_inexact, status);
/* ******************************* */
/* using float64 for approximation */
/* ******************************* */
- x = float32_to_float64(a STATUS_VAR);
- x = float64_mul(x, float64_ln2 STATUS_VAR);
+ x = float32_to_float64(a, status);
+ x = float64_mul(x, float64_ln2, status);
xn = x;
r = float64_one;
for (i = 0 ; i < 15 ; i++) {
float64 f;
- f = float64_mul(xn, float32_exp2_coefficients[i] STATUS_VAR);
- r = float64_add(r, f STATUS_VAR);
+ f = float64_mul(xn, float32_exp2_coefficients[i], status);
+ r = float64_add(r, f, status);
- xn = float64_mul(xn, x STATUS_VAR);
+ xn = float64_mul(xn, x, status);
}
return float64_to_float32(r, status);
int_fast16_t aExp;
uint32_t aSig, zSig, i;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
aSign = extractFloat32Sign( a );
normalizeFloat32Subnormal( aSig, &aExp, &aSig );
}
if ( aSign ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float32_default_nan;
}
if ( aExp == 0xFF ) {
- if ( aSig ) return propagateFloat32NaN( a, float32_zero STATUS_VAR );
+ if (aSig) {
+ return propagateFloat32NaN(a, float32_zero, status);
+ }
return a;
}
if ( zSign )
zSig = -zSig;
- return normalizeRoundAndPackFloat32( zSign, 0x85, zSig STATUS_VAR );
+ return normalizeRoundAndPackFloat32(zSign, 0x85, zSig, status);
}
/*----------------------------------------------------------------------------
int float32_eq(float32 a, float32 b, float_status *status)
{
uint32_t av, bv;
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
if ( ( ( extractFloat32Exp( a ) == 0xFF ) && extractFloat32Frac( a ) )
|| ( ( extractFloat32Exp( b ) == 0xFF ) && extractFloat32Frac( b ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
av = float32_val(a);
{
flag aSign, bSign;
uint32_t av, bv;
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
if ( ( ( extractFloat32Exp( a ) == 0xFF ) && extractFloat32Frac( a ) )
|| ( ( extractFloat32Exp( b ) == 0xFF ) && extractFloat32Frac( b ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
aSign = extractFloat32Sign( a );
{
flag aSign, bSign;
uint32_t av, bv;
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
if ( ( ( extractFloat32Exp( a ) == 0xFF ) && extractFloat32Frac( a ) )
|| ( ( extractFloat32Exp( b ) == 0xFF ) && extractFloat32Frac( b ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
aSign = extractFloat32Sign( a );
int float32_unordered(float32 a, float32 b, float_status *status)
{
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
if ( ( ( extractFloat32Exp( a ) == 0xFF ) && extractFloat32Frac( a ) )
|| ( ( extractFloat32Exp( b ) == 0xFF ) && extractFloat32Frac( b ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 1;
}
return 0;
int float32_eq_quiet(float32 a, float32 b, float_status *status)
{
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
if ( ( ( extractFloat32Exp( a ) == 0xFF ) && extractFloat32Frac( a ) )
|| ( ( extractFloat32Exp( b ) == 0xFF ) && extractFloat32Frac( b ) )
) {
if ( float32_is_signaling_nan( a ) || float32_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 0;
}
{
flag aSign, bSign;
uint32_t av, bv;
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
if ( ( ( extractFloat32Exp( a ) == 0xFF ) && extractFloat32Frac( a ) )
|| ( ( extractFloat32Exp( b ) == 0xFF ) && extractFloat32Frac( b ) )
) {
if ( float32_is_signaling_nan( a ) || float32_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 0;
}
{
flag aSign, bSign;
uint32_t av, bv;
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
if ( ( ( extractFloat32Exp( a ) == 0xFF ) && extractFloat32Frac( a ) )
|| ( ( extractFloat32Exp( b ) == 0xFF ) && extractFloat32Frac( b ) )
) {
if ( float32_is_signaling_nan( a ) || float32_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 0;
}
int float32_unordered_quiet(float32 a, float32 b, float_status *status)
{
- a = float32_squash_input_denormal(a STATUS_VAR);
- b = float32_squash_input_denormal(b STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
+ b = float32_squash_input_denormal(b, status);
if ( ( ( extractFloat32Exp( a ) == 0xFF ) && extractFloat32Frac( a ) )
|| ( ( extractFloat32Exp( b ) == 0xFF ) && extractFloat32Frac( b ) )
) {
if ( float32_is_signaling_nan( a ) || float32_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 1;
}
flag aSign;
int_fast16_t aExp, shiftCount;
uint64_t aSig;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
if ( aExp ) aSig |= LIT64( 0x0010000000000000 );
shiftCount = 0x42C - aExp;
if ( 0 < shiftCount ) shift64RightJamming( aSig, shiftCount, &aSig );
- return roundAndPackInt32( aSign, aSig STATUS_VAR );
+ return roundAndPackInt32(aSign, aSig, status);
}
int_fast16_t aExp, shiftCount;
uint64_t aSig, savedASig;
int32_t z;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
if ( aSign ) z = - z;
if ( ( z < 0 ) ^ aSign ) {
invalid:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return aSign ? (int32_t) 0x80000000 : 0x7FFFFFFF;
}
if ( ( aSig<<shiftCount ) != savedASig ) {
}
if ( ( (int16_t)z < 0 ) ^ aSign ) {
invalid:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return aSign ? (int32_t) 0xffff8000 : 0x7FFF;
}
if ( ( aSig<<shiftCount ) != savedASig ) {
flag aSign;
int_fast16_t aExp, shiftCount;
uint64_t aSig, aSigExtra;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
shiftCount = 0x433 - aExp;
if ( shiftCount <= 0 ) {
if ( 0x43E < aExp ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
if ( ! aSign
|| ( ( aExp == 0x7FF )
&& ( aSig != LIT64( 0x0010000000000000 ) ) )
else {
shift64ExtraRightJamming( aSig, 0, shiftCount, &aSig, &aSigExtra );
}
- return roundAndPackInt64( aSign, aSig, aSigExtra STATUS_VAR );
+ return roundAndPackInt64(aSign, aSig, aSigExtra, status);
}
int_fast16_t aExp, shiftCount;
uint64_t aSig;
int64 z;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
if ( 0 <= shiftCount ) {
if ( 0x43E <= aExp ) {
if ( float64_val(a) != LIT64( 0xC3E0000000000000 ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
if ( ! aSign
|| ( ( aExp == 0x7FF )
&& ( aSig != LIT64( 0x0010000000000000 ) ) )
int_fast16_t aExp;
uint64_t aSig;
uint32_t zSig;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
aSign = extractFloat64Sign( a );
if ( aExp == 0x7FF ) {
- if ( aSig ) return commonNaNToFloat32( float64ToCommonNaN( a STATUS_VAR ) STATUS_VAR );
+ if (aSig) {
+ return commonNaNToFloat32(float64ToCommonNaN(a, status), status);
+ }
return packFloat32( aSign, 0xFF, 0 );
}
shift64RightJamming( aSig, 22, &aSig );
zSig |= 0x40000000;
aExp -= 0x381;
}
- return roundAndPackFloat32( aSign, aExp, zSig STATUS_VAR );
+ return roundAndPackFloat32(aSign, aExp, zSig, status);
}
if (zExp > maxexp || (zExp == maxexp && rounding_bumps_exp)) {
if (ieee) {
- float_raise(float_flag_overflow | float_flag_inexact STATUS_VAR);
+ float_raise(float_flag_overflow | float_flag_inexact, status);
return packFloat16(zSign, 0x1f, 0);
} else {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return packFloat16(zSign, 0x1f, 0x3ff);
}
}
|| (!rounding_bumps_exp);
}
if (zSig & mask) {
- float_raise(float_flag_inexact STATUS_VAR);
+ float_raise(float_flag_inexact, status);
if (is_tiny) {
- float_raise(float_flag_underflow STATUS_VAR);
+ float_raise(float_flag_underflow, status);
}
}
if (aExp == 0x1f && ieee) {
if (aSig) {
- return commonNaNToFloat32(float16ToCommonNaN(a STATUS_VAR) STATUS_VAR);
+ return commonNaNToFloat32(float16ToCommonNaN(a, status), status);
}
return packFloat32(aSign, 0xff, 0);
}
int_fast16_t aExp;
uint32_t aSig;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
if (aSig) {
/* Input is a NaN */
if (!ieee) {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return packFloat16(aSign, 0, 0);
}
return commonNaNToFloat16(
- float32ToCommonNaN(a STATUS_VAR) STATUS_VAR);
+ float32ToCommonNaN(a, status), status);
}
/* Infinity */
if (!ieee) {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return packFloat16(aSign, 0x1f, 0x3ff);
}
return packFloat16(aSign, 0x1f, 0);
aSig |= 0x00800000;
aExp -= 0x71;
- return roundAndPackFloat16(aSign, aExp, aSig, ieee STATUS_VAR);
+ return roundAndPackFloat16(aSign, aExp, aSig, ieee, status);
}
float64 float16_to_float64(float16 a, flag ieee, float_status *status)
if (aExp == 0x1f && ieee) {
if (aSig) {
return commonNaNToFloat64(
- float16ToCommonNaN(a STATUS_VAR) STATUS_VAR);
+ float16ToCommonNaN(a, status), status);
}
return packFloat64(aSign, 0x7ff, 0);
}
uint64_t aSig;
uint32_t zSig;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aSig = extractFloat64Frac(a);
aExp = extractFloat64Exp(a);
if (aSig) {
/* Input is a NaN */
if (!ieee) {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return packFloat16(aSign, 0, 0);
}
return commonNaNToFloat16(
- float64ToCommonNaN(a STATUS_VAR) STATUS_VAR);
+ float64ToCommonNaN(a, status), status);
}
/* Infinity */
if (!ieee) {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return packFloat16(aSign, 0x1f, 0x3ff);
}
return packFloat16(aSign, 0x1f, 0);
zSig |= 0x00800000;
aExp -= 0x3F1;
- return roundAndPackFloat16(aSign, aExp, zSig, ieee STATUS_VAR);
+ return roundAndPackFloat16(aSign, aExp, zSig, ieee, status);
}
/*----------------------------------------------------------------------------
int_fast16_t aExp;
uint64_t aSig;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
aSign = extractFloat64Sign( a );
if ( aExp == 0x7FF ) {
- if ( aSig ) return commonNaNToFloatx80( float64ToCommonNaN( a STATUS_VAR ) STATUS_VAR );
+ if (aSig) {
+ return commonNaNToFloatx80(float64ToCommonNaN(a, status), status);
+ }
return packFloatx80( aSign, 0x7FFF, LIT64( 0x8000000000000000 ) );
}
if ( aExp == 0 ) {
int_fast16_t aExp;
uint64_t aSig, zSig0, zSig1;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
aSign = extractFloat64Sign( a );
if ( aExp == 0x7FF ) {
- if ( aSig ) return commonNaNToFloat128( float64ToCommonNaN( a STATUS_VAR ) STATUS_VAR );
+ if (aSig) {
+ return commonNaNToFloat128(float64ToCommonNaN(a, status), status);
+ }
return packFloat128( aSign, 0x7FFF, 0, 0 );
}
if ( aExp == 0 ) {
int_fast16_t aExp;
uint64_t lastBitMask, roundBitsMask;
uint64_t z;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aExp = extractFloat64Exp( a );
if ( 0x433 <= aExp ) {
if ( ( aExp == 0x7FF ) && extractFloat64Frac( a ) ) {
- return propagateFloat64NaN( a, a STATUS_VAR );
+ return propagateFloat64NaN(a, a, status);
}
return a;
}
float64 res;
oldmode = STATUS(float_rounding_mode);
STATUS(float_rounding_mode) = float_round_to_zero;
- res = float64_round_to_int(a STATUS_VAR);
+ res = float64_round_to_int(a, status);
STATUS(float_rounding_mode) = oldmode;
return res;
}
bSig <<= 9;
if ( 0 < expDiff ) {
if ( aExp == 0x7FF ) {
- if ( aSig ) return propagateFloat64NaN( a, b STATUS_VAR );
+ if (aSig) {
+ return propagateFloat64NaN(a, b, status);
+ }
return a;
}
if ( bExp == 0 ) {
}
else if ( expDiff < 0 ) {
if ( bExp == 0x7FF ) {
- if ( bSig ) return propagateFloat64NaN( a, b STATUS_VAR );
+ if (bSig) {
+ return propagateFloat64NaN(a, b, status);
+ }
return packFloat64( zSign, 0x7FF, 0 );
}
if ( aExp == 0 ) {
}
else {
if ( aExp == 0x7FF ) {
- if ( aSig | bSig ) return propagateFloat64NaN( a, b STATUS_VAR );
+ if (aSig | bSig) {
+ return propagateFloat64NaN(a, b, status);
+ }
return a;
}
if ( aExp == 0 ) {
if (STATUS(flush_to_zero)) {
if (aSig | bSig) {
- float_raise(float_flag_output_denormal STATUS_VAR);
+ float_raise(float_flag_output_denormal, status);
}
return packFloat64(zSign, 0, 0);
}
++zExp;
}
roundAndPack:
- return roundAndPackFloat64( zSign, zExp, zSig STATUS_VAR );
+ return roundAndPackFloat64(zSign, zExp, zSig, status);
}
if ( 0 < expDiff ) goto aExpBigger;
if ( expDiff < 0 ) goto bExpBigger;
if ( aExp == 0x7FF ) {
- if ( aSig | bSig ) return propagateFloat64NaN( a, b STATUS_VAR );
- float_raise( float_flag_invalid STATUS_VAR);
+ if (aSig | bSig) {
+ return propagateFloat64NaN(a, b, status);
+ }
+ float_raise(float_flag_invalid, status);
return float64_default_nan;
}
if ( aExp == 0 ) {
return packFloat64( STATUS(float_rounding_mode) == float_round_down, 0, 0 );
bExpBigger:
if ( bExp == 0x7FF ) {
- if ( bSig ) return propagateFloat64NaN( a, b STATUS_VAR );
+ if (bSig) {
+ return propagateFloat64NaN(a, b, status);
+ }
return packFloat64( zSign ^ 1, 0x7FF, 0 );
}
if ( aExp == 0 ) {
goto normalizeRoundAndPack;
aExpBigger:
if ( aExp == 0x7FF ) {
- if ( aSig ) return propagateFloat64NaN( a, b STATUS_VAR );
+ if (aSig) {
+ return propagateFloat64NaN(a, b, status);
+ }
return a;
}
if ( bExp == 0 ) {
zExp = aExp;
normalizeRoundAndPack:
--zExp;
- return normalizeRoundAndPackFloat64( zSign, zExp, zSig STATUS_VAR );
+ return normalizeRoundAndPackFloat64(zSign, zExp, zSig, status);
}
float64 float64_add(float64 a, float64 b, float_status *status)
{
flag aSign, bSign;
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
aSign = extractFloat64Sign( a );
bSign = extractFloat64Sign( b );
if ( aSign == bSign ) {
- return addFloat64Sigs( a, b, aSign STATUS_VAR );
+ return addFloat64Sigs(a, b, aSign, status);
}
else {
- return subFloat64Sigs( a, b, aSign STATUS_VAR );
+ return subFloat64Sigs(a, b, aSign, status);
}
}
float64 float64_sub(float64 a, float64 b, float_status *status)
{
flag aSign, bSign;
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
aSign = extractFloat64Sign( a );
bSign = extractFloat64Sign( b );
if ( aSign == bSign ) {
- return subFloat64Sigs( a, b, aSign STATUS_VAR );
+ return subFloat64Sigs(a, b, aSign, status);
}
else {
- return addFloat64Sigs( a, b, aSign STATUS_VAR );
+ return addFloat64Sigs(a, b, aSign, status);
}
}
int_fast16_t aExp, bExp, zExp;
uint64_t aSig, bSig, zSig0, zSig1;
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
zSign = aSign ^ bSign;
if ( aExp == 0x7FF ) {
if ( aSig || ( ( bExp == 0x7FF ) && bSig ) ) {
- return propagateFloat64NaN( a, b STATUS_VAR );
+ return propagateFloat64NaN(a, b, status);
}
if ( ( bExp | bSig ) == 0 ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float64_default_nan;
}
return packFloat64( zSign, 0x7FF, 0 );
}
if ( bExp == 0x7FF ) {
- if ( bSig ) return propagateFloat64NaN( a, b STATUS_VAR );
+ if (bSig) {
+ return propagateFloat64NaN(a, b, status);
+ }
if ( ( aExp | aSig ) == 0 ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float64_default_nan;
}
return packFloat64( zSign, 0x7FF, 0 );
zSig0 <<= 1;
--zExp;
}
- return roundAndPackFloat64( zSign, zExp, zSig0 STATUS_VAR );
+ return roundAndPackFloat64(zSign, zExp, zSig0, status);
}
uint64_t aSig, bSig, zSig;
uint64_t rem0, rem1;
uint64_t term0, term1;
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
bSign = extractFloat64Sign( b );
zSign = aSign ^ bSign;
if ( aExp == 0x7FF ) {
- if ( aSig ) return propagateFloat64NaN( a, b STATUS_VAR );
+ if (aSig) {
+ return propagateFloat64NaN(a, b, status);
+ }
if ( bExp == 0x7FF ) {
- if ( bSig ) return propagateFloat64NaN( a, b STATUS_VAR );
- float_raise( float_flag_invalid STATUS_VAR);
+ if (bSig) {
+ return propagateFloat64NaN(a, b, status);
+ }
+ float_raise(float_flag_invalid, status);
return float64_default_nan;
}
return packFloat64( zSign, 0x7FF, 0 );
}
if ( bExp == 0x7FF ) {
- if ( bSig ) return propagateFloat64NaN( a, b STATUS_VAR );
+ if (bSig) {
+ return propagateFloat64NaN(a, b, status);
+ }
return packFloat64( zSign, 0, 0 );
}
if ( bExp == 0 ) {
if ( bSig == 0 ) {
if ( ( aExp | aSig ) == 0 ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float64_default_nan;
}
- float_raise( float_flag_divbyzero STATUS_VAR);
+ float_raise(float_flag_divbyzero, status);
return packFloat64( zSign, 0x7FF, 0 );
}
normalizeFloat64Subnormal( bSig, &bExp, &bSig );
}
zSig |= ( rem1 != 0 );
}
- return roundAndPackFloat64( zSign, zExp, zSig STATUS_VAR );
+ return roundAndPackFloat64(zSign, zExp, zSig, status);
}
uint64_t q, alternateASig;
int64_t sigMean;
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
aSign = extractFloat64Sign( a );
bExp = extractFloat64Exp( b );
if ( aExp == 0x7FF ) {
if ( aSig || ( ( bExp == 0x7FF ) && bSig ) ) {
- return propagateFloat64NaN( a, b STATUS_VAR );
+ return propagateFloat64NaN(a, b, status);
}
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float64_default_nan;
}
if ( bExp == 0x7FF ) {
- if ( bSig ) return propagateFloat64NaN( a, b STATUS_VAR );
+ if (bSig) {
+ return propagateFloat64NaN(a, b, status);
+ }
return a;
}
if ( bExp == 0 ) {
if ( bSig == 0 ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float64_default_nan;
}
normalizeFloat64Subnormal( bSig, &bExp, &bSig );
}
zSign = ( (int64_t) aSig < 0 );
if ( zSign ) aSig = - aSig;
- return normalizeRoundAndPackFloat64( aSign ^ zSign, bExp, aSig STATUS_VAR );
+ return normalizeRoundAndPackFloat64(aSign ^ zSign, bExp, aSig, status);
}
int shiftcount;
flag signflip, infzero;
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
- c = float64_squash_input_denormal(c STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
+ c = float64_squash_input_denormal(c, status);
aSig = extractFloat64Frac(a);
aExp = extractFloat64Exp(a);
aSign = extractFloat64Sign(a);
if (((aExp == 0x7ff) && aSig) ||
((bExp == 0x7ff) && bSig) ||
((cExp == 0x7ff) && cSig)) {
- return propagateFloat64MulAddNaN(a, b, c, infzero STATUS_VAR);
+ return propagateFloat64MulAddNaN(a, b, c, infzero, status);
}
if (infzero) {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float64_default_nan;
}
if (cExp == 0x7ff) {
if (pInf && (pSign ^ cSign)) {
/* addition of opposite-signed infinities => InvalidOperation */
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float64_default_nan;
}
/* Otherwise generate an infinity of the same sign */
}
/* Exact zero plus a denorm */
if (STATUS(flush_to_zero)) {
- float_raise(float_flag_output_denormal STATUS_VAR);
+ float_raise(float_flag_output_denormal, status);
return packFloat64(cSign ^ signflip, 0, 0);
}
}
*/
cExp -= 2;
cSig = (cSig | 0x0010000000000000ULL) << 10;
- return roundAndPackFloat64(cSign ^ signflip, cExp, cSig STATUS_VAR);
+ return roundAndPackFloat64(cSign ^ signflip, cExp, cSig, status);
}
return packFloat64(cSign ^ signflip, cExp, cSig);
}
pExp--;
}
return roundAndPackFloat64(zSign, pExp - 1,
- pSig1 STATUS_VAR);
+ pSig1, status);
}
normalizeFloat64Subnormal(cSig, &cExp, &cSig);
}
if (flags & float_muladd_halve_result) {
zExp--;
}
- return roundAndPackFloat64(zSign, zExp, zSig1 STATUS_VAR);
+ return roundAndPackFloat64(zSign, zExp, zSig1, status);
} else {
/* Subtraction */
if (expDiff > 0) {
if (flags & float_muladd_halve_result) {
zExp--;
}
- return roundAndPackFloat64(zSign, zExp, zSig0 STATUS_VAR);
+ return roundAndPackFloat64(zSign, zExp, zSig0, status);
}
}
int_fast16_t aExp, zExp;
uint64_t aSig, zSig, doubleZSig;
uint64_t rem0, rem1, term0, term1;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
aSign = extractFloat64Sign( a );
if ( aExp == 0x7FF ) {
- if ( aSig ) return propagateFloat64NaN( a, a STATUS_VAR );
+ if (aSig) {
+ return propagateFloat64NaN(a, a, status);
+ }
if ( ! aSign ) return a;
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float64_default_nan;
}
if ( aSign ) {
if ( ( aExp | aSig ) == 0 ) return a;
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float64_default_nan;
}
if ( aExp == 0 ) {
}
zSig |= ( ( rem0 | rem1 ) != 0 );
}
- return roundAndPackFloat64( 0, zExp, zSig STATUS_VAR );
+ return roundAndPackFloat64(0, zExp, zSig, status);
}
flag aSign, zSign;
int_fast16_t aExp;
uint64_t aSig, aSig0, aSig1, zSig, i;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
normalizeFloat64Subnormal( aSig, &aExp, &aSig );
}
if ( aSign ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return float64_default_nan;
}
if ( aExp == 0x7FF ) {
- if ( aSig ) return propagateFloat64NaN( a, float64_zero STATUS_VAR );
+ if (aSig) {
+ return propagateFloat64NaN(a, float64_zero, status);
+ }
return a;
}
if ( zSign )
zSig = -zSig;
- return normalizeRoundAndPackFloat64( zSign, 0x408, zSig STATUS_VAR );
+ return normalizeRoundAndPackFloat64(zSign, 0x408, zSig, status);
}
/*----------------------------------------------------------------------------
int float64_eq(float64 a, float64 b, float_status *status)
{
uint64_t av, bv;
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
if ( ( ( extractFloat64Exp( a ) == 0x7FF ) && extractFloat64Frac( a ) )
|| ( ( extractFloat64Exp( b ) == 0x7FF ) && extractFloat64Frac( b ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
av = float64_val(a);
{
flag aSign, bSign;
uint64_t av, bv;
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
if ( ( ( extractFloat64Exp( a ) == 0x7FF ) && extractFloat64Frac( a ) )
|| ( ( extractFloat64Exp( b ) == 0x7FF ) && extractFloat64Frac( b ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
aSign = extractFloat64Sign( a );
flag aSign, bSign;
uint64_t av, bv;
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
if ( ( ( extractFloat64Exp( a ) == 0x7FF ) && extractFloat64Frac( a ) )
|| ( ( extractFloat64Exp( b ) == 0x7FF ) && extractFloat64Frac( b ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
aSign = extractFloat64Sign( a );
int float64_unordered(float64 a, float64 b, float_status *status)
{
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
if ( ( ( extractFloat64Exp( a ) == 0x7FF ) && extractFloat64Frac( a ) )
|| ( ( extractFloat64Exp( b ) == 0x7FF ) && extractFloat64Frac( b ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 1;
}
return 0;
int float64_eq_quiet(float64 a, float64 b, float_status *status)
{
uint64_t av, bv;
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
if ( ( ( extractFloat64Exp( a ) == 0x7FF ) && extractFloat64Frac( a ) )
|| ( ( extractFloat64Exp( b ) == 0x7FF ) && extractFloat64Frac( b ) )
) {
if ( float64_is_signaling_nan( a ) || float64_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 0;
}
{
flag aSign, bSign;
uint64_t av, bv;
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
if ( ( ( extractFloat64Exp( a ) == 0x7FF ) && extractFloat64Frac( a ) )
|| ( ( extractFloat64Exp( b ) == 0x7FF ) && extractFloat64Frac( b ) )
) {
if ( float64_is_signaling_nan( a ) || float64_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 0;
}
{
flag aSign, bSign;
uint64_t av, bv;
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
if ( ( ( extractFloat64Exp( a ) == 0x7FF ) && extractFloat64Frac( a ) )
|| ( ( extractFloat64Exp( b ) == 0x7FF ) && extractFloat64Frac( b ) )
) {
if ( float64_is_signaling_nan( a ) || float64_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 0;
}
int float64_unordered_quiet(float64 a, float64 b, float_status *status)
{
- a = float64_squash_input_denormal(a STATUS_VAR);
- b = float64_squash_input_denormal(b STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
+ b = float64_squash_input_denormal(b, status);
if ( ( ( extractFloat64Exp( a ) == 0x7FF ) && extractFloat64Frac( a ) )
|| ( ( extractFloat64Exp( b ) == 0x7FF ) && extractFloat64Frac( b ) )
) {
if ( float64_is_signaling_nan( a ) || float64_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 1;
}
shiftCount = 0x4037 - aExp;
if ( shiftCount <= 0 ) shiftCount = 1;
shift64RightJamming( aSig, shiftCount, &aSig );
- return roundAndPackInt32( aSign, aSig STATUS_VAR );
+ return roundAndPackInt32(aSign, aSig, status);
}
if ( aSign ) z = - z;
if ( ( z < 0 ) ^ aSign ) {
invalid:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return aSign ? (int32_t) 0x80000000 : 0x7FFFFFFF;
}
if ( ( aSig<<shiftCount ) != savedASig ) {
shiftCount = 0x403E - aExp;
if ( shiftCount <= 0 ) {
if ( shiftCount ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
if ( ! aSign
|| ( ( aExp == 0x7FFF )
&& ( aSig != LIT64( 0x8000000000000000 ) ) )
else {
shift64ExtraRightJamming( aSig, 0, shiftCount, &aSig, &aSigExtra );
}
- return roundAndPackInt64( aSign, aSig, aSigExtra STATUS_VAR );
+ return roundAndPackInt64(aSign, aSig, aSigExtra, status);
}
if ( 0 <= shiftCount ) {
aSig &= LIT64( 0x7FFFFFFFFFFFFFFF );
if ( ( a.high != 0xC03E ) || aSig ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
if ( ! aSign || ( ( aExp == 0x7FFF ) && aSig ) ) {
return LIT64( 0x7FFFFFFFFFFFFFFF );
}
aSign = extractFloatx80Sign( a );
if ( aExp == 0x7FFF ) {
if ( (uint64_t) ( aSig<<1 ) ) {
- return commonNaNToFloat32( floatx80ToCommonNaN( a STATUS_VAR ) STATUS_VAR );
+ return commonNaNToFloat32(floatx80ToCommonNaN(a, status), status);
}
return packFloat32( aSign, 0xFF, 0 );
}
shift64RightJamming( aSig, 33, &aSig );
if ( aExp || aSig ) aExp -= 0x3F81;
- return roundAndPackFloat32( aSign, aExp, aSig STATUS_VAR );
+ return roundAndPackFloat32(aSign, aExp, aSig, status);
}
aSign = extractFloatx80Sign( a );
if ( aExp == 0x7FFF ) {
if ( (uint64_t) ( aSig<<1 ) ) {
- return commonNaNToFloat64( floatx80ToCommonNaN( a STATUS_VAR ) STATUS_VAR );
+ return commonNaNToFloat64(floatx80ToCommonNaN(a, status), status);
}
return packFloat64( aSign, 0x7FF, 0 );
}
shift64RightJamming( aSig, 1, &zSig );
if ( aExp || aSig ) aExp -= 0x3C01;
- return roundAndPackFloat64( aSign, aExp, zSig STATUS_VAR );
+ return roundAndPackFloat64(aSign, aExp, zSig, status);
}
aExp = extractFloatx80Exp( a );
aSign = extractFloatx80Sign( a );
if ( ( aExp == 0x7FFF ) && (uint64_t) ( aSig<<1 ) ) {
- return commonNaNToFloat128( floatx80ToCommonNaN( a STATUS_VAR ) STATUS_VAR );
+ return commonNaNToFloat128(floatx80ToCommonNaN(a, status), status);
}
shift128Right( aSig<<1, 0, 16, &zSig0, &zSig1 );
return packFloat128( aSign, aExp, zSig0, zSig1 );
aExp = extractFloatx80Exp( a );
if ( 0x403E <= aExp ) {
if ( ( aExp == 0x7FFF ) && (uint64_t) ( extractFloatx80Frac( a )<<1 ) ) {
- return propagateFloatx80NaN( a, a STATUS_VAR );
+ return propagateFloatx80NaN(a, a, status);
}
return a;
}
expDiff = aExp - bExp;
if ( 0 < expDiff ) {
if ( aExp == 0x7FFF ) {
- if ( (uint64_t) ( aSig<<1 ) ) return propagateFloatx80NaN( a, b STATUS_VAR );
+ if ((uint64_t)(aSig << 1)) {
+ return propagateFloatx80NaN(a, b, status);
+ }
return a;
}
if ( bExp == 0 ) --expDiff;
}
else if ( expDiff < 0 ) {
if ( bExp == 0x7FFF ) {
- if ( (uint64_t) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b STATUS_VAR );
+ if ((uint64_t)(bSig << 1)) {
+ return propagateFloatx80NaN(a, b, status);
+ }
return packFloatx80( zSign, 0x7FFF, LIT64( 0x8000000000000000 ) );
}
if ( aExp == 0 ) ++expDiff;
else {
if ( aExp == 0x7FFF ) {
if ( (uint64_t) ( ( aSig | bSig )<<1 ) ) {
- return propagateFloatx80NaN( a, b STATUS_VAR );
+ return propagateFloatx80NaN(a, b, status);
}
return a;
}
zSig0 |= LIT64( 0x8000000000000000 );
++zExp;
roundAndPack:
- return
- roundAndPackFloatx80(
- STATUS(floatx80_rounding_precision), zSign, zExp, zSig0, zSig1 STATUS_VAR );
-
+ return roundAndPackFloatx80(STATUS(floatx80_rounding_precision),
+ zSign, zExp, zSig0, zSig1, status);
}
/*----------------------------------------------------------------------------
if ( expDiff < 0 ) goto bExpBigger;
if ( aExp == 0x7FFF ) {
if ( (uint64_t) ( ( aSig | bSig )<<1 ) ) {
- return propagateFloatx80NaN( a, b STATUS_VAR );
+ return propagateFloatx80NaN(a, b, status);
}
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
z.low = floatx80_default_nan_low;
z.high = floatx80_default_nan_high;
return z;
return packFloatx80( STATUS(float_rounding_mode) == float_round_down, 0, 0 );
bExpBigger:
if ( bExp == 0x7FFF ) {
- if ( (uint64_t) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b STATUS_VAR );
+ if ((uint64_t)(bSig << 1)) {
+ return propagateFloatx80NaN(a, b, status);
+ }
return packFloatx80( zSign ^ 1, 0x7FFF, LIT64( 0x8000000000000000 ) );
}
if ( aExp == 0 ) ++expDiff;
goto normalizeRoundAndPack;
aExpBigger:
if ( aExp == 0x7FFF ) {
- if ( (uint64_t) ( aSig<<1 ) ) return propagateFloatx80NaN( a, b STATUS_VAR );
+ if ((uint64_t)(aSig << 1)) {
+ return propagateFloatx80NaN(a, b, status);
+ }
return a;
}
if ( bExp == 0 ) --expDiff;
sub128( aSig, 0, bSig, zSig1, &zSig0, &zSig1 );
zExp = aExp;
normalizeRoundAndPack:
- return
- normalizeRoundAndPackFloatx80(
- STATUS(floatx80_rounding_precision), zSign, zExp, zSig0, zSig1 STATUS_VAR );
-
+ return normalizeRoundAndPackFloatx80(STATUS(floatx80_rounding_precision),
+ zSign, zExp, zSig0, zSig1, status);
}
/*----------------------------------------------------------------------------
aSign = extractFloatx80Sign( a );
bSign = extractFloatx80Sign( b );
if ( aSign == bSign ) {
- return addFloatx80Sigs( a, b, aSign STATUS_VAR );
+ return addFloatx80Sigs(a, b, aSign, status);
}
else {
- return subFloatx80Sigs( a, b, aSign STATUS_VAR );
+ return subFloatx80Sigs(a, b, aSign, status);
}
}
aSign = extractFloatx80Sign( a );
bSign = extractFloatx80Sign( b );
if ( aSign == bSign ) {
- return subFloatx80Sigs( a, b, aSign STATUS_VAR );
+ return subFloatx80Sigs(a, b, aSign, status);
}
else {
- return addFloatx80Sigs( a, b, aSign STATUS_VAR );
+ return addFloatx80Sigs(a, b, aSign, status);
}
}
if ( aExp == 0x7FFF ) {
if ( (uint64_t) ( aSig<<1 )
|| ( ( bExp == 0x7FFF ) && (uint64_t) ( bSig<<1 ) ) ) {
- return propagateFloatx80NaN( a, b STATUS_VAR );
+ return propagateFloatx80NaN(a, b, status);
}
if ( ( bExp | bSig ) == 0 ) goto invalid;
return packFloatx80( zSign, 0x7FFF, LIT64( 0x8000000000000000 ) );
}
if ( bExp == 0x7FFF ) {
- if ( (uint64_t) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b STATUS_VAR );
+ if ((uint64_t)(bSig << 1)) {
+ return propagateFloatx80NaN(a, b, status);
+ }
if ( ( aExp | aSig ) == 0 ) {
invalid:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
z.low = floatx80_default_nan_low;
z.high = floatx80_default_nan_high;
return z;
shortShift128Left( zSig0, zSig1, 1, &zSig0, &zSig1 );
--zExp;
}
- return
- roundAndPackFloatx80(
- STATUS(floatx80_rounding_precision), zSign, zExp, zSig0, zSig1 STATUS_VAR );
-
+ return roundAndPackFloatx80(STATUS(floatx80_rounding_precision),
+ zSign, zExp, zSig0, zSig1, status);
}
/*----------------------------------------------------------------------------
bSign = extractFloatx80Sign( b );
zSign = aSign ^ bSign;
if ( aExp == 0x7FFF ) {
- if ( (uint64_t) ( aSig<<1 ) ) return propagateFloatx80NaN( a, b STATUS_VAR );
+ if ((uint64_t)(aSig << 1)) {
+ return propagateFloatx80NaN(a, b, status);
+ }
if ( bExp == 0x7FFF ) {
- if ( (uint64_t) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b STATUS_VAR );
+ if ((uint64_t)(bSig << 1)) {
+ return propagateFloatx80NaN(a, b, status);
+ }
goto invalid;
}
return packFloatx80( zSign, 0x7FFF, LIT64( 0x8000000000000000 ) );
}
if ( bExp == 0x7FFF ) {
- if ( (uint64_t) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b STATUS_VAR );
+ if ((uint64_t)(bSig << 1)) {
+ return propagateFloatx80NaN(a, b, status);
+ }
return packFloatx80( zSign, 0, 0 );
}
if ( bExp == 0 ) {
if ( bSig == 0 ) {
if ( ( aExp | aSig ) == 0 ) {
invalid:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
z.low = floatx80_default_nan_low;
z.high = floatx80_default_nan_high;
return z;
}
- float_raise( float_flag_divbyzero STATUS_VAR);
+ float_raise(float_flag_divbyzero, status);
return packFloatx80( zSign, 0x7FFF, LIT64( 0x8000000000000000 ) );
}
normalizeFloatx80Subnormal( bSig, &bExp, &bSig );
}
zSig1 |= ( ( rem1 | rem2 ) != 0 );
}
- return
- roundAndPackFloatx80(
- STATUS(floatx80_rounding_precision), zSign, zExp, zSig0, zSig1 STATUS_VAR );
-
+ return roundAndPackFloatx80(STATUS(floatx80_rounding_precision),
+ zSign, zExp, zSig0, zSig1, status);
}
/*----------------------------------------------------------------------------
if ( aExp == 0x7FFF ) {
if ( (uint64_t) ( aSig0<<1 )
|| ( ( bExp == 0x7FFF ) && (uint64_t) ( bSig<<1 ) ) ) {
- return propagateFloatx80NaN( a, b STATUS_VAR );
+ return propagateFloatx80NaN(a, b, status);
}
goto invalid;
}
if ( bExp == 0x7FFF ) {
- if ( (uint64_t) ( bSig<<1 ) ) return propagateFloatx80NaN( a, b STATUS_VAR );
+ if ((uint64_t)(bSig << 1)) {
+ return propagateFloatx80NaN(a, b, status);
+ }
return a;
}
if ( bExp == 0 ) {
if ( bSig == 0 ) {
invalid:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
z.low = floatx80_default_nan_low;
z.high = floatx80_default_nan_high;
return z;
}
return
normalizeRoundAndPackFloatx80(
- 80, zSign, bExp + expDiff, aSig0, aSig1 STATUS_VAR );
+ 80, zSign, bExp + expDiff, aSig0, aSig1, status);
}
aExp = extractFloatx80Exp( a );
aSign = extractFloatx80Sign( a );
if ( aExp == 0x7FFF ) {
- if ( (uint64_t) ( aSig0<<1 ) ) return propagateFloatx80NaN( a, a STATUS_VAR );
+ if ((uint64_t)(aSig0 << 1)) {
+ return propagateFloatx80NaN(a, a, status);
+ }
if ( ! aSign ) return a;
goto invalid;
}
if ( aSign ) {
if ( ( aExp | aSig0 ) == 0 ) return a;
invalid:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
z.low = floatx80_default_nan_low;
z.high = floatx80_default_nan_high;
return z;
zSig0 |= doubleZSig0;
return
roundAndPackFloatx80(
- STATUS(floatx80_rounding_precision), 0, zExp, zSig0, zSig1 STATUS_VAR );
+ STATUS(floatx80_rounding_precision), 0, zExp, zSig0, zSig1, status);
}
|| ( ( extractFloatx80Exp( b ) == 0x7FFF )
&& (uint64_t) ( extractFloatx80Frac( b )<<1 ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
return
|| ( ( extractFloatx80Exp( b ) == 0x7FFF )
&& (uint64_t) ( extractFloatx80Frac( b )<<1 ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
aSign = extractFloatx80Sign( a );
|| ( ( extractFloatx80Exp( b ) == 0x7FFF )
&& (uint64_t) ( extractFloatx80Frac( b )<<1 ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
aSign = extractFloatx80Sign( a );
|| ( ( extractFloatx80Exp( b ) == 0x7FFF )
&& (uint64_t) ( extractFloatx80Frac( b )<<1 ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 1;
}
return 0;
) {
if ( floatx80_is_signaling_nan( a )
|| floatx80_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 0;
}
) {
if ( floatx80_is_signaling_nan( a )
|| floatx80_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 0;
}
) {
if ( floatx80_is_signaling_nan( a )
|| floatx80_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 0;
}
) {
if ( floatx80_is_signaling_nan( a )
|| floatx80_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 1;
}
aSig0 |= ( aSig1 != 0 );
shiftCount = 0x4028 - aExp;
if ( 0 < shiftCount ) shift64RightJamming( aSig0, shiftCount, &aSig0 );
- return roundAndPackInt32( aSign, aSig0 STATUS_VAR );
+ return roundAndPackInt32(aSign, aSig0, status);
}
if ( aSign ) z = - z;
if ( ( z < 0 ) ^ aSign ) {
invalid:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return aSign ? (int32_t) 0x80000000 : 0x7FFFFFFF;
}
if ( ( aSig0<<shiftCount ) != savedASig ) {
shiftCount = 0x402F - aExp;
if ( shiftCount <= 0 ) {
if ( 0x403E < aExp ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
if ( ! aSign
|| ( ( aExp == 0x7FFF )
&& ( aSig1 || ( aSig0 != LIT64( 0x0001000000000000 ) ) )
else {
shift64ExtraRightJamming( aSig0, aSig1, shiftCount, &aSig0, &aSig1 );
}
- return roundAndPackInt64( aSign, aSig0, aSig1 STATUS_VAR );
+ return roundAndPackInt64(aSign, aSig0, aSig1, status);
}
if ( aSig1 ) STATUS(float_exception_flags) |= float_flag_inexact;
}
else {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
if ( ! aSign || ( ( aExp == 0x7FFF ) && ( aSig0 | aSig1 ) ) ) {
return LIT64( 0x7FFFFFFFFFFFFFFF );
}
aSign = extractFloat128Sign( a );
if ( aExp == 0x7FFF ) {
if ( aSig0 | aSig1 ) {
- return commonNaNToFloat32( float128ToCommonNaN( a STATUS_VAR ) STATUS_VAR );
+ return commonNaNToFloat32(float128ToCommonNaN(a, status), status);
}
return packFloat32( aSign, 0xFF, 0 );
}
zSig |= 0x40000000;
aExp -= 0x3F81;
}
- return roundAndPackFloat32( aSign, aExp, zSig STATUS_VAR );
+ return roundAndPackFloat32(aSign, aExp, zSig, status);
}
aSign = extractFloat128Sign( a );
if ( aExp == 0x7FFF ) {
if ( aSig0 | aSig1 ) {
- return commonNaNToFloat64( float128ToCommonNaN( a STATUS_VAR ) STATUS_VAR );
+ return commonNaNToFloat64(float128ToCommonNaN(a, status), status);
}
return packFloat64( aSign, 0x7FF, 0 );
}
aSig0 |= LIT64( 0x4000000000000000 );
aExp -= 0x3C01;
}
- return roundAndPackFloat64( aSign, aExp, aSig0 STATUS_VAR );
+ return roundAndPackFloat64(aSign, aExp, aSig0, status);
}
aSign = extractFloat128Sign( a );
if ( aExp == 0x7FFF ) {
if ( aSig0 | aSig1 ) {
- return commonNaNToFloatx80( float128ToCommonNaN( a STATUS_VAR ) STATUS_VAR );
+ return commonNaNToFloatx80(float128ToCommonNaN(a, status), status);
}
return packFloatx80( aSign, 0x7FFF, LIT64( 0x8000000000000000 ) );
}
aSig0 |= LIT64( 0x0001000000000000 );
}
shortShift128Left( aSig0, aSig1, 15, &aSig0, &aSig1 );
- return roundAndPackFloatx80( 80, aSign, aExp, aSig0, aSig1 STATUS_VAR );
+ return roundAndPackFloatx80(80, aSign, aExp, aSig0, aSig1, status);
}
if ( ( aExp == 0x7FFF )
&& ( extractFloat128Frac0( a ) | extractFloat128Frac1( a ) )
) {
- return propagateFloat128NaN( a, a STATUS_VAR );
+ return propagateFloat128NaN(a, a, status);
}
return a;
}
expDiff = aExp - bExp;
if ( 0 < expDiff ) {
if ( aExp == 0x7FFF ) {
- if ( aSig0 | aSig1 ) return propagateFloat128NaN( a, b STATUS_VAR );
+ if (aSig0 | aSig1) {
+ return propagateFloat128NaN(a, b, status);
+ }
return a;
}
if ( bExp == 0 ) {
}
else if ( expDiff < 0 ) {
if ( bExp == 0x7FFF ) {
- if ( bSig0 | bSig1 ) return propagateFloat128NaN( a, b STATUS_VAR );
+ if (bSig0 | bSig1) {
+ return propagateFloat128NaN(a, b, status);
+ }
return packFloat128( zSign, 0x7FFF, 0, 0 );
}
if ( aExp == 0 ) {
else {
if ( aExp == 0x7FFF ) {
if ( aSig0 | aSig1 | bSig0 | bSig1 ) {
- return propagateFloat128NaN( a, b STATUS_VAR );
+ return propagateFloat128NaN(a, b, status);
}
return a;
}
if ( aExp == 0 ) {
if (STATUS(flush_to_zero)) {
if (zSig0 | zSig1) {
- float_raise(float_flag_output_denormal STATUS_VAR);
+ float_raise(float_flag_output_denormal, status);
}
return packFloat128(zSign, 0, 0, 0);
}
shift128ExtraRightJamming(
zSig0, zSig1, zSig2, 1, &zSig0, &zSig1, &zSig2 );
roundAndPack:
- return roundAndPackFloat128( zSign, zExp, zSig0, zSig1, zSig2 STATUS_VAR );
+ return roundAndPackFloat128(zSign, zExp, zSig0, zSig1, zSig2, status);
}
if ( expDiff < 0 ) goto bExpBigger;
if ( aExp == 0x7FFF ) {
if ( aSig0 | aSig1 | bSig0 | bSig1 ) {
- return propagateFloat128NaN( a, b STATUS_VAR );
+ return propagateFloat128NaN(a, b, status);
}
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
z.low = float128_default_nan_low;
z.high = float128_default_nan_high;
return z;
return packFloat128( STATUS(float_rounding_mode) == float_round_down, 0, 0, 0 );
bExpBigger:
if ( bExp == 0x7FFF ) {
- if ( bSig0 | bSig1 ) return propagateFloat128NaN( a, b STATUS_VAR );
+ if (bSig0 | bSig1) {
+ return propagateFloat128NaN(a, b, status);
+ }
return packFloat128( zSign ^ 1, 0x7FFF, 0, 0 );
}
if ( aExp == 0 ) {
goto normalizeRoundAndPack;
aExpBigger:
if ( aExp == 0x7FFF ) {
- if ( aSig0 | aSig1 ) return propagateFloat128NaN( a, b STATUS_VAR );
+ if (aSig0 | aSig1) {
+ return propagateFloat128NaN(a, b, status);
+ }
return a;
}
if ( bExp == 0 ) {
zExp = aExp;
normalizeRoundAndPack:
--zExp;
- return normalizeRoundAndPackFloat128( zSign, zExp - 14, zSig0, zSig1 STATUS_VAR );
+ return normalizeRoundAndPackFloat128(zSign, zExp - 14, zSig0, zSig1,
+ status);
}
aSign = extractFloat128Sign( a );
bSign = extractFloat128Sign( b );
if ( aSign == bSign ) {
- return addFloat128Sigs( a, b, aSign STATUS_VAR );
+ return addFloat128Sigs(a, b, aSign, status);
}
else {
- return subFloat128Sigs( a, b, aSign STATUS_VAR );
+ return subFloat128Sigs(a, b, aSign, status);
}
}
aSign = extractFloat128Sign( a );
bSign = extractFloat128Sign( b );
if ( aSign == bSign ) {
- return subFloat128Sigs( a, b, aSign STATUS_VAR );
+ return subFloat128Sigs(a, b, aSign, status);
}
else {
- return addFloat128Sigs( a, b, aSign STATUS_VAR );
+ return addFloat128Sigs(a, b, aSign, status);
}
}
if ( aExp == 0x7FFF ) {
if ( ( aSig0 | aSig1 )
|| ( ( bExp == 0x7FFF ) && ( bSig0 | bSig1 ) ) ) {
- return propagateFloat128NaN( a, b STATUS_VAR );
+ return propagateFloat128NaN(a, b, status);
}
if ( ( bExp | bSig0 | bSig1 ) == 0 ) goto invalid;
return packFloat128( zSign, 0x7FFF, 0, 0 );
}
if ( bExp == 0x7FFF ) {
- if ( bSig0 | bSig1 ) return propagateFloat128NaN( a, b STATUS_VAR );
+ if (bSig0 | bSig1) {
+ return propagateFloat128NaN(a, b, status);
+ }
if ( ( aExp | aSig0 | aSig1 ) == 0 ) {
invalid:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
z.low = float128_default_nan_low;
z.high = float128_default_nan_high;
return z;
zSig0, zSig1, zSig2, 1, &zSig0, &zSig1, &zSig2 );
++zExp;
}
- return roundAndPackFloat128( zSign, zExp, zSig0, zSig1, zSig2 STATUS_VAR );
+ return roundAndPackFloat128(zSign, zExp, zSig0, zSig1, zSig2, status);
}
bSign = extractFloat128Sign( b );
zSign = aSign ^ bSign;
if ( aExp == 0x7FFF ) {
- if ( aSig0 | aSig1 ) return propagateFloat128NaN( a, b STATUS_VAR );
+ if (aSig0 | aSig1) {
+ return propagateFloat128NaN(a, b, status);
+ }
if ( bExp == 0x7FFF ) {
- if ( bSig0 | bSig1 ) return propagateFloat128NaN( a, b STATUS_VAR );
+ if (bSig0 | bSig1) {
+ return propagateFloat128NaN(a, b, status);
+ }
goto invalid;
}
return packFloat128( zSign, 0x7FFF, 0, 0 );
}
if ( bExp == 0x7FFF ) {
- if ( bSig0 | bSig1 ) return propagateFloat128NaN( a, b STATUS_VAR );
+ if (bSig0 | bSig1) {
+ return propagateFloat128NaN(a, b, status);
+ }
return packFloat128( zSign, 0, 0, 0 );
}
if ( bExp == 0 ) {
if ( ( bSig0 | bSig1 ) == 0 ) {
if ( ( aExp | aSig0 | aSig1 ) == 0 ) {
invalid:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
z.low = float128_default_nan_low;
z.high = float128_default_nan_high;
return z;
}
- float_raise( float_flag_divbyzero STATUS_VAR);
+ float_raise(float_flag_divbyzero, status);
return packFloat128( zSign, 0x7FFF, 0, 0 );
}
normalizeFloat128Subnormal( bSig0, bSig1, &bExp, &bSig0, &bSig1 );
zSig1 |= ( ( rem1 | rem2 | rem3 ) != 0 );
}
shift128ExtraRightJamming( zSig0, zSig1, 0, 15, &zSig0, &zSig1, &zSig2 );
- return roundAndPackFloat128( zSign, zExp, zSig0, zSig1, zSig2 STATUS_VAR );
+ return roundAndPackFloat128(zSign, zExp, zSig0, zSig1, zSig2, status);
}
if ( aExp == 0x7FFF ) {
if ( ( aSig0 | aSig1 )
|| ( ( bExp == 0x7FFF ) && ( bSig0 | bSig1 ) ) ) {
- return propagateFloat128NaN( a, b STATUS_VAR );
+ return propagateFloat128NaN(a, b, status);
}
goto invalid;
}
if ( bExp == 0x7FFF ) {
- if ( bSig0 | bSig1 ) return propagateFloat128NaN( a, b STATUS_VAR );
+ if (bSig0 | bSig1) {
+ return propagateFloat128NaN(a, b, status);
+ }
return a;
}
if ( bExp == 0 ) {
if ( ( bSig0 | bSig1 ) == 0 ) {
invalid:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
z.low = float128_default_nan_low;
z.high = float128_default_nan_high;
return z;
}
zSign = ( (int64_t) aSig0 < 0 );
if ( zSign ) sub128( 0, 0, aSig0, aSig1, &aSig0, &aSig1 );
- return
- normalizeRoundAndPackFloat128( aSign ^ zSign, bExp - 4, aSig0, aSig1 STATUS_VAR );
-
+ return normalizeRoundAndPackFloat128(aSign ^ zSign, bExp - 4, aSig0, aSig1,
+ status);
}
/*----------------------------------------------------------------------------
aExp = extractFloat128Exp( a );
aSign = extractFloat128Sign( a );
if ( aExp == 0x7FFF ) {
- if ( aSig0 | aSig1 ) return propagateFloat128NaN( a, a STATUS_VAR );
+ if (aSig0 | aSig1) {
+ return propagateFloat128NaN(a, a, status);
+ }
if ( ! aSign ) return a;
goto invalid;
}
if ( aSign ) {
if ( ( aExp | aSig0 | aSig1 ) == 0 ) return a;
invalid:
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
z.low = float128_default_nan_low;
z.high = float128_default_nan_high;
return z;
zSig1 |= ( ( rem1 | rem2 | rem3 ) != 0 );
}
shift128ExtraRightJamming( zSig0, zSig1, 0, 14, &zSig0, &zSig1, &zSig2 );
- return roundAndPackFloat128( 0, zExp, zSig0, zSig1, zSig2 STATUS_VAR );
+ return roundAndPackFloat128(0, zExp, zSig0, zSig1, zSig2, status);
}
|| ( ( extractFloat128Exp( b ) == 0x7FFF )
&& ( extractFloat128Frac0( b ) | extractFloat128Frac1( b ) ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
return
|| ( ( extractFloat128Exp( b ) == 0x7FFF )
&& ( extractFloat128Frac0( b ) | extractFloat128Frac1( b ) ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
aSign = extractFloat128Sign( a );
|| ( ( extractFloat128Exp( b ) == 0x7FFF )
&& ( extractFloat128Frac0( b ) | extractFloat128Frac1( b ) ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 0;
}
aSign = extractFloat128Sign( a );
|| ( ( extractFloat128Exp( b ) == 0x7FFF )
&& ( extractFloat128Frac0( b ) | extractFloat128Frac1( b ) ) )
) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return 1;
}
return 0;
) {
if ( float128_is_signaling_nan( a )
|| float128_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 0;
}
) {
if ( float128_is_signaling_nan( a )
|| float128_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 0;
}
) {
if ( float128_is_signaling_nan( a )
|| float128_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 0;
}
) {
if ( float128_is_signaling_nan( a )
|| float128_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return 1;
}
/* misc functions */
float32 uint32_to_float32(uint32_t a, float_status *status)
{
- return int64_to_float32(a STATUS_VAR);
+ return int64_to_float32(a, status);
}
float64 uint32_to_float64(uint32_t a, float_status *status)
{
- return int64_to_float64(a STATUS_VAR);
+ return int64_to_float64(a, status);
}
uint32 float32_to_uint32(float32 a, float_status *status)
uint32 res;
int old_exc_flags = get_float_exception_flags(status);
- v = float32_to_int64(a STATUS_VAR);
+ v = float32_to_int64(a, status);
if (v < 0) {
res = 0;
} else if (v > 0xffffffff) {
return v;
}
set_float_exception_flags(old_exc_flags, status);
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return res;
}
uint32 res;
int old_exc_flags = get_float_exception_flags(status);
- v = float32_to_int64_round_to_zero(a STATUS_VAR);
+ v = float32_to_int64_round_to_zero(a, status);
if (v < 0) {
res = 0;
} else if (v > 0xffffffff) {
return v;
}
set_float_exception_flags(old_exc_flags, status);
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return res;
}
int_fast16_t res;
int old_exc_flags = get_float_exception_flags(status);
- v = float32_to_int32(a STATUS_VAR);
+ v = float32_to_int32(a, status);
if (v < -0x8000) {
res = -0x8000;
} else if (v > 0x7fff) {
}
set_float_exception_flags(old_exc_flags, status);
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return res;
}
uint_fast16_t res;
int old_exc_flags = get_float_exception_flags(status);
- v = float32_to_int32(a STATUS_VAR);
+ v = float32_to_int32(a, status);
if (v < 0) {
res = 0;
} else if (v > 0xffff) {
}
set_float_exception_flags(old_exc_flags, status);
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return res;
}
uint_fast16_t res;
int old_exc_flags = get_float_exception_flags(status);
- v = float32_to_int64_round_to_zero(a STATUS_VAR);
+ v = float32_to_int64_round_to_zero(a, status);
if (v < 0) {
res = 0;
} else if (v > 0xffff) {
return v;
}
set_float_exception_flags(old_exc_flags, status);
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return res;
}
uint32 res;
int old_exc_flags = get_float_exception_flags(status);
- v = float64_to_uint64(a STATUS_VAR);
+ v = float64_to_uint64(a, status);
if (v > 0xffffffff) {
res = 0xffffffff;
} else {
return v;
}
set_float_exception_flags(old_exc_flags, status);
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return res;
}
uint32 res;
int old_exc_flags = get_float_exception_flags(status);
- v = float64_to_uint64_round_to_zero(a STATUS_VAR);
+ v = float64_to_uint64_round_to_zero(a, status);
if (v > 0xffffffff) {
res = 0xffffffff;
} else {
return v;
}
set_float_exception_flags(old_exc_flags, status);
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return res;
}
int_fast16_t res;
int old_exc_flags = get_float_exception_flags(status);
- v = float64_to_int32(a STATUS_VAR);
+ v = float64_to_int32(a, status);
if (v < -0x8000) {
res = -0x8000;
} else if (v > 0x7fff) {
}
set_float_exception_flags(old_exc_flags, status);
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return res;
}
uint_fast16_t res;
int old_exc_flags = get_float_exception_flags(status);
- v = float64_to_int32(a STATUS_VAR);
+ v = float64_to_int32(a, status);
if (v < 0) {
res = 0;
} else if (v > 0xffff) {
}
set_float_exception_flags(old_exc_flags, status);
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return res;
}
uint_fast16_t res;
int old_exc_flags = get_float_exception_flags(status);
- v = float64_to_int64_round_to_zero(a STATUS_VAR);
+ v = float64_to_int64_round_to_zero(a, status);
if (v < 0) {
res = 0;
} else if (v > 0xffff) {
return v;
}
set_float_exception_flags(old_exc_flags, status);
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return res;
}
flag aSign;
int_fast16_t aExp, shiftCount;
uint64_t aSig, aSigExtra;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aSig = extractFloat64Frac(a);
aExp = extractFloat64Exp(a);
aSign = extractFloat64Sign(a);
if (aSign && (aExp > 1022)) {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
if (float64_is_any_nan(a)) {
return LIT64(0xFFFFFFFFFFFFFFFF);
} else {
shiftCount = 0x433 - aExp;
if (shiftCount <= 0) {
if (0x43E < aExp) {
- float_raise(float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
return LIT64(0xFFFFFFFFFFFFFFFF);
}
aSigExtra = 0;
} else {
shift64ExtraRightJamming(aSig, 0, shiftCount, &aSig, &aSigExtra);
}
- return roundAndPackUint64(aSign, aSig, aSigExtra STATUS_VAR);
+ return roundAndPackUint64(aSign, aSig, aSigExtra, status);
}
uint64_t float64_to_uint64_round_to_zero(float64 a, float_status *status)
{
signed char current_rounding_mode = STATUS(float_rounding_mode);
- set_float_rounding_mode(float_round_to_zero STATUS_VAR);
- int64_t v = float64_to_uint64(a STATUS_VAR);
- set_float_rounding_mode(current_rounding_mode STATUS_VAR);
+ set_float_rounding_mode(float_round_to_zero, status);
+ int64_t v = float64_to_uint64(a, status);
+ set_float_rounding_mode(current_rounding_mode, status);
return v;
}
{ \
flag aSign, bSign; \
uint ## s ## _t av, bv; \
- a = float ## s ## _squash_input_denormal(a STATUS_VAR); \
- b = float ## s ## _squash_input_denormal(b STATUS_VAR); \
+ a = float ## s ## _squash_input_denormal(a, status); \
+ b = float ## s ## _squash_input_denormal(b, status); \
\
if (( ( extractFloat ## s ## Exp( a ) == nan_exp ) && \
extractFloat ## s ## Frac( a ) ) || \
if (!is_quiet || \
float ## s ## _is_signaling_nan( a ) || \
float ## s ## _is_signaling_nan( b ) ) { \
- float_raise( float_flag_invalid STATUS_VAR); \
+ float_raise(float_flag_invalid, status); \
} \
return float_relation_unordered; \
} \
\
int float ## s ## _compare(float ## s a, float ## s b, float_status *status) \
{ \
- return float ## s ## _compare_internal(a, b, 0 STATUS_VAR); \
+ return float ## s ## _compare_internal(a, b, 0, status); \
} \
\
int float ## s ## _compare_quiet(float ## s a, float ## s b, \
float_status *status) \
{ \
- return float ## s ## _compare_internal(a, b, 1 STATUS_VAR); \
+ return float ## s ## _compare_internal(a, b, 1, status); \
}
COMPARE(32, 0xff)
if (!is_quiet ||
floatx80_is_signaling_nan( a ) ||
floatx80_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return float_relation_unordered;
}
int floatx80_compare(floatx80 a, floatx80 b, float_status *status)
{
- return floatx80_compare_internal(a, b, 0 STATUS_VAR);
+ return floatx80_compare_internal(a, b, 0, status);
}
int floatx80_compare_quiet(floatx80 a, floatx80 b, float_status *status)
{
- return floatx80_compare_internal(a, b, 1 STATUS_VAR);
+ return floatx80_compare_internal(a, b, 1, status);
}
static inline int float128_compare_internal(float128 a, float128 b,
if (!is_quiet ||
float128_is_signaling_nan( a ) ||
float128_is_signaling_nan( b ) ) {
- float_raise( float_flag_invalid STATUS_VAR);
+ float_raise(float_flag_invalid, status);
}
return float_relation_unordered;
}
int float128_compare(float128 a, float128 b, float_status *status)
{
- return float128_compare_internal(a, b, 0 STATUS_VAR);
+ return float128_compare_internal(a, b, 0, status);
}
int float128_compare_quiet(float128 a, float128 b, float_status *status)
{
- return float128_compare_internal(a, b, 1 STATUS_VAR);
+ return float128_compare_internal(a, b, 1, status);
}
/* min() and max() functions. These can't be implemented as
{ \
flag aSign, bSign; \
uint ## s ## _t av, bv, aav, abv; \
- a = float ## s ## _squash_input_denormal(a STATUS_VAR); \
- b = float ## s ## _squash_input_denormal(b STATUS_VAR); \
+ a = float ## s ## _squash_input_denormal(a, status); \
+ b = float ## s ## _squash_input_denormal(b, status); \
if (float ## s ## _is_any_nan(a) || \
float ## s ## _is_any_nan(b)) { \
if (isieee) { \
return a; \
} \
} \
- return propagateFloat ## s ## NaN(a, b STATUS_VAR); \
+ return propagateFloat ## s ## NaN(a, b, status); \
} \
aSign = extractFloat ## s ## Sign(a); \
bSign = extractFloat ## s ## Sign(b); \
float ## s float ## s ## _min(float ## s a, float ## s b, \
float_status *status) \
{ \
- return float ## s ## _minmax(a, b, 1, 0, 0 STATUS_VAR); \
+ return float ## s ## _minmax(a, b, 1, 0, 0, status); \
} \
\
float ## s float ## s ## _max(float ## s a, float ## s b, \
float_status *status) \
{ \
- return float ## s ## _minmax(a, b, 0, 0, 0 STATUS_VAR); \
+ return float ## s ## _minmax(a, b, 0, 0, 0, status); \
} \
\
float ## s float ## s ## _minnum(float ## s a, float ## s b, \
float_status *status) \
{ \
- return float ## s ## _minmax(a, b, 1, 1, 0 STATUS_VAR); \
+ return float ## s ## _minmax(a, b, 1, 1, 0, status); \
} \
\
float ## s float ## s ## _maxnum(float ## s a, float ## s b, \
float_status *status) \
{ \
- return float ## s ## _minmax(a, b, 0, 1, 0 STATUS_VAR); \
+ return float ## s ## _minmax(a, b, 0, 1, 0, status); \
} \
\
float ## s float ## s ## _minnummag(float ## s a, float ## s b, \
float_status *status) \
{ \
- return float ## s ## _minmax(a, b, 1, 1, 1 STATUS_VAR); \
+ return float ## s ## _minmax(a, b, 1, 1, 1, status); \
} \
\
float ## s float ## s ## _maxnummag(float ## s a, float ## s b, \
float_status *status) \
{ \
- return float ## s ## _minmax(a, b, 0, 1, 1 STATUS_VAR); \
+ return float ## s ## _minmax(a, b, 0, 1, 1, status); \
}
MINMAX(32)
int16_t aExp;
uint32_t aSig;
- a = float32_squash_input_denormal(a STATUS_VAR);
+ a = float32_squash_input_denormal(a, status);
aSig = extractFloat32Frac( a );
aExp = extractFloat32Exp( a );
aSign = extractFloat32Sign( a );
if ( aExp == 0xFF ) {
if ( aSig ) {
- return propagateFloat32NaN( a, a STATUS_VAR );
+ return propagateFloat32NaN(a, a, status);
}
return a;
}
aExp += n - 1;
aSig <<= 7;
- return normalizeRoundAndPackFloat32( aSign, aExp, aSig STATUS_VAR );
+ return normalizeRoundAndPackFloat32(aSign, aExp, aSig, status);
}
float64 float64_scalbn(float64 a, int n, float_status *status)
int16_t aExp;
uint64_t aSig;
- a = float64_squash_input_denormal(a STATUS_VAR);
+ a = float64_squash_input_denormal(a, status);
aSig = extractFloat64Frac( a );
aExp = extractFloat64Exp( a );
aSign = extractFloat64Sign( a );
if ( aExp == 0x7FF ) {
if ( aSig ) {
- return propagateFloat64NaN( a, a STATUS_VAR );
+ return propagateFloat64NaN(a, a, status);
}
return a;
}
aExp += n - 1;
aSig <<= 10;
- return normalizeRoundAndPackFloat64( aSign, aExp, aSig STATUS_VAR );
+ return normalizeRoundAndPackFloat64(aSign, aExp, aSig, status);
}
floatx80 floatx80_scalbn(floatx80 a, int n, float_status *status)
if ( aExp == 0x7FFF ) {
if ( aSig<<1 ) {
- return propagateFloatx80NaN( a, a STATUS_VAR );
+ return propagateFloatx80NaN(a, a, status);
}
return a;
}
aExp += n;
return normalizeRoundAndPackFloatx80( STATUS(floatx80_rounding_precision),
- aSign, aExp, aSig, 0 STATUS_VAR );
+ aSign, aExp, aSig, 0, status);
}
float128 float128_scalbn(float128 a, int n, float_status *status)
aSign = extractFloat128Sign( a );
if ( aExp == 0x7FFF ) {
if ( aSig0 | aSig1 ) {
- return propagateFloat128NaN( a, a STATUS_VAR );
+ return propagateFloat128NaN(a, a, status);
}
return a;
}
aExp += n - 1;
return normalizeRoundAndPackFloat128( aSign, aExp, aSig0, aSig1
- STATUS_VAR );
+ , status);
}