\left(a + \left(b + \left(c + d\right)\right)\right) \cdot 2
\mathsf{expm1}\left(\mathsf{log1p}\left(a + \left(\left(b + c\right) + d\right)\right)\right) \cdot 2double code(double a, double b, double c, double d) {
return ((double) (((double) (a + ((double) (b + ((double) (c + d)))))) * 2.0));
}
double code(double a, double b, double c, double d) {
return ((double) (((double) expm1(((double) log1p(((double) (a + ((double) (((double) (b + c)) + d)))))))) * 2.0));
}




Bits error versus a




Bits error versus b




Bits error versus c




Bits error versus d
Results
| Original | 3.6 |
|---|---|
| Target | 3.8 |
| Herbie | 3.0 |
Initial program 3.6
rmApplied associate-+r+2.7
rmApplied expm1-log1p-u3.0
Final simplification3.0
herbie shell --seed 2020121 +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))