1 - \frac{\left(1 - x\right) \cdot y}{y + 1}\begin{array}{l}
\mathbf{if}\;y \leq -13978.14375007929 \lor \neg \left(y \leq 13879.373838587384\right):\\
\;\;\;\;\left(\frac{x}{y \cdot y} + \left(x + \left(\frac{1}{y} + \frac{1}{{y}^{3}}\right)\right)\right) - \left(\frac{x}{{y}^{3}} + \left(\frac{x}{y} + \frac{1}{y \cdot y}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;1 - \frac{y}{y + 1} \cdot \left(1 - x\right)\\
\end{array}(FPCore (x y) :precision binary64 (- 1.0 (/ (* (- 1.0 x) y) (+ y 1.0))))
(FPCore (x y)
:precision binary64
(if (or (<= y -13978.14375007929) (not (<= y 13879.373838587384)))
(-
(+ (/ x (* y y)) (+ x (+ (/ 1.0 y) (/ 1.0 (pow y 3.0)))))
(+ (/ x (pow y 3.0)) (+ (/ x y) (/ 1.0 (* y y)))))
(- 1.0 (* (/ y (+ y 1.0)) (- 1.0 x)))))double code(double x, double y) {
return 1.0 - (((1.0 - x) * y) / (y + 1.0));
}
double code(double x, double y) {
double tmp;
if ((y <= -13978.14375007929) || !(y <= 13879.373838587384)) {
tmp = ((x / (y * y)) + (x + ((1.0 / y) + (1.0 / pow(y, 3.0))))) - ((x / pow(y, 3.0)) + ((x / y) + (1.0 / (y * y))));
} else {
tmp = 1.0 - ((y / (y + 1.0)) * (1.0 - x));
}
return tmp;
}




Bits error versus x




Bits error versus y
Results
| Original | 22.4 |
|---|---|
| Target | 0.3 |
| Herbie | 0.0 |
if y < -13978.1437500792908 or 13879.373838587384 < y Initial program 45.6
Taylor expanded around inf 0.0
Simplified0.0
if -13978.1437500792908 < y < 13879.373838587384Initial program 0.1
Taylor expanded around 0 0.1
Simplified0.1
Final simplification0.0
herbie shell --seed 2021023
(FPCore (x y)
:name "Diagrams.Trail:splitAtParam from diagrams-lib-1.3.0.3, D"
:precision binary64
:herbie-target
(if (< y -3693.8482788297247) (- (/ 1.0 y) (- (/ x y) x)) (if (< y 6799310503.41891) (- 1.0 (/ (* (- 1.0 x) y) (+ y 1.0))) (- (/ 1.0 y) (- (/ x y) x))))
(- 1.0 (/ (* (- 1.0 x) y) (+ y 1.0))))