\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 r377998 = x;
double r377999 = 2.0;
double r378000 = r377998 * r377999;
double r378001 = y;
double r378002 = z;
double r378003 = r378001 * r378002;
double r378004 = t;
double r378005 = r378004 * r378002;
double r378006 = r378003 - r378005;
double r378007 = r378000 / r378006;
return r378007;
}
double f(double x, double y, double z, double t) {
double r378008 = z;
double r378009 = -6.948586446622268e+96;
bool r378010 = r378008 <= r378009;
double r378011 = x;
double r378012 = 2.0;
double r378013 = r378011 * r378012;
double r378014 = r378013 / r378008;
double r378015 = y;
double r378016 = t;
double r378017 = r378015 - r378016;
double r378018 = r378014 / r378017;
double r378019 = 1.381776730178542e+41;
bool r378020 = r378008 <= r378019;
double r378021 = r378012 / r378008;
double r378022 = r378017 / r378021;
double r378023 = r378011 / r378022;
double r378024 = cbrt(r378008);
double r378025 = r378024 * r378024;
double r378026 = r378011 / r378025;
double r378027 = r378012 / r378024;
double r378028 = r378017 / r378027;
double r378029 = r378026 / r378028;
double r378030 = r378020 ? r378023 : r378029;
double r378031 = r378010 ? r378018 : r378030;
return r378031;
}




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