\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 r68386 = x;
double r68387 = y;
double r68388 = r68386 - r68387;
double r68389 = r68386 + r68387;
double r68390 = r68388 * r68389;
double r68391 = r68386 * r68386;
double r68392 = r68387 * r68387;
double r68393 = r68391 + r68392;
double r68394 = r68390 / r68393;
return r68394;
}
double f(double x, double y) {
double r68395 = x;
double r68396 = y;
double r68397 = r68395 - r68396;
double r68398 = hypot(r68395, r68396);
double r68399 = r68397 / r68398;
double r68400 = r68395 + r68396;
double r68401 = r68398 / r68400;
double r68402 = r68399 / r68401;
return r68402;
}




Bits error versus x




Bits error versus y
Results
| Original | 19.9 |
|---|---|
| Target | 0.0 |
| Herbie | 0.0 |
Initial program 19.9
rmApplied add-sqr-sqrt19.9
Applied associate-/r*19.9
Simplified19.9
rmApplied hypot-def0.0
rmApplied associate-/r/0.0
Applied associate-/l*0.0
Final simplification0.0
herbie shell --seed 2020003 +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))))