\left(a + \left(b + \left(c + d\right)\right)\right) \cdot 2
\mathsf{expm1}\left(\sqrt[3]{{\left(\mathsf{log1p}\left(\mathsf{log1p}\left(\mathsf{expm1}\left(a + \left(\left(b + c\right) + d\right)\right)\right)\right)\right)}^{3}}\right) \cdot 2double code(double a, double b, double c, double d) {
return ((a + (b + (c + d))) * 2.0);
}
double code(double a, double b, double c, double d) {
return (expm1(cbrt(pow(log1p(log1p(expm1((a + ((b + c) + d))))), 3.0))) * 2.0);
}




Bits error versus a




Bits error versus b




Bits error versus c




Bits error versus d
Results
| Original | 3.7 |
|---|---|
| Target | 3.9 |
| Herbie | 3.3 |
Initial program 3.7
rmApplied associate-+r+2.8
rmApplied expm1-log1p-u3.2
rmApplied add-cbrt-cube3.3
Simplified3.3
rmApplied log1p-expm1-u3.3
Final simplification3.3
herbie shell --seed 2020079 +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))