\frac{x}{\left(y - z\right) \cdot \left(t - z\right)}\sqrt[3]{\frac{\sqrt[3]{x}}{t - z}} \cdot \left(\sqrt[3]{x} \cdot \left(\frac{\sqrt[3]{x}}{y - z} \cdot \left(\sqrt[3]{\frac{\sqrt[3]{x}}{t - z}} \cdot \left(\sqrt[3]{\sqrt[3]{\frac{\sqrt[3]{x}}{t - z}}} \cdot \left(\sqrt[3]{\sqrt[3]{\frac{\sqrt[3]{x}}{t - z}}} \cdot \sqrt[3]{\sqrt[3]{\frac{\sqrt[3]{x}}{t - z}}}\right)\right)\right)\right)\right)(FPCore (x y z t) :precision binary64 (/ x (* (- y z) (- t z))))
(FPCore (x y z t)
:precision binary64
(*
(cbrt (/ (cbrt x) (- t z)))
(*
(cbrt x)
(*
(/ (cbrt x) (- y z))
(*
(cbrt (/ (cbrt x) (- t z)))
(*
(cbrt (cbrt (/ (cbrt x) (- t z))))
(*
(cbrt (cbrt (/ (cbrt x) (- t z))))
(cbrt (cbrt (/ (cbrt x) (- t z)))))))))))double code(double x, double y, double z, double t) {
return (x / ((double) (((double) (y - z)) * ((double) (t - z)))));
}
double code(double x, double y, double z, double t) {
return ((double) (((double) cbrt((((double) cbrt(x)) / ((double) (t - z))))) * ((double) (((double) cbrt(x)) * ((double) ((((double) cbrt(x)) / ((double) (y - z))) * ((double) (((double) cbrt((((double) cbrt(x)) / ((double) (t - z))))) * ((double) (((double) cbrt(((double) cbrt((((double) cbrt(x)) / ((double) (t - z))))))) * ((double) (((double) cbrt(((double) cbrt((((double) cbrt(x)) / ((double) (t - z))))))) * ((double) cbrt(((double) cbrt((((double) cbrt(x)) / ((double) (t - z)))))))))))))))))));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.9 |
|---|---|
| Target | 8.8 |
| Herbie | 1.8 |
Initial program Error: 7.9 bits
rmApplied add-cube-cbrtError: 8.4 bits
Applied times-fracError: 1.6 bits
SimplifiedError: 1.6 bits
rmApplied add-cube-cbrtError: 1.8 bits
Applied associate-*r*Error: 1.8 bits
SimplifiedError: 1.6 bits
rmApplied add-cube-cbrtError: 1.8 bits
Final simplificationError: 1.8 bits
herbie shell --seed 2020203
(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.0 (* (- y z) (- t z)))))
(/ x (* (- y z) (- t z))))