\sin x \cdot \frac{\sinh y}{y}\sin x \cdot \mathsf{fma}\left(\frac{1}{6}, {y}^{2}, \mathsf{fma}\left(\frac{1}{120}, {y}^{4}, 1\right)\right)double code(double x, double y) {
return ((double) (((double) sin(x)) * ((double) (((double) sinh(y)) / y))));
}
double code(double x, double y) {
return ((double) (((double) sin(x)) * ((double) fma(0.16666666666666666, ((double) pow(y, 2.0)), ((double) fma(0.008333333333333333, ((double) pow(y, 4.0)), 1.0))))));
}



Bits error versus x



Bits error versus y
Results
Initial program 0.0
Taylor expanded around 0 0.6
Simplified0.6
Final simplification0.6
herbie shell --seed 2020121 +o rules:numerics
(FPCore (x y)
:name "Linear.Quaternion:$ccos from linear-1.19.1.3"
:precision binary64
(* (sin x) (/ (sinh y) y)))