\frac{x}{x \cdot x + 1}\begin{array}{l}
\mathbf{if}\;x \leq -20496344444626.17:\\
\;\;\;\;\left(\frac{1}{x} + \frac{1}{{x}^{5}}\right) - {\left(\frac{1}{x}\right)}^{3}\\
\mathbf{elif}\;x \leq 456.806429230397:\\
\;\;\;\;\frac{x}{{x}^{4} + -1} \cdot \left(-1 + x \cdot x\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{1}{x} + \frac{1}{{x}^{5}}\right) - {\left(\frac{1}{x}\right)}^{3}\\
\end{array}(FPCore (x) :precision binary64 (/ x (+ (* x x) 1.0)))
(FPCore (x)
:precision binary64
(if (<= x -20496344444626.17)
(- (+ (/ 1.0 x) (/ 1.0 (pow x 5.0))) (pow (/ 1.0 x) 3.0))
(if (<= x 456.806429230397)
(* (/ x (+ (pow x 4.0) -1.0)) (+ -1.0 (* x x)))
(- (+ (/ 1.0 x) (/ 1.0 (pow x 5.0))) (pow (/ 1.0 x) 3.0)))))double code(double x) {
return x / ((x * x) + 1.0);
}
double code(double x) {
double tmp;
if (x <= -20496344444626.17) {
tmp = ((1.0 / x) + (1.0 / pow(x, 5.0))) - pow((1.0 / x), 3.0);
} else if (x <= 456.806429230397) {
tmp = (x / (pow(x, 4.0) + -1.0)) * (-1.0 + (x * x));
} else {
tmp = ((1.0 / x) + (1.0 / pow(x, 5.0))) - pow((1.0 / x), 3.0);
}
return tmp;
}




Bits error versus x
Results
| Original | 14.9 |
|---|---|
| Target | 0.1 |
| Herbie | 0.0 |
if x < -20496344444626.172 or 456.806429230396986 < x Initial program 30.4
Taylor expanded around inf 0.0
Simplified0.0
if -20496344444626.172 < x < 456.806429230396986Initial program 0.0
rmApplied flip-+_binary640.0
Applied associate-/r/_binary640.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2020273
(FPCore (x)
:name "x / (x^2 + 1)"
:precision binary64
:herbie-target
(/ 1.0 (+ x (/ 1.0 x)))
(/ x (+ (* x x) 1.0)))