double f(double a, double b, double c, double d) {
double r12124106 = a;
double r12124107 = b;
double r12124108 = c;
double r12124109 = d;
double r12124110 = r12124108 + r12124109;
double r12124111 = r12124107 + r12124110;
double r12124112 = r12124106 + r12124111;
double r12124113 = 2.0;
double r12124114 = r12124112 * r12124113;
return r12124114;
}
double f(double a, double b, double c, double d) {
double r12124115 = 2.0;
double r12124116 = a;
double r12124117 = c;
double r12124118 = r12124116 + r12124117;
double r12124119 = b;
double r12124120 = d;
double r12124121 = r12124119 + r12124120;
double r12124122 = r12124118 + r12124121;
double r12124123 = exp(r12124122);
double r12124124 = log(r12124123);
double r12124125 = r12124119 + r12124117;
double r12124126 = r12124125 + r12124120;
double r12124127 = r12124126 + r12124116;
double r12124128 = r12124124 * r12124127;
double r12124129 = r12124128 * r12124127;
double r12124130 = cbrt(r12124129);
double r12124131 = r12124115 * r12124130;
return r12124131;
}
\left(a + \left(b + \left(c + d\right)\right)\right) \cdot 2
2 \cdot \sqrt[3]{\left(\log \left(e^{\left(a + c\right) + \left(b + d\right)}\right) \cdot \left(\left(\left(b + c\right) + d\right) + a\right)\right) \cdot \left(\left(\left(b + c\right) + d\right) + a\right)}



Bits error versus a




Bits error versus b




Bits error versus c




Bits error versus d
| Original | 3.7 |
|---|---|
| Target | 3.8 |
| Herbie | 2.7 |
Initial program 3.7
rmApplied associate-+r+2.7
rmApplied add-cbrt-cube2.9
rmApplied associate-+r+2.7
rmApplied add-log-exp2.7
Applied add-log-exp2.7
Applied sum-log2.8
Simplified2.7
Final simplification2.7
herbie shell --seed 2019101 +o rules:numerics
(FPCore (a b c d)
:name "Expression, p6"
: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))