\frac{x \cdot 100}{x + y}\mathsf{log1p}\left(\mathsf{expm1}\left(\frac{x}{y + x} \cdot 100\right)\right)double f(double x, double y) {
double r728626 = x;
double r728627 = 100.0;
double r728628 = r728626 * r728627;
double r728629 = y;
double r728630 = r728626 + r728629;
double r728631 = r728628 / r728630;
return r728631;
}
double f(double x, double y) {
double r728632 = x;
double r728633 = y;
double r728634 = r728633 + r728632;
double r728635 = r728632 / r728634;
double r728636 = 100.0;
double r728637 = r728635 * r728636;
double r728638 = expm1(r728637);
double r728639 = log1p(r728638);
return r728639;
}




Bits error versus x




Bits error versus y
Results
| Original | 0.4 |
|---|---|
| Target | 0.2 |
| Herbie | 0.2 |
Initial program 0.4
rmApplied log1p-expm1-u0.5
Simplified0.2
Final simplification0.2
herbie shell --seed 2019196 +o rules:numerics
(FPCore (x y)
:name "Development.Shake.Progress:message from shake-0.15.5"
:herbie-target
(* (/ x 1.0) (/ 100.0 (+ x y)))
(/ (* x 100.0) (+ x y)))