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}\;\left(a + b \cdot c\right) \cdot c \le -3.1861915954316597 \cdot 10^{216}:\\
\;\;\;\;\mathsf{fma}\left(2, \mathsf{fma}\left(x, y, z \cdot t\right), 2 \cdot \left(\left(-\mathsf{fma}\left(c, b, a\right)\right) \cdot \left(c \cdot i\right)\right)\right)\\
\mathbf{elif}\;\left(a + b \cdot c\right) \cdot c \le 4.53656211072419482 \cdot 10^{219}:\\
\;\;\;\;2 \cdot \left(\left(x \cdot y + z \cdot t\right) - \left(\left(a + b \cdot c\right) \cdot c\right) \cdot i\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(2, \mathsf{fma}\left(x, y, z \cdot t\right), {\left(\sqrt{\sqrt{2}}\right)}^{3} \cdot \left(\sqrt{\sqrt{2}} \cdot \left(\left(-\mathsf{fma}\left(c, b, a\right)\right) \cdot \left(c \cdot i\right)\right)\right)\right)\\
\end{array}double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return (2.0 * (((x * y) + (z * t)) - (((a + (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 ((((a + (b * c)) * c) <= -3.18619159543166e+216)) {
VAR = fma(2.0, fma(x, y, (z * t)), (2.0 * (-fma(c, b, a) * (c * i))));
} else {
double VAR_1;
if ((((a + (b * c)) * c) <= 4.536562110724195e+219)) {
VAR_1 = (2.0 * (((x * y) + (z * t)) - (((a + (b * c)) * c) * i)));
} else {
VAR_1 = fma(2.0, fma(x, y, (z * t)), (pow(sqrt(sqrt(2.0)), 3.0) * (sqrt(sqrt(2.0)) * (-fma(c, b, a) * (c * i)))));
}
VAR = VAR_1;
}
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 | 5.9 |
|---|---|
| Target | 1.7 |
| Herbie | 1.3 |
if (* (+ a (* b c)) c) < -3.18619159543166e+216Initial program 33.2
Simplified5.4
if -3.18619159543166e+216 < (* (+ a (* b c)) c) < 4.536562110724195e+219Initial program 0.3
if 4.536562110724195e+219 < (* (+ a (* b c)) c) Initial program 32.9
Simplified6.3
rmApplied add-sqr-sqrt6.9
Applied associate-*l*6.7
rmApplied add-sqr-sqrt6.7
Applied sqrt-prod6.3
Applied associate-*l*6.4
rmApplied associate-*r*6.4
Simplified6.5
Final simplification1.3
herbie shell --seed 2020100 +o rules:numerics
(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 (- (+ (* x y) (* z t)) (* (+ a (* b c)) (* c i))))
(* 2 (- (+ (* x y) (* z t)) (* (* (+ a (* b c)) c) i))))