\frac{x}{\left(y - z\right) \cdot \left(t - z\right)}\left(\frac{\sqrt[3]{x}}{\sqrt[3]{t - z}} \cdot \frac{\sqrt[3]{x}}{\sqrt[3]{t - z}}\right) \cdot \left(\frac{1}{\sqrt[3]{y - z} \cdot \sqrt[3]{y - z}} \cdot \frac{\frac{\sqrt[3]{x}}{\sqrt[3]{y - z}}}{\sqrt[3]{t - z}}\right)double f(double x, double y, double z, double t) {
double r41480589 = x;
double r41480590 = y;
double r41480591 = z;
double r41480592 = r41480590 - r41480591;
double r41480593 = t;
double r41480594 = r41480593 - r41480591;
double r41480595 = r41480592 * r41480594;
double r41480596 = r41480589 / r41480595;
return r41480596;
}
double f(double x, double y, double z, double t) {
double r41480597 = x;
double r41480598 = cbrt(r41480597);
double r41480599 = t;
double r41480600 = z;
double r41480601 = r41480599 - r41480600;
double r41480602 = cbrt(r41480601);
double r41480603 = r41480598 / r41480602;
double r41480604 = r41480603 * r41480603;
double r41480605 = 1.0;
double r41480606 = y;
double r41480607 = r41480606 - r41480600;
double r41480608 = cbrt(r41480607);
double r41480609 = r41480608 * r41480608;
double r41480610 = r41480605 / r41480609;
double r41480611 = r41480598 / r41480608;
double r41480612 = r41480611 / r41480602;
double r41480613 = r41480610 * r41480612;
double r41480614 = r41480604 * r41480613;
return r41480614;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.6 |
|---|---|
| Target | 8.3 |
| Herbie | 1.2 |
Initial program 7.6
rmApplied associate-/r*2.1
rmApplied add-cube-cbrt2.7
Applied *-un-lft-identity2.7
Applied add-cube-cbrt2.8
Applied times-frac2.8
Applied times-frac1.2
Simplified1.3
rmApplied *-un-lft-identity1.3
Applied cbrt-prod1.3
Applied add-cube-cbrt1.4
Applied *-un-lft-identity1.4
Applied cbrt-prod1.4
Applied times-frac1.4
Applied times-frac1.2
Simplified1.2
Final simplification1.2
herbie shell --seed 2019171
(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))))