\frac{x \cdot \left(y - z\right)}{t - z}\left(x \cdot \frac{\sqrt[3]{y - z} \cdot \sqrt[3]{y - z}}{\sqrt[3]{t - z} \cdot \sqrt[3]{t - z}}\right) \cdot \frac{\sqrt[3]{y - z}}{\sqrt[3]{t - z}}(FPCore (x y z t) :precision binary64 (/ (* x (- y z)) (- t z)))
(FPCore (x y z t) :precision binary64 (* (* x (/ (* (cbrt (- y z)) (cbrt (- y z))) (* (cbrt (- t z)) (cbrt (- t z))))) (/ (cbrt (- y z)) (cbrt (- t 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(y - z) * cbrt(y - z)) / (cbrt(t - z) * cbrt(t - z)))) * (cbrt(y - z) / cbrt(t - z));
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 11.7 |
|---|---|
| Target | 2.2 |
| Herbie | 1.0 |
Initial program 11.7
rmApplied *-un-lft-identity_binary64_1610511.7
Applied times-frac_binary64_161112.2
Simplified2.2
rmApplied add-cube-cbrt_binary64_161403.2
Applied add-cube-cbrt_binary64_161402.9
Applied times-frac_binary64_161112.9
Applied associate-*r*_binary64_160451.0
Final simplification1.0
herbie shell --seed 2021098
(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)))