\frac{\left(x \cdot 2\right) \cdot y}{x - y}
{\left(\frac{-0.5}{x} + \frac{0.5}{y}\right)}^{-1}
(FPCore (x y) :precision binary64 (/ (* (* x 2.0) y) (- x y)))
(FPCore (x y) :precision binary64 (pow (+ (/ -0.5 x) (/ 0.5 y)) -1.0))
double code(double x, double y) {
return ((x * 2.0) * y) / (x - y);
}
double code(double x, double y) {
return pow(((-0.5 / x) + (0.5 / y)), -1.0);
}




Bits error versus x




Bits error versus y
Results
| Original | 15.7 |
|---|---|
| Target | 0.3 |
| Herbie | 0.2 |
Initial program 15.7
Simplified7.4
Applied egg-rr7.6
Taylor expanded in x around 0 0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2022130
(FPCore (x y)
:name "Linear.Projection:perspective from linear-1.19.1.3, B"
:precision binary64
:herbie-target
(if (< x -1.7210442634149447e+81) (* (/ (* 2.0 x) (- x y)) y) (if (< x 83645045635564430.0) (/ (* x 2.0) (/ (- x y) y)) (* (/ (* 2.0 x) (- x y)) y)))
(/ (* (* x 2.0) y) (- x y)))