\frac{x \cdot 2}{y \cdot z - t \cdot z}\begin{array}{l}
\mathbf{if}\;y \le -2.558225734480629833036083443548553390336 \cdot 10^{208}:\\
\;\;\;\;\left(x \cdot 2\right) \cdot \frac{\frac{1}{z}}{y - t}\\
\mathbf{elif}\;y \le 6.171057382867539840170668000327560127887 \cdot 10^{-299}:\\
\;\;\;\;\frac{\frac{x \cdot 2}{z}}{y - t}\\
\mathbf{elif}\;y \le 1.951130067970849161540972989726809538744 \cdot 10^{111}:\\
\;\;\;\;\frac{x \cdot 2}{y \cdot z + z \cdot \left(-t\right)}\\
\mathbf{elif}\;y \le 1.538759424307254056311250070293585704887 \cdot 10^{171}:\\
\;\;\;\;\frac{\frac{x \cdot 2}{z}}{y - t}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{2 \cdot x}{y - t}}{z}\\
\end{array}double f(double x, double y, double z, double t) {
double r362913 = x;
double r362914 = 2.0;
double r362915 = r362913 * r362914;
double r362916 = y;
double r362917 = z;
double r362918 = r362916 * r362917;
double r362919 = t;
double r362920 = r362919 * r362917;
double r362921 = r362918 - r362920;
double r362922 = r362915 / r362921;
return r362922;
}
double f(double x, double y, double z, double t) {
double r362923 = y;
double r362924 = -2.55822573448063e+208;
bool r362925 = r362923 <= r362924;
double r362926 = x;
double r362927 = 2.0;
double r362928 = r362926 * r362927;
double r362929 = 1.0;
double r362930 = z;
double r362931 = r362929 / r362930;
double r362932 = t;
double r362933 = r362923 - r362932;
double r362934 = r362931 / r362933;
double r362935 = r362928 * r362934;
double r362936 = 6.17105738286754e-299;
bool r362937 = r362923 <= r362936;
double r362938 = r362928 / r362930;
double r362939 = r362938 / r362933;
double r362940 = 1.9511300679708492e+111;
bool r362941 = r362923 <= r362940;
double r362942 = r362923 * r362930;
double r362943 = -r362932;
double r362944 = r362930 * r362943;
double r362945 = r362942 + r362944;
double r362946 = r362928 / r362945;
double r362947 = 1.538759424307254e+171;
bool r362948 = r362923 <= r362947;
double r362949 = r362927 * r362926;
double r362950 = r362949 / r362933;
double r362951 = r362950 / r362930;
double r362952 = r362948 ? r362939 : r362951;
double r362953 = r362941 ? r362946 : r362952;
double r362954 = r362937 ? r362939 : r362953;
double r362955 = r362925 ? r362935 : r362954;
return r362955;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.9 |
|---|---|
| Target | 2.2 |
| Herbie | 5.5 |
if y < -2.55822573448063e+208Initial program 10.7
Simplified6.9
rmApplied associate-/r*6.8
rmApplied *-un-lft-identity6.8
Applied div-inv6.9
Applied times-frac6.7
Simplified6.7
if -2.55822573448063e+208 < y < 6.17105738286754e-299 or 1.9511300679708492e+111 < y < 1.538759424307254e+171Initial program 6.4
Simplified5.6
rmApplied associate-/r*5.2
if 6.17105738286754e-299 < y < 1.9511300679708492e+111Initial program 5.3
Simplified5.0
rmApplied sub-neg5.0
Applied distribute-lft-in5.3
Simplified5.3
if 1.538759424307254e+171 < y Initial program 10.5
Simplified7.0
rmApplied associate-/r*8.4
rmApplied clear-num8.9
rmApplied associate-/r/7.5
Applied associate-/r*7.1
Simplified6.7
Final simplification5.5
herbie shell --seed 2019326
(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))))