\frac{x}{y + x}\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{x}{y + x}\right)\right)double f(double x, double y) {
double r234650 = x;
double r234651 = y;
double r234652 = r234651 + r234650;
double r234653 = r234650 / r234652;
return r234653;
}
double f(double x, double y) {
double r234654 = x;
double r234655 = y;
double r234656 = r234655 + r234654;
double r234657 = r234654 / r234656;
double r234658 = log1p(r234657);
double r234659 = expm1(r234658);
return r234659;
}



Bits error versus x



Bits error versus y
Results
Initial program 0.0
rmApplied expm1-log1p-u0.0
Final simplification0.0
herbie shell --seed 2020083 +o rules:numerics
(FPCore (x y)
:name "AI.Clustering.Hierarchical.Internal:average from clustering-0.2.1, B"
:precision binary64
(/ x (+ y x)))