\left(a + \left(b + \left(c + d\right)\right)\right) \cdot 2
\left(a + \left(\left(b + c\right) + d\right)\right) \cdot 2
double f(double a, double b, double c, double d) {
double r99207 = a;
double r99208 = b;
double r99209 = c;
double r99210 = d;
double r99211 = r99209 + r99210;
double r99212 = r99208 + r99211;
double r99213 = r99207 + r99212;
double r99214 = 2.0;
double r99215 = r99213 * r99214;
return r99215;
}
double f(double a, double b, double c, double d) {
double r99216 = a;
double r99217 = b;
double r99218 = c;
double r99219 = r99217 + r99218;
double r99220 = d;
double r99221 = r99219 + r99220;
double r99222 = r99216 + r99221;
double r99223 = 2.0;
double r99224 = r99222 * r99223;
return r99224;
}




Bits error versus a




Bits error versus b




Bits error versus c




Bits error versus d
Results
| Original | 3.7 |
|---|---|
| Target | 3.8 |
| Herbie | 2.8 |
Initial program 3.7
rmApplied associate-+r+2.8
rmApplied log1p-expm1-u2.8
rmApplied log1p-expm12.8
Final simplification2.8
herbie shell --seed 2020001 +o rules:numerics
(FPCore (a b c d)
:name "Expression, p6"
:precision binary64
:pre (and (<= -14 a -13) (<= -3 b -2) (<= 3 c 3.5) (<= 12.5 d 13.5))
:herbie-target
(+ (* (+ a b) 2) (* (+ c d) 2))
(* (+ a (+ b (+ c d))) 2))