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) (b1 * b2)) <= -2.223137056513538e+153)) {
VAR = ((double) (((double) (a2 * ((double) (a1 / b2)))) / b1));
} else {
double VAR_1;
if ((((double) (b1 * b2)) <= -1.5147515431627668e-177)) {
VAR_1 = ((double) (a1 * ((double) (a2 / ((double) (b1 * b2))))));
} else {
double VAR_2;
if ((((double) (b1 * b2)) <= 1.1272910860476818e-228)) {
VAR_2 = ((double) (((double) (a2 * ((double) (a1 / b1)))) * ((double) (1.0 / b2))));
} else {
double VAR_3;
if ((((double) (b1 * b2)) <= 1.3195895724947344e+60)) {
VAR_3 = ((double) (a1 * ((double) (a2 / ((double) (b1 * b2))))));
} else {
double VAR_4;
if ((((double) (b1 * b2)) <= 1.45640378524597e+229)) {
VAR_4 = ((double) (1.0 / ((double) (b2 * ((double) (b1 / ((double) (a2 * a1))))))));
} else {
VAR_4 = ((double) (a1 / ((double) (b1 * ((double) (b2 / a2))))));
}
VAR_3 = VAR_4;
}
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.2 |
|---|---|
| Target | 11.5 |
| Herbie | 6.3 |
if (* b1 b2) < -2.22313705651353812e153Initial program 14.2
rmApplied associate-/r*8.7
Simplified6.1
rmApplied div-inv6.1
rmApplied associate-*r/8.7
Applied associate-*l/8.5
Simplified5.8
if -2.22313705651353812e153 < (* b1 b2) < -1.51475154316276679e-177 or 1.1272910860476818e-228 < (* b1 b2) < 1.31958957249473445e60Initial program 3.9
Simplified3.7
if -1.51475154316276679e-177 < (* b1 b2) < 1.1272910860476818e-228Initial program 33.5
rmApplied associate-/r*18.1
Simplified9.5
rmApplied div-inv9.6
if 1.31958957249473445e60 < (* b1 b2) < 1.45640378524596996e229Initial program 6.6
rmApplied clear-num7.1
Simplified11.7
if 1.45640378524596996e229 < (* b1 b2) Initial program 15.4
rmApplied associate-/l*15.4
Simplified7.0
Final simplification6.3
herbie shell --seed 2020192
(FPCore (a1 a2 b1 b2)
:name "Quotient of products"
:precision binary64
:herbie-target
(* (/ a1 b1) (/ a2 b2))
(/ (* a1 a2) (* b1 b2)))