\frac{\left(x \cdot 2\right) \cdot y}{x - y}\begin{array}{l}
\mathbf{if}\;y \le -29856493142547368 \lor \neg \left(y \le 2.1640309147455017 \cdot 10^{-67}\right):\\
\;\;\;\;\frac{x \cdot 2}{\frac{x - y}{y}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x - y} \cdot \left(y \cdot 2\right)\\
\end{array}double f(double x, double y) {
double r2048961 = x;
double r2048962 = 2.0;
double r2048963 = r2048961 * r2048962;
double r2048964 = y;
double r2048965 = r2048963 * r2048964;
double r2048966 = r2048961 - r2048964;
double r2048967 = r2048965 / r2048966;
return r2048967;
}
double f(double x, double y) {
double r2048968 = y;
double r2048969 = -2.985649314254737e+16;
bool r2048970 = r2048968 <= r2048969;
double r2048971 = 2.1640309147455017e-67;
bool r2048972 = r2048968 <= r2048971;
double r2048973 = !r2048972;
bool r2048974 = r2048970 || r2048973;
double r2048975 = x;
double r2048976 = 2.0;
double r2048977 = r2048975 * r2048976;
double r2048978 = r2048975 - r2048968;
double r2048979 = r2048978 / r2048968;
double r2048980 = r2048977 / r2048979;
double r2048981 = r2048975 / r2048978;
double r2048982 = r2048968 * r2048976;
double r2048983 = r2048981 * r2048982;
double r2048984 = r2048974 ? r2048980 : r2048983;
return r2048984;
}




Bits error versus x




Bits error versus y
Results
| Original | 15.3 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
if y < -2.985649314254737e+16 or 2.1640309147455017e-67 < y Initial program 15.3
rmApplied associate-/l*0.4
if -2.985649314254737e+16 < y < 2.1640309147455017e-67Initial program 15.3
rmApplied associate-/l*16.1
rmApplied div-inv16.2
Applied times-frac0.2
Simplified0.0
Final simplification0.2
herbie shell --seed 2020047 +o rules:numerics
(FPCore (x y)
:name "Linear.Projection:perspective from linear-1.19.1.3, B"
:precision binary64
:herbie-target
(if (< x -1.7210442634149447e+81) (* (/ (* 2 x) (- x y)) y) (if (< x 83645045635564432) (/ (* x 2) (/ (- x y) y)) (* (/ (* 2 x) (- x y)) y)))
(/ (* (* x 2) y) (- x y)))