\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 -8.812482183219169527609542821052580809101 \cdot 10^{214} \lor \neg \left(\left(y - z\right) \cdot \left(t - z\right) \le 21364.92791914778717909939587116241455078\right):\\
\;\;\;\;\frac{\frac{x}{y - z}}{t - z}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{1}{\left(y - z\right) \cdot \left(t - z\right)}\\
\end{array}double f(double x, double y, double z, double t) {
double r1012548 = x;
double r1012549 = y;
double r1012550 = z;
double r1012551 = r1012549 - r1012550;
double r1012552 = t;
double r1012553 = r1012552 - r1012550;
double r1012554 = r1012551 * r1012553;
double r1012555 = r1012548 / r1012554;
return r1012555;
}
double f(double x, double y, double z, double t) {
double r1012556 = y;
double r1012557 = z;
double r1012558 = r1012556 - r1012557;
double r1012559 = t;
double r1012560 = r1012559 - r1012557;
double r1012561 = r1012558 * r1012560;
double r1012562 = -8.81248218321917e+214;
bool r1012563 = r1012561 <= r1012562;
double r1012564 = 21364.927919147787;
bool r1012565 = r1012561 <= r1012564;
double r1012566 = !r1012565;
bool r1012567 = r1012563 || r1012566;
double r1012568 = x;
double r1012569 = r1012568 / r1012558;
double r1012570 = r1012569 / r1012560;
double r1012571 = 1.0;
double r1012572 = r1012571 / r1012561;
double r1012573 = r1012568 * r1012572;
double r1012574 = r1012567 ? r1012570 : r1012573;
return r1012574;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.3 |
|---|---|
| Target | 8.1 |
| Herbie | 1.4 |
if (* (- y z) (- t z)) < -8.81248218321917e+214 or 21364.927919147787 < (* (- y z) (- t z)) Initial program 9.5
rmApplied associate-/r*1.1
rmApplied clear-num1.7
rmApplied *-un-lft-identity1.7
Applied *-un-lft-identity1.7
Applied times-frac1.7
Applied *-un-lft-identity1.7
Applied times-frac1.7
Applied add-cube-cbrt1.7
Applied times-frac1.7
Simplified1.7
Simplified1.1
if -8.81248218321917e+214 < (* (- y z) (- t z)) < 21364.927919147787Initial program 1.9
rmApplied div-inv2.3
Final simplification1.4
herbie shell --seed 2019350 +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))))