r \cdot \frac{\sin b}{\cos \left(a + b\right)}\frac{r \cdot \sin b}{\cos b \cdot \cos a - \sin a \cdot \sin b}double f(double r, double a, double b) {
double r25444 = r;
double r25445 = b;
double r25446 = sin(r25445);
double r25447 = a;
double r25448 = r25447 + r25445;
double r25449 = cos(r25448);
double r25450 = r25446 / r25449;
double r25451 = r25444 * r25450;
return r25451;
}
double f(double r, double a, double b) {
double r25452 = r;
double r25453 = b;
double r25454 = sin(r25453);
double r25455 = r25452 * r25454;
double r25456 = cos(r25453);
double r25457 = a;
double r25458 = cos(r25457);
double r25459 = r25456 * r25458;
double r25460 = sin(r25457);
double r25461 = r25460 * r25454;
double r25462 = r25459 - r25461;
double r25463 = r25455 / r25462;
return r25463;
}



Bits error versus r



Bits error versus a



Bits error versus b
Results
Initial program 14.7
rmApplied cos-sum0.3
rmApplied add-cbrt-cube0.4
Applied add-cbrt-cube0.4
Applied cbrt-unprod0.4
Simplified0.4
rmApplied *-un-lft-identity0.4
Applied associate-*l*0.4
Simplified0.3
Final simplification0.3
herbie shell --seed 2019325
(FPCore (r a b)
:name "r*sin(b)/cos(a+b), B"
:precision binary64
(* r (/ (sin b) (cos (+ a b)))))