riscv-pk/softfloat/extF80M_sqrt.c

61     uint64_t rem64;
76     rem64 = aSPtr->signif;
80         if ( rem64 & UINT64_C( 0x7FFFFFFFFFFFFFFF ) ) {
85         rem64 = UINT64_C( 0x8000000000000000 );
91     if ( ! (rem64 & UINT64_C( 0x8000000000000000 )) ) {
92         if ( ! rem64 ) {
96         expA += softfloat_normExtF80SigM( &rem64 );
103     softfloat_shortShiftLeft64To96M( rem64, 30 - expA, rem96 );
104     sig32A = rem64>>32;
108     rem64 =
111     rem96[indexWord( 3, 2 )] = rem64>>32;
112     rem96[indexWord( 3, 1 )] = rem64;
115     q = ((uint32_t) (rem64>>2) * (uint64_t) recipSqrt32)>>32;
132     rem64 = (uint64_t) rem32<<32 | rem[indexWord( 4, 2 )];
135     q = (((uint32_t) (rem64>>2) * (uint64_t) recipSqrt32)>>32) + 2;
136     if ( rem64>>34 ) q += recipSqrt32;
175     zSPtr->signif  = rem64;