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



Bits error versus x



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