double code(double a1, double a2, double b1, double b2) {
return ((double) (((double) (a1 * a2)) / ((double) (b1 * b2))));
}
double code(double a1, double a2, double b1, double b2) {
double VAR;
if ((((double) (((double) (a1 * a2)) / ((double) (b1 * b2)))) <= -inf.0)) {
VAR = ((double) (((double) (a1 / b1)) * ((double) (a2 / b2))));
} else {
double VAR_1;
if ((((double) (((double) (a1 * a2)) / ((double) (b1 * b2)))) <= -2.886617978924377e-298)) {
VAR_1 = ((double) (((double) (a1 * a2)) / ((double) (b1 * b2))));
} else {
double VAR_2;
if ((((double) (((double) (a1 * a2)) / ((double) (b1 * b2)))) <= -0.0)) {
VAR_2 = ((double) (((double) (a1 * ((double) (a2 / b2)))) / b1));
} else {
double VAR_3;
if ((((double) (((double) (a1 * a2)) / ((double) (b1 * b2)))) <= 4.446807209367256e+306)) {
VAR_3 = ((double) (((double) (a1 * a2)) / ((double) (b1 * b2))));
} else {
VAR_3 = ((double) (((double) (a1 / b1)) * ((double) (a2 / b2))));
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus a1




Bits error versus a2




Bits error versus b1




Bits error versus b2
Results
| Original | 11.5 |
|---|---|
| Target | 11.7 |
| Herbie | 2.6 |
if (/ (* a1 a2) (* b1 b2)) < -inf.0 or 4.4468072093672558e306 < (/ (* a1 a2) (* b1 b2)) Initial program 63.7
Simplified41.4
rmApplied *-un-lft-identity41.4
Applied times-frac15.4
Applied associate-*r*9.4
Simplified9.3
if -inf.0 < (/ (* a1 a2) (* b1 b2)) < -2.8866179789243771e-298 or -0.0 < (/ (* a1 a2) (* b1 b2)) < 4.4468072093672558e306Initial program 0.7
if -2.8866179789243771e-298 < (/ (* a1 a2) (* b1 b2)) < -0.0Initial program 14.0
Simplified7.7
rmApplied add-cube-cbrt7.9
Applied times-frac3.9
Applied associate-*r*3.1
Simplified3.1
rmApplied associate-*r/3.1
Applied associate-*r/4.9
Applied associate-*l/5.4
Simplified3.9
Final simplification2.6
herbie shell --seed 2020180
(FPCore (a1 a2 b1 b2)
:name "Quotient of products"
:precision binary64
:herbie-target
(* (/ a1 b1) (/ a2 b2))
(/ (* a1 a2) (* b1 b2)))