2 \cdot \left(\left(x \cdot y + z \cdot t\right) - \left(\left(a + b \cdot c\right) \cdot c\right) \cdot i\right)
\begin{array}{l}
\mathbf{if}\;b \le -1.68010719244830644 \cdot 10^{-120} \lor \neg \left(b \le 2.51425795026807799 \cdot 10^{-228}\right):\\
\;\;\;\;2 \cdot \left(\left(x \cdot y + z \cdot t\right) - \sqrt[3]{a + b \cdot c} \cdot \left(\sqrt[3]{a + b \cdot c} \cdot \left(\sqrt[3]{a + b \cdot c} \cdot \left(c \cdot i\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;2 \cdot \left(\left(\sqrt[3]{x \cdot y + z \cdot t} \cdot \sqrt[3]{x \cdot y + z \cdot t}\right) \cdot \sqrt[3]{x \cdot y + z \cdot t} - \left(\left(a + b \cdot c\right) \cdot c\right) \cdot i\right)\\
\end{array}double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return ((double) (2.0 * ((double) (((double) (((double) (x * y)) + ((double) (z * t)))) - ((double) (((double) (((double) (a + ((double) (b * c)))) * c)) * i))))));
}
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double VAR;
if (((b <= -1.6801071924483064e-120) || !(b <= 2.514257950268078e-228))) {
VAR = ((double) (2.0 * ((double) (((double) (((double) (x * y)) + ((double) (z * t)))) - ((double) (((double) cbrt(((double) (a + ((double) (b * c)))))) * ((double) (((double) cbrt(((double) (a + ((double) (b * c)))))) * ((double) (((double) cbrt(((double) (a + ((double) (b * c)))))) * ((double) (c * i))))))))))));
} else {
VAR = ((double) (2.0 * ((double) (((double) (((double) (((double) cbrt(((double) (((double) (x * y)) + ((double) (z * t)))))) * ((double) cbrt(((double) (((double) (x * y)) + ((double) (z * t)))))))) * ((double) cbrt(((double) (((double) (x * y)) + ((double) (z * t)))))))) - ((double) (((double) (((double) (a + ((double) (b * c)))) * c)) * i))))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b




Bits error versus c




Bits error versus i
Results
| Original | 6.0 |
|---|---|
| Target | 2.0 |
| Herbie | 3.1 |
if b < -1.68010719244830644e-120 or 2.51425795026807799e-228 < b Initial program 6.6
rmApplied associate-*l*2.1
rmApplied add-cube-cbrt2.5
Applied associate-*l*2.5
rmApplied associate-*l*2.5
if -1.68010719244830644e-120 < b < 2.51425795026807799e-228Initial program 4.1
rmApplied add-cube-cbrt5.0
Final simplification3.1
herbie shell --seed 2020171
(FPCore (x y z t a b c i)
:name "Diagrams.ThreeD.Shapes:frustum from diagrams-lib-1.3.0.3, A"
:precision binary64
:herbie-target
(* 2.0 (- (+ (* x y) (* z t)) (* (+ a (* b c)) (* c i))))
(* 2.0 (- (+ (* x y) (* z t)) (* (* (+ a (* b c)) c) i))))