x \cdot \frac{\sin y}{y}x \cdot \frac{1}{\frac{y}{\sin y}}double code(double x, double y) {
return (x * (sin(y) / y));
}
double code(double x, double y) {
return (x * (1.0 / (y / sin(y))));
}



Bits error versus x



Bits error versus y
Results
Initial program 0.1
rmApplied clear-num0.2
Final simplification0.2
herbie shell --seed 2020102 +o rules:numerics
(FPCore (x y)
:name "Linear.Quaternion:$cexp from linear-1.19.1.3"
:precision binary64
(* x (/ (sin y) y)))