float64 float64_mul( float64 a, float64 b ) 
float64 float64_mul( struct roundingData *roundData, float64 a, float64 b ) 
float64 float64_mul( float64 a, float64 b ){    flag aSign, bSign, zSign;    int16 aExp, bExp, zExp;    bits64 aSig, bSig, zSig0, zSig1;    aSig = extractFloat64Frac( a );    aExp = extractFloat64Exp( a );    aSign = extractFloat64Sign( a );    bSig = extractFloat64Frac( b );    bExp = extractFloat64Exp( b );    bSign = extractFloat64Sign( b );    zSign = aSign ^ bSign;    if ( aExp == 0x7FF ) {        if ( aSig || ( ( bExp == 0x7FF ) && bSig ) ) {            return propagateFloat64NaN( a, b );        }        if ( ( bExp | bSig ) == 0 ) {            float_raise( float_flag_invalid );            return float64_default_nan;        }        return packFloat64( zSign, 0x7FF, 0 );    }    if ( bExp == 0x7FF ) {        if ( bSig ) return propagateFloat64NaN( a, b );        if ( ( aExp | aSig ) == 0 ) {            float_raise( float_flag_invalid );            return float64_default_nan;        }        return packFloat64( zSign, 0x7FF, 0 );    }    if ( aExp == 0 ) {        if ( aSig == 0 ) return packFloat64( zSign, 0, 0 );        normalizeFloat64Subnormal( aSig, &aExp, &aSig );    }    if ( bExp == 0 ) {        if ( bSig == 0 ) return packFloat64( zSign, 0, 0 );        normalizeFloat64Subnormal( bSig, &bExp, &bSig );    }    zExp = aExp + bExp - 0x3FF;    aSig = ( aSig | LIT64( 0x0010000000000000 ) )<<10;    bSig = ( bSig | LIT64( 0x0010000000000000 ) )<<11;    mul64To128( aSig, bSig, &zSig0, &zSig1 );    zSig0 |= ( zSig1 != 0 );    if ( 0 <= (sbits64) ( zSig0<<1 ) ) {        zSig0 <<= 1;        --zExp;    }    return roundAndPackFloat64( zSign, zExp, zSig0 );}
10457	/local/tlutelli/issta_data/temp/c/2005_temp/2005/10457/f148af2593ef76ac705d1cc6abe48f455c9912cc/softfloat.c/clean/arch/arm/nwfpe/softfloat.c