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




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 11.2 |
|---|---|
| Target | 2.1 |
| Herbie | 1.0 |
Initial program 11.2
rmApplied associate-/l*_binary64_163912.1
rmApplied add-cube-cbrt_binary64_164813.2
Applied add-cube-cbrt_binary64_164812.9
Applied times-frac_binary64_164522.9
Applied associate-/r*_binary64_163901.0
Final simplification1.0
herbie shell --seed 2021076
(FPCore (x y z t)
:name "Graphics.Rendering.Chart.Plot.AreaSpots:renderAreaSpots4D from Chart-1.5.3"
:precision binary64
:herbie-target
(/ x (/ (- t z) (- y z)))
(/ (* x (- y z)) (- t z)))