\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(2 \cdot y\right)\\
\end{array}double f(double x, double y) {
double r582799 = x;
double r582800 = 2.0;
double r582801 = r582799 * r582800;
double r582802 = y;
double r582803 = r582801 * r582802;
double r582804 = r582799 - r582802;
double r582805 = r582803 / r582804;
return r582805;
}
double f(double x, double y) {
double r582806 = y;
double r582807 = -2.985649314254737e+16;
bool r582808 = r582806 <= r582807;
double r582809 = 2.1640309147455017e-67;
bool r582810 = r582806 <= r582809;
double r582811 = !r582810;
bool r582812 = r582808 || r582811;
double r582813 = x;
double r582814 = 2.0;
double r582815 = r582813 * r582814;
double r582816 = r582813 - r582806;
double r582817 = r582816 / r582806;
double r582818 = r582815 / r582817;
double r582819 = r582813 / r582816;
double r582820 = r582814 * r582806;
double r582821 = r582819 * r582820;
double r582822 = r582812 ? r582818 : r582821;
return r582822;
}




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)))