\frac{x}{\left(y - z\right) \cdot \left(t - z\right)}\begin{array}{l}
\mathbf{if}\;x \le -2.304871333744858729950181693700096214694 \cdot 10^{-57}:\\
\;\;\;\;\frac{\frac{1}{\frac{t - z}{x}}}{y - z}\\
\mathbf{elif}\;x \le 2.08364774622770445821819618058878978416 \cdot 10^{72}:\\
\;\;\;\;\frac{x}{\left(y - z\right) \cdot \left(t - z\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{t - z}{\sqrt{x}}}{\frac{\sqrt{x}}{y - z}}}\\
\end{array}double f(double x, double y, double z, double t) {
double r34520053 = x;
double r34520054 = y;
double r34520055 = z;
double r34520056 = r34520054 - r34520055;
double r34520057 = t;
double r34520058 = r34520057 - r34520055;
double r34520059 = r34520056 * r34520058;
double r34520060 = r34520053 / r34520059;
return r34520060;
}
double f(double x, double y, double z, double t) {
double r34520061 = x;
double r34520062 = -2.3048713337448587e-57;
bool r34520063 = r34520061 <= r34520062;
double r34520064 = 1.0;
double r34520065 = t;
double r34520066 = z;
double r34520067 = r34520065 - r34520066;
double r34520068 = r34520067 / r34520061;
double r34520069 = r34520064 / r34520068;
double r34520070 = y;
double r34520071 = r34520070 - r34520066;
double r34520072 = r34520069 / r34520071;
double r34520073 = 2.0836477462277045e+72;
bool r34520074 = r34520061 <= r34520073;
double r34520075 = r34520071 * r34520067;
double r34520076 = r34520061 / r34520075;
double r34520077 = sqrt(r34520061);
double r34520078 = r34520067 / r34520077;
double r34520079 = r34520077 / r34520071;
double r34520080 = r34520078 / r34520079;
double r34520081 = r34520064 / r34520080;
double r34520082 = r34520074 ? r34520076 : r34520081;
double r34520083 = r34520063 ? r34520072 : r34520082;
return r34520083;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.4 |
|---|---|
| Target | 8.3 |
| Herbie | 1.7 |
if x < -2.3048713337448587e-57Initial program 12.2
rmApplied associate-/r*2.2
rmApplied clear-num3.0
rmApplied associate-/r/2.5
Applied associate-/r*1.9
if -2.3048713337448587e-57 < x < 2.0836477462277045e+72Initial program 1.5
rmApplied associate-/r*1.8
rmApplied div-inv1.9
Applied associate-/l*1.5
Simplified1.5
if 2.0836477462277045e+72 < x Initial program 17.7
rmApplied associate-/r*2.9
rmApplied clear-num3.7
rmApplied *-un-lft-identity3.7
Applied add-sqr-sqrt3.8
Applied times-frac3.9
Applied associate-/r*1.9
Simplified1.9
Final simplification1.7
herbie shell --seed 2019200
(FPCore (x y z t)
:name "Data.Random.Distribution.Triangular:triangularCDF from random-fu-0.2.6.2, B"
:herbie-target
(if (< (/ x (* (- y z) (- t z))) 0.0) (/ (/ x (- y z)) (- t z)) (* x (/ 1.0 (* (- y z) (- t z)))))
(/ x (* (- y z) (- t z))))