r \cdot \frac{\sin b}{\cos \left(a + b\right)}r \cdot \frac{\sin b}{\frac{\mathsf{fma}\left(\cos a, \cos b, \sin a \cdot \sin b\right) \cdot \left(\cos a \cdot \cos b - \sin a \cdot \sin b\right)}{\mathsf{fma}\left(\cos a, \cos b, \sin a \cdot \sin b\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 * (sin(b) / ((fma(cos(a), cos(b), (sin(a) * sin(b))) * ((cos(a) * cos(b)) - (sin(a) * sin(b)))) / fma(cos(a), cos(b), (sin(a) * sin(b))))));
}



Bits error versus r



Bits error versus a



Bits error versus b
Results
Initial program 14.9
rmApplied cos-sum0.3
rmApplied flip--0.4
Simplified0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020078 +o rules:numerics
(FPCore (r a b)
:name "r*sin(b)/cos(a+b), B"
:precision binary64
(* r (/ (sin b) (cos (+ a b)))))