\frac{x}{x \cdot x + 1}\begin{array}{l}
\mathbf{if}\;x \le -4182345296.2656183 \lor \neg \left(x \le 50461.0592015108195\right):\\
\;\;\;\;1 \cdot \left(\frac{1}{{x}^{5}} - \frac{1}{{x}^{3}}\right) + \frac{1}{x}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{1}{x \cdot x + 1}\\
\end{array}double code(double x) {
return (x / ((x * x) + 1.0));
}
double code(double x) {
double temp;
if (((x <= -4182345296.2656183) || !(x <= 50461.05920151082))) {
temp = ((1.0 * ((1.0 / pow(x, 5.0)) - (1.0 / pow(x, 3.0)))) + (1.0 / x));
} else {
temp = (x * (1.0 / ((x * x) + 1.0)));
}
return temp;
}




Bits error versus x
Results
| Original | 15.1 |
|---|---|
| Target | 0.1 |
| Herbie | 0.0 |
if x < -4182345296.2656183 or 50461.05920151082 < x Initial program 30.8
Taylor expanded around inf 0.0
Simplified0.0
if -4182345296.2656183 < x < 50461.05920151082Initial program 0.0
rmApplied div-inv0.0
Final simplification0.0
herbie shell --seed 2020057
(FPCore (x)
:name "x / (x^2 + 1)"
:precision binary64
:herbie-target
(/ 1 (+ x (/ 1 x)))
(/ x (+ (* x x) 1)))