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

/freebsd-13-stable/contrib/llvm-project/compiler-rt/lib/builtins/
H A Ddivsf3.c88 // [1, 2.0) and get a Q32 approximate reciprocal using a small minimax
89 // polynomial approximation: reciprocal = 3/4 + 1/sqrt(2) - b/2. This
92 uint32_t reciprocal = UINT32_C(0x7504f333) - q31b; local
94 // Now refine the reciprocal estimate using a Newton-Raphson iteration:
101 correction = -((uint64_t)reciprocal * q31b >> 32);
102 reciprocal = (uint64_t)reciprocal * correction >> 31;
103 correction = -((uint64_t)reciprocal * q31b >> 32);
104 reciprocal = (uint64_t)reciprocal * correctio
[all...]
H A Ddivdf3.c87 // [1, 2.0) and get a Q32 approximate reciprocal using a small minimax
88 // polynomial approximation: reciprocal = 3/4 + 1/sqrt(2) - b/2. This
94 // Now refine the reciprocal estimate using a Newton-Raphson iteration:
108 // The reciprocal may have overflowed to zero if the upper half of b is
110 // computation that follows, so we adjust the reciprocal down by one bit.
116 uint64_t correction, reciprocal; local
120 reciprocal = (uint64_t)recip32 * cHi + ((uint64_t)recip32 * cLo >> 32);
122 // Adjust the final 64-bit reciprocal estimate downward to ensure that it is
123 // strictly smaller than the infinitely precise exact reciprocal. Because
126 reciprocal
[all...]
H A Ddivtf3.c88 // [1, 2.0) and get a Q64 approximate reciprocal using a small minimax
89 // polynomial approximation: reciprocal = 3/4 + 1/sqrt(2) - b/2. This
95 // Now refine the reciprocal estimate using a Newton-Raphson iteration:
113 // The reciprocal may have overflowed to zero if the upper half of b is
115 // computation that follows, so we adjust the reciprocal down by one bit.
121 rep_t correction, reciprocal; local
137 reciprocal = r64cH + (r64cL >> 64);
139 // Adjust the final 128-bit reciprocal estimate downward to ensure that it
140 // is strictly smaller than the infinitely precise exact reciprocal. Because
143 reciprocal
[all...]
/freebsd-13-stable/contrib/llvm-project/llvm/lib/Support/
H A DAPFloat.cpp2663 /* Extra half-ulp lost in reciprocal of exponent. */
4002 IEEEFloat reciprocal(*semantics, 1ULL);
4003 if (reciprocal.divide(*this, rmNearestTiesToEven) != opOK)
4008 if (reciprocal.isDenormal())
4011 assert(reciprocal.isFiniteNonZero() &&
4012 reciprocal.significandLSB() == reciprocal.semantics->precision - 1);
4015 *inv = APFloat(reciprocal, *semantics);

Completed in 113 milliseconds