\frac{x}{x \cdot x + 1}\begin{array}{l}
\mathbf{if}\;x \leq -5.176580413214303 \cdot 10^{+23} \lor \neg \left(x \leq 56166425.491830505\right):\\
\;\;\;\;\frac{1}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{1 + x \cdot x}\\
\end{array}(FPCore (x) :precision binary64 (/ x (+ (* x x) 1.0)))
(FPCore (x) :precision binary64 (if (or (<= x -5.176580413214303e+23) (not (<= x 56166425.491830505))) (/ 1.0 x) (/ x (+ 1.0 (* x x)))))
double code(double x) {
return x / ((x * x) + 1.0);
}
double code(double x) {
double tmp;
if ((x <= -5.176580413214303e+23) || !(x <= 56166425.491830505)) {
tmp = 1.0 / x;
} else {
tmp = x / (1.0 + (x * x));
}
return tmp;
}







Bits error versus x
Results
| Original | 14.9 |
|---|---|
| Target | 0.1 |
| Herbie | 0.0 |
| Alternative 1 | |
|---|---|
| Error | 0.1 |
| Cost | 448 |
| Alternative 2 | |
|---|---|
| Error | 0.6 |
| Cost | 520 |
| Alternative 3 | |
|---|---|
| Error | 31.1 |
| Cost | 64 |


if x < -5.1765804132143032e23 or 56166425.4918305054 < x Initial program 31.4
Taylor expanded around inf 0.0
Simplified0.0
if -5.1765804132143032e23 < x < 56166425.4918305054Initial program 0.0
rmApplied *-un-lft-identity_binary64_14420.0
Applied *-un-lft-identity_binary64_14420.0
Applied times-frac_binary64_14480.0
Simplified0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2021065
(FPCore (x)
:name "x / (x^2 + 1)"
:precision binary64
:herbie-target
(/ 1.0 (+ x (/ 1.0 x)))
(/ x (+ (* x x) 1.0)))