\frac{\left(x \cdot 2\right) \cdot y}{x - y}\begin{array}{l}
\mathbf{if}\;y \le -1.93403941868859091 \cdot 10^{-36} \lor \neg \left(y \le 1.19064852231214522 \cdot 10^{-72}\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 r626673 = x;
double r626674 = 2.0;
double r626675 = r626673 * r626674;
double r626676 = y;
double r626677 = r626675 * r626676;
double r626678 = r626673 - r626676;
double r626679 = r626677 / r626678;
return r626679;
}
double f(double x, double y) {
double r626680 = y;
double r626681 = -1.934039418688591e-36;
bool r626682 = r626680 <= r626681;
double r626683 = 1.1906485223121452e-72;
bool r626684 = r626680 <= r626683;
double r626685 = !r626684;
bool r626686 = r626682 || r626685;
double r626687 = x;
double r626688 = 2.0;
double r626689 = r626687 * r626688;
double r626690 = r626687 - r626680;
double r626691 = r626690 / r626680;
double r626692 = r626689 / r626691;
double r626693 = r626687 / r626690;
double r626694 = r626680 * r626688;
double r626695 = r626693 * r626694;
double r626696 = r626686 ? r626692 : r626695;
return r626696;
}




Bits error versus x




Bits error versus y
Results
| Original | 15.3 |
|---|---|
| Target | 0.3 |
| Herbie | 0.3 |
if y < -1.934039418688591e-36 or 1.1906485223121452e-72 < y Initial program 14.4
rmApplied associate-/l*0.4
if -1.934039418688591e-36 < y < 1.1906485223121452e-72Initial program 16.6
rmApplied associate-/l*17.8
rmApplied div-inv17.9
Applied times-frac0.2
Simplified0.0
Final simplification0.3
herbie shell --seed 2020060
(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)))