\frac{x \cdot 2}{y \cdot z - t \cdot z}\begin{array}{l}
\mathbf{if}\;z \le -6.948586446622267662315685094511317832413 \cdot 10^{96}:\\
\;\;\;\;\frac{\frac{x \cdot 2}{z}}{y - t}\\
\mathbf{elif}\;z \le 138177673017854208086533707964897199915000:\\
\;\;\;\;\frac{x}{\frac{y - t}{\frac{2}{z}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{\sqrt[3]{z} \cdot \sqrt[3]{z}}}{\frac{y - t}{\frac{2}{\sqrt[3]{z}}}}\\
\end{array}double f(double x, double y, double z, double t) {
double r485613 = x;
double r485614 = 2.0;
double r485615 = r485613 * r485614;
double r485616 = y;
double r485617 = z;
double r485618 = r485616 * r485617;
double r485619 = t;
double r485620 = r485619 * r485617;
double r485621 = r485618 - r485620;
double r485622 = r485615 / r485621;
return r485622;
}
double f(double x, double y, double z, double t) {
double r485623 = z;
double r485624 = -6.948586446622268e+96;
bool r485625 = r485623 <= r485624;
double r485626 = x;
double r485627 = 2.0;
double r485628 = r485626 * r485627;
double r485629 = r485628 / r485623;
double r485630 = y;
double r485631 = t;
double r485632 = r485630 - r485631;
double r485633 = r485629 / r485632;
double r485634 = 1.381776730178542e+41;
bool r485635 = r485623 <= r485634;
double r485636 = r485627 / r485623;
double r485637 = r485632 / r485636;
double r485638 = r485626 / r485637;
double r485639 = cbrt(r485623);
double r485640 = r485639 * r485639;
double r485641 = r485626 / r485640;
double r485642 = r485627 / r485639;
double r485643 = r485632 / r485642;
double r485644 = r485641 / r485643;
double r485645 = r485635 ? r485638 : r485644;
double r485646 = r485625 ? r485633 : r485645;
return r485646;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.0 |
|---|---|
| Target | 2.2 |
| Herbie | 3.0 |
if z < -6.948586446622268e+96Initial program 14.7
Simplified11.1
rmApplied associate-/r*2.3
if -6.948586446622268e+96 < z < 1.381776730178542e+41Initial program 2.7
Simplified2.7
rmApplied associate-/r*8.5
rmApplied *-un-lft-identity8.5
Applied times-frac8.6
Applied associate-/l*2.8
if 1.381776730178542e+41 < z Initial program 11.7
Simplified9.4
rmApplied associate-/r*2.5
rmApplied add-cube-cbrt3.1
Applied times-frac3.1
Applied associate-/l*4.2
Final simplification3.0
herbie shell --seed 2019303
(FPCore (x y z t)
:name "Linear.Projection:infinitePerspective from linear-1.19.1.3, A"
:precision binary64
:herbie-target
(if (< (/ (* x 2) (- (* y z) (* t z))) -2.559141628295061e-13) (* (/ x (* (- y t) z)) 2) (if (< (/ (* x 2) (- (* y z) (* t z))) 1.045027827330126e-269) (/ (* (/ x z) 2) (- y t)) (* (/ x (* (- y t) z)) 2)))
(/ (* x 2) (- (* y z) (* t z))))