\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}\;z \leq 1.7667806396472788 \cdot 10^{-44}:\\
\;\;\;\;\left(x \cdot 2 - \left(9 \cdot y\right) \cdot \left(z \cdot t\right)\right) + a \cdot \left(27 \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot 2 - 9 \cdot \left(t \cdot \left(z \cdot y\right)\right)\right) + b \cdot \left(a \cdot 27\right)\\
\end{array}(FPCore (x y z t a b) :precision binary64 (+ (- (* x 2.0) (* (* (* y 9.0) z) t)) (* (* a 27.0) b)))
(FPCore (x y z t a b) :precision binary64 (if (<= z 1.7667806396472788e-44) (+ (- (* x 2.0) (* (* 9.0 y) (* z t))) (* a (* 27.0 b))) (+ (- (* x 2.0) (* 9.0 (* t (* z y)))) (* b (* a 27.0)))))
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 tmp;
if (z <= 1.7667806396472788e-44) {
tmp = ((x * 2.0) - ((9.0 * y) * (z * t))) + (a * (27.0 * b));
} else {
tmp = ((x * 2.0) - (9.0 * (t * (z * y)))) + (b * (a * 27.0));
}
return tmp;
}




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.1 |
|---|---|
| Target | 3.7 |
| Herbie | 0.6 |
if z < 1.76678063964727877e-44Initial program 3.6
rmApplied associate-*l*_binary64_191153.7
rmApplied associate-*l*_binary64_191150.7
Simplified0.7
rmApplied associate-*l*_binary64_191150.6
rmApplied associate-*r*_binary64_191140.7
Simplified0.7
if 1.76678063964727877e-44 < z Initial program 0.3
Taylor expanded around 0 0.3
Final simplification0.6
herbie shell --seed 2021076
(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.0) (* (* (* y 9.0) z) t)) (* a (* 27.0 b))) (+ (- (* x 2.0) (* 9.0 (* y (* t z)))) (* (* a 27.0) b)))
(+ (- (* x 2.0) (* (* (* y 9.0) z) t)) (* (* a 27.0) b)))