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)) <= -2.9710752032757204e+75)) {
VAR = ((double) (((double) (y * ((double) (x / ((double) (a * 2.0)))))) - ((double) (((double) (9.0 * ((double) (t / a)))) * ((double) (z / 2.0))))));
} else {
double VAR_1;
if ((((double) (x * y)) <= -1.441855701961825e-72)) {
VAR_1 = ((double) (((double) (((double) (x * y)) - ((double) (9.0 * ((double) (t * z)))))) / ((double) (a * 2.0))));
} else {
double VAR_2;
if ((((double) (x * y)) <= 3.19291229011511e-90)) {
VAR_2 = ((double) (((double) (y * ((double) (x / ((double) (a * 2.0)))))) - ((double) (9.0 * ((double) (t * ((double) (z / ((double) (a * 2.0))))))))));
} else {
double VAR_3;
if ((((double) (x * y)) <= 4.466868290700456e+162)) {
VAR_3 = ((double) (((double) (((double) (x * y)) - ((double) (9.0 * ((double) (t * z)))))) / ((double) (a * 2.0))));
} else {
VAR_3 = ((double) (((double) (y * ((double) (x / ((double) (a * 2.0)))))) - ((double) (((double) (9.0 * ((double) (t / a)))) * ((double) (z / 2.0))))));
}
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 | 7.7 |
|---|---|
| Target | 5.7 |
| Herbie | 5.3 |
if (* x y) < -2.9710752032757204e75 or 4.4668682907004558e162 < (* x y) Initial program 18.6
Simplified18.6
rmApplied div-sub18.6
Simplified7.9
Simplified4.3
rmApplied *-un-lft-identity4.3
Applied times-frac4.3
Applied associate-*r*4.2
Simplified4.2
if -2.9710752032757204e75 < (* x y) < -1.44185570196182512e-72 or 3.1929122901151098e-90 < (* x y) < 4.4668682907004558e162Initial program 3.8
Simplified3.7
rmApplied div-sub3.8
Simplified12.0
Simplified11.8
rmApplied associate-*r/12.0
Applied associate-*r/3.8
Applied sub-div3.7
Simplified3.8
if -1.44185570196182512e-72 < (* x y) < 3.1929122901151098e-90Initial program 5.1
Simplified5.1
rmApplied div-sub5.1
Simplified5.8
Simplified7.1
rmApplied associate-*l*7.1
Final simplification5.3
herbie shell --seed 2020190
(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)))