Searched refs:reciprocal (Results 1 - 3 of 3) sorted by relevance

/freebsd-10.0-release/contrib/compiler-rt/lib/
H A Ddivsf3.c79 // [1, 2.0) and get a Q32 approximate reciprocal using a small minimax
80 // polynomial approximation: reciprocal = 3/4 + 1/sqrt(2) - b/2. This
83 uint32_t reciprocal = UINT32_C(0x7504f333) - q31b; local
85 // Now refine the reciprocal estimate using a Newton-Raphson iteration:
93 correction = -((uint64_t)reciprocal * q31b >> 32);
94 reciprocal = (uint64_t)reciprocal * correction >> 31;
95 correction = -((uint64_t)reciprocal * q31b >> 32);
96 reciprocal = (uint64_t)reciprocal * correctio
[all...]
H A Ddivdf3.c79 // [1, 2.0) and get a Q32 approximate reciprocal using a small minimax
80 // polynomial approximation: reciprocal = 3/4 + 1/sqrt(2) - b/2. This
85 // Now refine the reciprocal estimate using a Newton-Raphson iteration:
109 uint64_t correction, reciprocal; local
113 reciprocal = (uint64_t)recip32*cHi + ((uint64_t)recip32*cLo >> 32);
116 // 64-bit reciprocal estimate downward to ensure that it is strictly smaller
117 // than the infinitely precise exact reciprocal. Because the computation
120 reciprocal -= 2;
122 // The numerical reciprocal is accurate to within 2^-56, lies in the
123 // interval [0.5, 1.0), and is strictly smaller than the true reciprocal
[all...]
/freebsd-10.0-release/contrib/llvm/lib/Support/
H A DAPFloat.cpp2346 /* Extra half-ulp lost in reciprocal of exponent. */
3576 APFloat reciprocal(*semantics, 1ULL);
3577 if (reciprocal.divide(*this, rmNearestTiesToEven) != opOK)
3582 if (reciprocal.significandMSB() + 1 < reciprocal.semantics->precision)
3585 assert(reciprocal.category == fcNormal &&
3586 reciprocal.significandLSB() == reciprocal.semantics->precision - 1);
3589 *inv = reciprocal;

Completed in 66 milliseconds