\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 r810847 = x;
double r810848 = y;
double r810849 = z;
double r810850 = r810848 - r810849;
double r810851 = t;
double r810852 = r810851 - r810849;
double r810853 = r810850 * r810852;
double r810854 = r810847 / r810853;
return r810854;
}
double f(double x, double y, double z, double t) {
double r810855 = x;
double r810856 = cbrt(r810855);
double r810857 = r810856 * r810856;
double r810858 = y;
double r810859 = z;
double r810860 = r810858 - r810859;
double r810861 = cbrt(r810860);
double r810862 = r810861 * r810861;
double r810863 = r810857 / r810862;
double r810864 = t;
double r810865 = r810864 - r810859;
double r810866 = cbrt(r810865);
double r810867 = r810866 * r810866;
double r810868 = r810863 / r810867;
double r810869 = r810856 / r810861;
double r810870 = r810869 / r810866;
double r810871 = r810868 * r810870;
return r810871;
}




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
(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))))