3073     bits64 aSig0, aSig1, bSig;  in floatx80_rem()  local
3077     aSig0 = extractFloatx80Frac( a );  in floatx80_rem()
3084         if (    (bits64) ( aSig0<<1 )  in floatx80_rem()
3106         if ( (bits64) ( aSig0<<1 ) == 0 ) return a;  in floatx80_rem()
3107         normalizeFloatx80Subnormal( aSig0, &aExp, &aSig0 );  in floatx80_rem()
3115         shift128Right( aSig0, 0, 1, &aSig0, &aSig1 );  in floatx80_rem()
3118     q = ( bSig <= aSig0 );  in floatx80_rem()
3119     if ( q ) aSig0 -= bSig;  in floatx80_rem()
3122         q = estimateDiv128To64( aSig0, aSig1, bSig );  in floatx80_rem()
3125         sub128( aSig0, aSig1, term0, term1, &aSig0, &aSig1 );  in floatx80_rem()
3126         shortShift128Left( aSig0, aSig1, 62, &aSig0, &aSig1 );  in floatx80_rem()
3131         q = estimateDiv128To64( aSig0, aSig1, bSig );  in floatx80_rem()
3135         sub128( aSig0, aSig1, term0, term1, &aSig0, &aSig1 );  in floatx80_rem()
3137         while ( le128( term0, term1, aSig0, aSig1 ) ) {  in floatx80_rem()
3139             sub128( aSig0, aSig1, term0, term1, &aSig0, &aSig1 );  in floatx80_rem()
3146     sub128( term0, term1, aSig0, aSig1, &alternateASig0, &alternateASig1 );  in floatx80_rem()
3147     if (    lt128( alternateASig0, alternateASig1, aSig0, aSig1 )  in floatx80_rem()
3148          || (    eq128( alternateASig0, alternateASig1, aSig0, aSig1 )  in floatx80_rem()
3151         aSig0 = alternateASig0;  in floatx80_rem()
3158             roundData, zSign, bExp + expDiff, aSig0, aSig1 );  in floatx80_rem()
3173     bits64 aSig0, aSig1, zSig0, zSig1;  in floatx80_sqrt()  local
3178     aSig0 = extractFloatx80Frac( a );  in floatx80_sqrt()
3182         if ( (bits64) ( aSig0<<1 ) ) return propagateFloatx80NaN( a, a );  in floatx80_sqrt()
3187         if ( ( aExp | aSig0 ) == 0 ) return a;  in floatx80_sqrt()
3196         if ( aSig0 == 0 ) return packFloatx80( 0, 0, 0 );  in floatx80_sqrt()
3197         normalizeFloatx80Subnormal( aSig0, &aExp, &aSig0 );  in floatx80_sqrt()
3200     zSig0 = estimateSqrt32( aExp, aSig0>>32 );  in floatx80_sqrt()
3203     shift128Right( aSig0, 0, ( aExp & 1 ) + 2, &aSig0, &aSig1 );  in floatx80_sqrt()
3204     zSig0 = estimateDiv128To64( aSig0, aSig1, zSig0 ) + zSig0 + 4;  in floatx80_sqrt()
3206     shortShift128Left( aSig0, aSig1, 2, &aSig0, &aSig1 );  in floatx80_sqrt()
3208     sub128( aSig0, aSig1, term0, term1, &rem0, &rem1 );  in floatx80_sqrt()