\frac{x \cdot y - \left(z \cdot 9\right) \cdot t}{a \cdot 2}\begin{array}{l}
\mathbf{if}\;x \cdot y = -inf.0:\\
\;\;\;\;0.5 \cdot \left(x \cdot \frac{y}{a}\right) - 4.5 \cdot \frac{t \cdot z}{a}\\
\mathbf{elif}\;x \cdot y \le 7.13360614779853435 \cdot 10^{-270}:\\
\;\;\;\;0.5 \cdot \frac{x \cdot y}{a} - \left(4.5 \cdot \frac{t}{\sqrt[3]{a} \cdot \sqrt[3]{a}}\right) \cdot \frac{z}{\sqrt[3]{a}}\\
\mathbf{elif}\;x \cdot y \le 3.18111439552301917 \cdot 10^{-7}:\\
\;\;\;\;\frac{x \cdot y - z \cdot \left(9 \cdot t\right)}{a \cdot 2}\\
\mathbf{elif}\;x \cdot y \le 4.83469013735259175 \cdot 10^{119}:\\
\;\;\;\;0.5 \cdot \frac{x \cdot y}{a} - 4.5 \cdot \frac{t}{\frac{a}{z}}\\
\mathbf{else}:\\
\;\;\;\;0.5 \cdot \left(x \cdot \frac{y}{a}\right) - 4.5 \cdot \frac{t \cdot z}{a}\\
\end{array}double code(double x, double y, double z, double t, double a) {
return ((double) (((double) (((double) (x * y)) - ((double) (((double) (z * 9.0)) * t)))) / ((double) (a * 2.0))));
}
double code(double x, double y, double z, double t, double a) {
double VAR;
if ((((double) (x * y)) <= -inf.0)) {
VAR = ((double) (((double) (0.5 * ((double) (x * ((double) (y / a)))))) - ((double) (4.5 * ((double) (((double) (t * z)) / a))))));
} else {
double VAR_1;
if ((((double) (x * y)) <= 7.133606147798534e-270)) {
VAR_1 = ((double) (((double) (0.5 * ((double) (((double) (x * y)) / a)))) - ((double) (((double) (4.5 * ((double) (t / ((double) (((double) cbrt(a)) * ((double) cbrt(a)))))))) * ((double) (z / ((double) cbrt(a))))))));
} else {
double VAR_2;
if ((((double) (x * y)) <= 3.181114395523019e-07)) {
VAR_2 = ((double) (((double) (((double) (x * y)) - ((double) (z * ((double) (9.0 * t)))))) / ((double) (a * 2.0))));
} else {
double VAR_3;
if ((((double) (x * y)) <= 4.834690137352592e+119)) {
VAR_3 = ((double) (((double) (0.5 * ((double) (((double) (x * y)) / a)))) - ((double) (4.5 * ((double) (t / ((double) (a / z))))))));
} else {
VAR_3 = ((double) (((double) (0.5 * ((double) (x * ((double) (y / a)))))) - ((double) (4.5 * ((double) (((double) (t * z)) / a))))));
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
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
Results
| Original | 8.2 |
|---|---|
| Target | 5.8 |
| Herbie | 4.8 |
if (* x y) < -inf.0 or 4.834690137352592e+119 < (* x y) Initial program 28.5
Taylor expanded around 0 28.3
rmApplied *-un-lft-identity28.3
Applied times-frac6.9
Simplified6.9
if -inf.0 < (* x y) < 7.133606147798534e-270Initial program 5.0
Taylor expanded around 0 4.9
rmApplied add-cube-cbrt5.4
Applied times-frac5.0
Applied associate-*r*5.0
if 7.133606147798534e-270 < (* x y) < 3.181114395523019e-07Initial program 3.6
rmApplied associate-*l*3.7
if 3.181114395523019e-07 < (* x y) < 4.834690137352592e+119Initial program 5.6
Taylor expanded around 0 5.4
rmApplied associate-/l*3.0
Final simplification4.8
herbie shell --seed 2020126
(FPCore (x y z t a)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, I"
:precision binary64
:herbie-target
(if (< a -2.090464557976709e+86) (- (* 0.5 (/ (* y x) a)) (* 4.5 (/ t (/ a z)))) (if (< a 2.144030707833976e+99) (/ (- (* x y) (* z (* 9.0 t))) (* a 2.0)) (- (* (/ y a) (* x 0.5)) (* (/ t a) (* z 4.5)))))
(/ (- (* x y) (* (* z 9.0) t)) (* a 2.0)))