--- /dev/null
+\r
+#include <stdbool.h>\r
+#include <stdint.h>\r
+#include "platform.h"\r
+#include "primitives.h"\r
+#include "internals.h"\r
+#include "specialize.h"\r
+#include "softfloat.h"\r
+\r
+float32_t\r
+ softfloat_mulAddF32(\r
+ int op, uint_fast32_t uiA, uint_fast32_t uiB, uint_fast32_t uiC )\r
+{\r
+ bool signA;\r
+ int_fast16_t expA;\r
+ uint_fast32_t sigA;\r
+ bool signB;\r
+ int_fast16_t expB;\r
+ uint_fast32_t sigB;\r
+ bool signC;\r
+ int_fast16_t expC;\r
+ uint_fast32_t sigC;\r
+ bool signProd;\r
+ uint_fast32_t magBits, uiZ;\r
+ struct exp16_sig32 normExpSig;\r
+ int_fast16_t expProd;\r
+ uint_fast64_t sigProd;\r
+ bool signZ;\r
+ int_fast16_t expZ;\r
+ uint_fast32_t sigZ;\r
+ int_fast16_t expDiff;\r
+ uint_fast64_t sigZ64, sigC64;\r
+ int shiftCount;\r
+ union ui32_f32 uZ;\r
+\r
+ signA = signF32UI( uiA );\r
+ expA = expF32UI( uiA );\r
+ sigA = fracF32UI( uiA );\r
+ signB = signF32UI( uiB );\r
+ expB = expF32UI( uiB );\r
+ sigB = fracF32UI( uiB );\r
+ signC = signF32UI( uiC ) ^ ( op == softfloat_mulAdd_subC );\r
+ expC = expF32UI( uiC );\r
+ sigC = fracF32UI( uiC );\r
+ signProd = signA ^ signB ^ ( op == softfloat_mulAdd_subProd );\r
+ if ( expA == 0xFF ) {\r
+ if ( sigA || ( ( expB == 0xFF ) && sigB ) ) goto propagateNaN_ABC;\r
+ magBits = expB | sigB;\r
+ goto infProdArg;\r
+ }\r
+ if ( expB == 0xFF ) {\r
+ if ( sigB ) goto propagateNaN_ABC;\r
+ magBits = expA | sigA;\r
+ goto infProdArg;\r
+ }\r
+ if ( expC == 0xFF ) {\r
+ if ( sigC ) {\r
+ uiZ = 0;\r
+ goto propagateNaN_ZC;\r
+ }\r
+ uiZ = uiC;\r
+ goto uiZ;\r
+ }\r
+ if ( ! expA ) {\r
+ if ( ! sigA ) goto zeroProd;\r
+ normExpSig = softfloat_normSubnormalF32Sig( sigA );\r
+ expA = normExpSig.exp;\r
+ sigA = normExpSig.sig;\r
+ }\r
+ if ( ! expB ) {\r
+ if ( ! sigB ) goto zeroProd;\r
+ normExpSig = softfloat_normSubnormalF32Sig( sigB );\r
+ expB = normExpSig.exp;\r
+ sigB = normExpSig.sig;\r
+ }\r
+ expProd = expA + expB - 0x7E;\r
+ sigA = ( sigA | 0x00800000 )<<7;\r
+ sigB = ( sigB | 0x00800000 )<<7;\r
+ sigProd = (uint_fast64_t) sigA * sigB;\r
+ if ( sigProd < UINT64_C( 0x2000000000000000 ) ) {\r
+ --expProd;\r
+ sigProd <<= 1;\r
+ }\r
+ signZ = signProd;\r
+ if ( ! expC ) {\r
+ if ( ! sigC ) {\r
+ expZ = expProd - 1;\r
+ sigZ = softfloat_shortShift64RightJam( sigProd, 31 );\r
+ goto roundPack;\r
+ }\r
+ normExpSig = softfloat_normSubnormalF32Sig( sigC );\r
+ expC = normExpSig.exp;\r
+ sigC = normExpSig.sig;\r
+ }\r
+ sigC = ( sigC | 0x00800000 )<<6;\r
+ expDiff = expProd - expC;\r
+ if ( signProd == signC ) {\r
+ if ( expDiff <= 0 ) {\r
+ expZ = expC;\r
+ sigZ = sigC + softfloat_shift64RightJam( sigProd, 32 - expDiff );\r
+ } else {\r
+ expZ = expProd;\r
+ sigZ64 =\r
+ sigProd\r
+ + softfloat_shift64RightJam(\r
+ (uint_fast64_t) sigC<<32, expDiff );\r
+ sigZ = softfloat_shortShift64RightJam( sigZ64, 32 );\r
+ }\r
+ if ( sigZ < 0x40000000 ) {\r
+ --expZ;\r
+ sigZ <<= 1;\r
+ }\r
+ } else {\r
+/*** OPTIMIZE BETTER? ***/\r
+ sigC64 = (uint_fast64_t) sigC<<32;\r
+ if ( expDiff < 0 ) {\r
+ signZ = signC;\r
+ expZ = expC;\r
+ sigZ64 = sigC64 - softfloat_shift64RightJam( sigProd, - expDiff );\r
+ } else if ( ! expDiff ) {\r
+ expZ = expProd;\r
+ sigZ64 = sigProd - sigC64;\r
+ if ( ! sigZ64 ) goto completeCancellation;\r
+ if ( sigZ64 & UINT64_C( 0x8000000000000000 ) ) {\r
+ signZ ^= 1;\r
+ sigZ64 = - sigZ64;\r
+ }\r
+ } else {\r
+ expZ = expProd;\r
+ sigZ64 = sigProd - softfloat_shift64RightJam( sigC64, expDiff );\r
+ }\r
+ shiftCount = softfloat_countLeadingZeros64( sigZ64 ) - 1;\r
+ expZ -= shiftCount;\r
+ shiftCount -= 32;\r
+ if ( shiftCount < 0 ) {\r
+ sigZ = softfloat_shortShift64RightJam( sigZ64, - shiftCount );\r
+ } else {\r
+ sigZ = (uint_fast32_t) sigZ64<<shiftCount;\r
+ }\r
+ }\r
+ roundPack:\r
+ return softfloat_roundPackToF32( signZ, expZ, sigZ );\r
+ propagateNaN_ABC:\r
+ uiZ = softfloat_propagateNaNF32UI( uiA, uiB );\r
+ goto propagateNaN_ZC;\r
+ infProdArg:\r
+ if ( magBits ) {\r
+ uiZ = packToF32UI( signProd, 0xFF, 0 );\r
+ if ( expC != 0xFF ) goto uiZ;\r
+ if ( sigC ) goto propagateNaN_ZC;\r
+ if ( signProd == signC ) goto uiZ;\r
+ }\r
+ invalid:\r
+ softfloat_raiseFlags( softfloat_flag_invalid );\r
+ uiZ = defaultNaNF32UI;\r
+ propagateNaN_ZC:\r
+ uiZ = softfloat_propagateNaNF32UI( uiZ, uiC );\r
+ goto uiZ;\r
+ zeroProd:\r
+ uiZ = uiC;\r
+ if ( ! ( expC | sigC ) && ( signProd != signC ) ) {\r
+ completeCancellation:\r
+ uiZ =\r
+ packToF32UI( softfloat_roundingMode == softfloat_round_min, 0, 0 );\r
+ }\r
+ uiZ:\r
+ uZ.ui = uiZ;\r
+ return uZ.f;\r
+\r
+}\r
+\r