\frac{x}{x \cdot x + 1}\begin{array}{l}
\mathbf{if}\;x \le -7.4175050023885899 \cdot 10^{46} \lor \neg \left(x \le 458.065142995027088\right):\\
\;\;\;\;\mathsf{fma}\left(1, \frac{1}{{x}^{5}} - \frac{1}{{x}^{3}}, \frac{1}{x}\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{1}{x \cdot x + 1}\\
\end{array}double code(double x) {
return ((double) (x / ((double) (((double) (x * x)) + 1.0))));
}
double code(double x) {
double VAR;
if (((x <= -7.41750500238859e+46) || !(x <= 458.0651429950271))) {
VAR = ((double) fma(1.0, ((double) (((double) (1.0 / ((double) pow(x, 5.0)))) - ((double) (1.0 / ((double) pow(x, 3.0)))))), ((double) (1.0 / x))));
} else {
VAR = ((double) (x * ((double) (1.0 / ((double) (((double) (x * x)) + 1.0))))));
}
return VAR;
}




Bits error versus x
Results
| Original | 14.3 |
|---|---|
| Target | 0.1 |
| Herbie | 0.0 |
if x < -7.41750500238859e+46 or 458.0651429950271 < x Initial program 31.9
Taylor expanded around inf 0.0
Simplified0.0
if -7.41750500238859e+46 < x < 458.0651429950271Initial program 0.0
rmApplied div-inv0.0
Final simplification0.0
herbie shell --seed 2020113 +o rules:numerics
(FPCore (x)
:name "x / (x^2 + 1)"
:precision binary64
:herbie-target
(/ 1 (+ x (/ 1 x)))
(/ x (+ (* x x) 1)))