\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}\;t \le -2.04995563102749789 \cdot 10^{-48}:\\
\;\;\;\;\left(2 \cdot x - 9 \cdot \left(t \cdot \left(z \cdot y\right)\right)\right) + \sqrt{27} \cdot \left(\sqrt{27} \cdot \left(a \cdot b\right)\right)\\
\mathbf{elif}\;t \le 2.6373343599119812 \cdot 10^{78}:\\
\;\;\;\;\left(2 \cdot x - 9 \cdot \left(\left(t \cdot z\right) \cdot y\right)\right) + 27 \cdot \left(a \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(a, 27 \cdot b, x \cdot 2 - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right)\\
\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 VAR;
if ((t <= -2.049955631027498e-48)) {
VAR = (((2.0 * x) - (9.0 * (t * (z * y)))) + (sqrt(27.0) * (sqrt(27.0) * (a * b))));
} else {
double VAR_1;
if ((t <= 2.6373343599119812e+78)) {
VAR_1 = (((2.0 * x) - (9.0 * ((t * z) * y))) + (27.0 * (a * b)));
} else {
VAR_1 = fma(a, (27.0 * b), ((x * 2.0) - (((y * 9.0) * z) * t)));
}
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
Results
| Original | 3.8 |
|---|---|
| Target | 2.7 |
| Herbie | 0.8 |
if t < -2.049955631027498e-48Initial program 1.1
Taylor expanded around inf 1.0
Taylor expanded around 0 0.9
rmApplied add-sqr-sqrt0.9
Applied associate-*l*1.0
if -2.049955631027498e-48 < t < 2.6373343599119812e+78Initial program 5.6
Taylor expanded around inf 5.5
Taylor expanded around 0 5.4
rmApplied associate-*r*0.6
if 2.6373343599119812e+78 < t Initial program 1.2
Simplified1.2
Final simplification0.8
herbie shell --seed 2020078 +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)))