\cosh x \cdot \frac{\sin y}{y}\cosh x \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\frac{\sin y}{y}\right)\right)double f(double x, double y) {
double r445755 = x;
double r445756 = cosh(r445755);
double r445757 = y;
double r445758 = sin(r445757);
double r445759 = r445758 / r445757;
double r445760 = r445756 * r445759;
return r445760;
}
double f(double x, double y) {
double r445761 = x;
double r445762 = cosh(r445761);
double r445763 = y;
double r445764 = sin(r445763);
double r445765 = r445764 / r445763;
double r445766 = expm1(r445765);
double r445767 = log1p(r445766);
double r445768 = r445762 * r445767;
return r445768;
}




Bits error versus x




Bits error versus y
Results
| Original | 0.1 |
|---|---|
| Target | 0.1 |
| Herbie | 0.1 |
Initial program 0.1
rmApplied log1p-expm1-u0.1
Final simplification0.1
herbie shell --seed 2019323 +o rules:numerics
(FPCore (x y)
:name "Linear.Quaternion:$csinh from linear-1.19.1.3"
:precision binary64
:herbie-target
(/ (* (cosh x) (sin y)) y)
(* (cosh x) (/ (sin y) y)))