r \cdot \frac{\sin b}{\cos \left(a + b\right)}\frac{r}{\mathsf{fma}\left(\frac{\cos a}{\sin b}, \cos b, -\sin a\right)}double code(double r, double a, double b) {
return (r * (sin(b) / cos((a + b))));
}
double code(double r, double a, double b) {
return (r / fma((cos(a) / sin(b)), cos(b), -sin(a)));
}



Bits error versus r



Bits error versus a



Bits error versus b
Results
Initial program 15.1
rmApplied cos-sum0.3
rmApplied associate-*r/0.3
rmApplied associate-/l*0.4
Taylor expanded around inf 0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020100 +o rules:numerics
(FPCore (r a b)
:name "r*sin(b)/cos(a+b), B"
:precision binary64
(* r (/ (sin b) (cos (+ a b)))))