r \cdot \frac{\sin b}{\cos \left(a + b\right)}\frac{r}{\mathsf{fma}\left(\cos b \cdot \cos a, \frac{1}{\sin b}, -\sin a\right)}double f(double r, double a, double b) {
double r25316 = r;
double r25317 = b;
double r25318 = sin(r25317);
double r25319 = a;
double r25320 = r25319 + r25317;
double r25321 = cos(r25320);
double r25322 = r25318 / r25321;
double r25323 = r25316 * r25322;
return r25323;
}
double f(double r, double a, double b) {
double r25324 = r;
double r25325 = b;
double r25326 = cos(r25325);
double r25327 = a;
double r25328 = cos(r25327);
double r25329 = r25326 * r25328;
double r25330 = 1.0;
double r25331 = sin(r25325);
double r25332 = r25330 / r25331;
double r25333 = sin(r25327);
double r25334 = -r25333;
double r25335 = fma(r25329, r25332, r25334);
double r25336 = r25324 / r25335;
return r25336;
}



Bits error versus r



Bits error versus a



Bits error versus b
Initial program 14.7
rmApplied cos-sum0.3
rmApplied associate-*r/0.3
rmApplied associate-/l*0.4
Simplified0.4
rmApplied div-inv0.4
Applied fma-neg0.4
Final simplification0.4
herbie shell --seed 2019323 +o rules:numerics
(FPCore (r a b)
:name "r*sin(b)/cos(a+b), B"
:precision binary64
(* r (/ (sin b) (cos (+ a b)))))