\frac{x}{\left(y - z\right) \cdot \left(t - z\right)}\begin{array}{l}
\mathbf{if}\;\left(y - z\right) \cdot \left(t - z\right) \le -6.059807205138159821074577788383065783998 \cdot 10^{153} \lor \neg \left(\left(y - z\right) \cdot \left(t - z\right) \le 2.837428587682230126049873099227392657567 \cdot 10^{298}\right):\\
\;\;\;\;\frac{\frac{x}{y - z}}{t - z}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\left(y - z\right) \cdot \left(t - z\right)}\\
\end{array}double f(double x, double y, double z, double t) {
double r548820 = x;
double r548821 = y;
double r548822 = z;
double r548823 = r548821 - r548822;
double r548824 = t;
double r548825 = r548824 - r548822;
double r548826 = r548823 * r548825;
double r548827 = r548820 / r548826;
return r548827;
}
double f(double x, double y, double z, double t) {
double r548828 = y;
double r548829 = z;
double r548830 = r548828 - r548829;
double r548831 = t;
double r548832 = r548831 - r548829;
double r548833 = r548830 * r548832;
double r548834 = -6.05980720513816e+153;
bool r548835 = r548833 <= r548834;
double r548836 = 2.83742858768223e+298;
bool r548837 = r548833 <= r548836;
double r548838 = !r548837;
bool r548839 = r548835 || r548838;
double r548840 = x;
double r548841 = r548840 / r548830;
double r548842 = r548841 / r548832;
double r548843 = r548840 / r548833;
double r548844 = r548839 ? r548842 : r548843;
return r548844;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.4 |
|---|---|
| Target | 8.3 |
| Herbie | 0.8 |
if (* (- y z) (- t z)) < -6.05980720513816e+153 or 2.83742858768223e+298 < (* (- y z) (- t z)) Initial program 13.8
rmApplied associate-/r*0.3
if -6.05980720513816e+153 < (* (- y z) (- t z)) < 2.83742858768223e+298Initial program 1.2
Final simplification0.8
herbie shell --seed 2019303 +o rules:numerics
(FPCore (x y z t)
:name "Data.Random.Distribution.Triangular:triangularCDF from random-fu-0.2.6.2, B"
:precision binary64
:herbie-target
(if (< (/ x (* (- y z) (- t z))) 0.0) (/ (/ x (- y z)) (- t z)) (* x (/ 1 (* (- y z) (- t z)))))
(/ x (* (- y z) (- t z))))