\cosh x \cdot \frac{\sin y}{y}\mathsf{expm1}\left(\mathsf{log1p}\left(\cosh x \cdot \frac{\sin y}{y}\right)\right)double code(double x, double y) {
return (cosh(x) * (sin(y) / y));
}
double code(double x, double y) {
return expm1(log1p((cosh(x) * (sin(y) / y))));
}




Bits error versus x




Bits error versus y
Results
| Original | 0.1 |
|---|---|
| Target | 0.1 |
| Herbie | 0.2 |
Initial program 0.1
rmApplied expm1-log1p-u0.2
Final simplification0.2
herbie shell --seed 2020057 +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)))