\left(x \cdot 2 - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right) + \left(a \cdot 27\right) \cdot b
\begin{array}{l}
\mathbf{if}\;\left(x \cdot 2 - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right) + \left(a \cdot 27\right) \cdot b = -\infty \lor \neg \left(\left(x \cdot 2 - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right) + \left(a \cdot 27\right) \cdot b \le 3.7125537497480025 \cdot 10^{275}\right):\\
\;\;\;\;\mathsf{fma}\left(a, 27 \cdot b, x \cdot 2 - 9 \cdot \left(\left(t \cdot y\right) \cdot z\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot 2 - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right) + \left(a \cdot 27\right) \cdot b\\
\end{array}double code(double x, double y, double z, double t, double a, double b) {
return (((x * 2.0) - (((y * 9.0) * z) * t)) + ((a * 27.0) * b));
}
double code(double x, double y, double z, double t, double a, double b) {
double temp;
if ((((((x * 2.0) - (((y * 9.0) * z) * t)) + ((a * 27.0) * b)) <= -inf.0) || !((((x * 2.0) - (((y * 9.0) * z) * t)) + ((a * 27.0) * b)) <= 3.7125537497480025e+275))) {
temp = fma(a, (27.0 * b), ((x * 2.0) - (9.0 * ((t * y) * z))));
} else {
temp = (((x * 2.0) - (((y * 9.0) * z) * t)) + ((a * 27.0) * b));
}
return temp;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b
Results
| Original | 3.5 |
|---|---|
| Target | 2.6 |
| Herbie | 0.7 |
if (+ (- (* x 2.0) (* (* (* y 9.0) z) t)) (* (* a 27.0) b)) < -inf.0 or 3.7125537497480025e+275 < (+ (- (* x 2.0) (* (* (* y 9.0) z) t)) (* (* a 27.0) b)) Initial program 31.9
Simplified31.0
rmApplied *-commutative31.0
Applied associate-*l*30.5
Applied associate-*l*30.2
Simplified30.2
rmApplied *-commutative30.2
Applied associate-*r*3.0
if -inf.0 < (+ (- (* x 2.0) (* (* (* y 9.0) z) t)) (* (* a 27.0) b)) < 3.7125537497480025e+275Initial program 0.4
Final simplification0.7
herbie shell --seed 2020066 +o rules:numerics
(FPCore (x y z t a b)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, A"
:precision binary64
:herbie-target
(if (< y 7.590524218811189e-161) (+ (- (* x 2) (* (* (* y 9) z) t)) (* a (* 27 b))) (+ (- (* x 2) (* 9 (* y (* t z)))) (* (* a 27) b)))
(+ (- (* x 2) (* (* (* y 9) z) t)) (* (* a 27) b)))