\frac{x}{\left(y - z\right) \cdot \left(t - z\right)}\frac{\sqrt[3]{x}}{\sqrt[3]{y - z} \cdot \sqrt[3]{y - z}} \cdot \left(\frac{\sqrt[3]{x}}{\sqrt[3]{y - z}} \cdot \frac{\sqrt[3]{x}}{t - z}\right)double f(double x, double y, double z, double t) {
double r802824 = x;
double r802825 = y;
double r802826 = z;
double r802827 = r802825 - r802826;
double r802828 = t;
double r802829 = r802828 - r802826;
double r802830 = r802827 * r802829;
double r802831 = r802824 / r802830;
return r802831;
}
double f(double x, double y, double z, double t) {
double r802832 = x;
double r802833 = cbrt(r802832);
double r802834 = y;
double r802835 = z;
double r802836 = r802834 - r802835;
double r802837 = cbrt(r802836);
double r802838 = r802837 * r802837;
double r802839 = r802833 / r802838;
double r802840 = r802833 / r802837;
double r802841 = t;
double r802842 = r802841 - r802835;
double r802843 = r802833 / r802842;
double r802844 = r802840 * r802843;
double r802845 = r802839 * r802844;
return r802845;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.4 |
|---|---|
| Target | 8.0 |
| Herbie | 1.2 |
Initial program 7.4
rmApplied add-cube-cbrt7.9
Applied times-frac1.8
rmApplied add-cube-cbrt2.0
Applied times-frac1.9
Applied associate-*l*1.2
Final simplification1.2
herbie shell --seed 2020033 +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))))