\left(\left(\left(e + d\right) + c\right) + b\right) + a
\left(e + d\right) + \mathsf{expm1}\left(\mathsf{log1p}\left(c + \left(b + a\right)\right)\right)double f(double a, double b, double c, double d, double e) {
double r87955 = e;
double r87956 = d;
double r87957 = r87955 + r87956;
double r87958 = c;
double r87959 = r87957 + r87958;
double r87960 = b;
double r87961 = r87959 + r87960;
double r87962 = a;
double r87963 = r87961 + r87962;
return r87963;
}
double f(double a, double b, double c, double d, double e) {
double r87964 = e;
double r87965 = d;
double r87966 = r87964 + r87965;
double r87967 = c;
double r87968 = b;
double r87969 = a;
double r87970 = r87968 + r87969;
double r87971 = r87967 + r87970;
double r87972 = log1p(r87971);
double r87973 = expm1(r87972);
double r87974 = r87966 + r87973;
return r87974;
}




Bits error versus a




Bits error versus b




Bits error versus c




Bits error versus d




Bits error versus e
Results
| Original | 0.4 |
|---|---|
| Target | 0.2 |
| Herbie | 0.3 |
Initial program 0.4
rmApplied associate-+l+0.3
rmApplied associate-+l+0.3
rmApplied expm1-log1p-u0.3
Final simplification0.3
herbie shell --seed 2020047 +o rules:numerics
(FPCore (a b c d e)
:name "Expression 1, p15"
:precision binary64
:pre (<= 1 a 2 b 4 c 8 d 16 e 32)
:herbie-target
(+ (+ d (+ c (+ a b))) e)
(+ (+ (+ (+ e d) c) b) a))