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




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.6 |
|---|---|
| Target | 8.2 |
| Herbie | 1.1 |
Initial program 7.6
Applied associate-/r*_binary641.9
Applied *-un-lft-identity_binary641.9
Applied add-cube-cbrt_binary642.5
Applied add-cube-cbrt_binary642.7
Applied times-frac_binary642.7
Applied times-frac_binary641.1
Simplified1.1
Final simplification1.1
herbie shell --seed 2021313
(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))))