\cos x \cdot \frac{\sinh y}{y}\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{\cos x \cdot \sinh y}{y}\right)\right)double f(double x, double y) {
double r129438 = x;
double r129439 = cos(r129438);
double r129440 = y;
double r129441 = sinh(r129440);
double r129442 = r129441 / r129440;
double r129443 = r129439 * r129442;
return r129443;
}
double f(double x, double y) {
double r129444 = x;
double r129445 = cos(r129444);
double r129446 = y;
double r129447 = sinh(r129446);
double r129448 = r129445 * r129447;
double r129449 = r129448 / r129446;
double r129450 = log1p(r129449);
double r129451 = expm1(r129450);
return r129451;
}



Bits error versus x



Bits error versus y
Results
Initial program 0.0
rmApplied expm1-log1p-u0.3
rmApplied associate-*r/0.3
Final simplification0.3
herbie shell --seed 2020043 +o rules:numerics
(FPCore (x y)
:name "Linear.Quaternion:$csin from linear-1.19.1.3"
:precision binary64
(* (cos x) (/ (sinh y) y)))