1 - \frac{\left(1 - x\right) \cdot y}{y + 1}\begin{array}{l}
\mathbf{if}\;y \leq -417680595896.5533:\\
\;\;\;\;x + \frac{1}{y}\\
\mathbf{elif}\;y \leq 89866444.08949669:\\
\;\;\;\;1 - \left(y \cdot y + \left(1 - y\right)\right) \cdot \frac{y \cdot \left(1 - x\right)}{1 + {y}^{3}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{1 - x}{y}\\
\end{array}(FPCore (x y) :precision binary64 (- 1.0 (/ (* (- 1.0 x) y) (+ y 1.0))))
(FPCore (x y)
:precision binary64
(if (<= y -417680595896.5533)
(+ x (/ 1.0 y))
(if (<= y 89866444.08949669)
(- 1.0 (* (+ (* y y) (- 1.0 y)) (/ (* y (- 1.0 x)) (+ 1.0 (pow y 3.0)))))
(+ x (/ (- 1.0 x) y)))))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 <= -417680595896.5533) {
tmp = x + (1.0 / y);
} else if (y <= 89866444.08949669) {
tmp = 1.0 - (((y * y) + (1.0 - y)) * ((y * (1.0 - x)) / (1.0 + pow(y, 3.0))));
} else {
tmp = x + ((1.0 - x) / y);
}
return tmp;
}












Bits error versus x












Bits error versus y
Results
| Original | 22.0 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
| Alternative 1 | |
|---|---|
| Error | 0.2 |
| Cost | 8008 |
| Alternative 2 | |
|---|---|
| Error | 0.5 |
| Cost | 1346 |
| Alternative 3 | |
|---|---|
| Error | 0.8 |
| Cost | 904 |
| Alternative 4 | |
|---|---|
| Error | 1.0 |
| Cost | 776 |
| Alternative 5 | |
|---|---|
| Error | 9.1 |
| Cost | 776 |
| Alternative 6 | |
|---|---|
| Error | 9.3 |
| Cost | 648 |
| Alternative 7 | |
|---|---|
| Error | 16.8 |
| Cost | 706 |
| Alternative 8 | |
|---|---|
| Error | 38.9 |
| Cost | 64 |

if y < -417680595896.553284Initial program 46.1
rmApplied flip-+_binary64_1846650.1
Applied associate-/r/_binary64_1843850.2
Simplified50.2
Taylor expanded around inf 0.0
Simplified0.0
Taylor expanded around 0 0.1
Simplified0.1
if -417680595896.553284 < y < 89866444.0894966871Initial program 0.2
rmApplied flip3-+_binary64_184950.3
Applied associate-/r/_binary64_184380.3
Simplified0.3
Simplified0.3
if 89866444.0894966871 < y Initial program 45.9
rmApplied flip-+_binary64_1846650.6
Applied associate-/r/_binary64_1843850.6
Simplified50.6
Taylor expanded around inf 0.1
Simplified0.1
Simplified0.1
Final simplification0.2
herbie shell --seed 2021044
(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))))