\frac{\left(x - y\right) \cdot \left(x + y\right)}{x \cdot x + y \cdot y}\frac{\frac{x - y}{\mathsf{hypot}\left(x, y\right)}}{\frac{\mathsf{hypot}\left(x, y\right)}{x + y}}double f(double x, double y) {
double r114636 = x;
double r114637 = y;
double r114638 = r114636 - r114637;
double r114639 = r114636 + r114637;
double r114640 = r114638 * r114639;
double r114641 = r114636 * r114636;
double r114642 = r114637 * r114637;
double r114643 = r114641 + r114642;
double r114644 = r114640 / r114643;
return r114644;
}
double f(double x, double y) {
double r114645 = x;
double r114646 = y;
double r114647 = r114645 - r114646;
double r114648 = hypot(r114645, r114646);
double r114649 = r114647 / r114648;
double r114650 = r114645 + r114646;
double r114651 = r114648 / r114650;
double r114652 = r114649 / r114651;
return r114652;
}




Bits error versus x




Bits error versus y
Results
| Original | 20.2 |
|---|---|
| Target | 0.1 |
| Herbie | 0.0 |
Initial program 20.2
rmApplied add-sqr-sqrt20.2
Applied associate-/r*20.3
Simplified20.3
rmApplied associate-/r/20.3
Applied associate-/l*20.3
Simplified0.0
Final simplification0.0
herbie shell --seed 2020062 +o rules:numerics
(FPCore (x y)
:name "Kahan p9 Example"
:precision binary64
:pre (and (< 0.0 x 1) (< y 1))
:herbie-target
(if (< 0.5 (fabs (/ x y)) 2) (/ (* (- x y) (+ x y)) (+ (* x x) (* y y))) (- 1 (/ 2 (+ 1 (* (/ x y) (/ x y))))))
(/ (* (- x y) (+ x y)) (+ (* x x) (* y y))))