\frac{x}{\left(y - z\right) \cdot \left(t - z\right)}\frac{\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{\sqrt[3]{y - z} \cdot \sqrt[3]{y - z}}}{\sqrt[3]{t - z} \cdot \sqrt[3]{t - z}} \cdot \frac{\frac{\sqrt[3]{x}}{\sqrt[3]{y - z}}}{\sqrt[3]{t - z}}double f(double x, double y, double z, double t) {
double r785593 = x;
double r785594 = y;
double r785595 = z;
double r785596 = r785594 - r785595;
double r785597 = t;
double r785598 = r785597 - r785595;
double r785599 = r785596 * r785598;
double r785600 = r785593 / r785599;
return r785600;
}
double f(double x, double y, double z, double t) {
double r785601 = x;
double r785602 = cbrt(r785601);
double r785603 = r785602 * r785602;
double r785604 = y;
double r785605 = z;
double r785606 = r785604 - r785605;
double r785607 = cbrt(r785606);
double r785608 = r785607 * r785607;
double r785609 = r785603 / r785608;
double r785610 = t;
double r785611 = r785610 - r785605;
double r785612 = cbrt(r785611);
double r785613 = r785612 * r785612;
double r785614 = r785609 / r785613;
double r785615 = r785602 / r785607;
double r785616 = r785615 / r785612;
double r785617 = r785614 * r785616;
return r785617;
}




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.3 |
Initial program 7.3
rmApplied associate-/r*2.1
rmApplied add-cube-cbrt2.7
Applied add-cube-cbrt3.0
Applied add-cube-cbrt3.1
Applied times-frac3.1
Applied times-frac1.3
Final simplification1.3
herbie shell --seed 2020062 +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))))