\frac{x \cdot \left(y - z\right)}{t - z}
\begin{array}{l}
t_1 := \sqrt[3]{y - z}\\
t_2 := \sqrt[3]{t - z}\\
\left(x \cdot \frac{t_1 \cdot t_1}{t_2 \cdot t_2}\right) \cdot \frac{t_1}{t_2}
\end{array}
(FPCore (x y z t) :precision binary64 (/ (* x (- y z)) (- t z)))
(FPCore (x y z t) :precision binary64 (let* ((t_1 (cbrt (- y z))) (t_2 (cbrt (- t z)))) (* (* x (/ (* t_1 t_1) (* t_2 t_2))) (/ t_1 t_2))))
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) {
double t_1 = cbrt(y - z);
double t_2 = cbrt(t - z);
return (x * ((t_1 * t_1) / (t_2 * t_2))) * (t_1 / t_2);
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 11.9 |
|---|---|
| Target | 2.1 |
| Herbie | 1.0 |
Initial program 11.9
rmApplied *-un-lft-identity_binary6411.9
Applied times-frac_binary642.2
Simplified2.2
rmApplied add-cube-cbrt_binary643.2
Applied add-cube-cbrt_binary642.9
Applied times-frac_binary642.9
Applied associate-*r*_binary641.0
Final simplification1.0
herbie shell --seed 2021204
(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)))