\frac{r \cdot \sin b}{\cos \left(a + b\right)}r \cdot \frac{\sin b}{\mathsf{fma}\left(\left(\cos a\right), \left(\cos b\right), \left(\left(-\sin b\right) \cdot \sin a\right)\right)}double f(double r, double a, double b) {
double r1195026 = r;
double r1195027 = b;
double r1195028 = sin(r1195027);
double r1195029 = r1195026 * r1195028;
double r1195030 = a;
double r1195031 = r1195030 + r1195027;
double r1195032 = cos(r1195031);
double r1195033 = r1195029 / r1195032;
return r1195033;
}
double f(double r, double a, double b) {
double r1195034 = r;
double r1195035 = b;
double r1195036 = sin(r1195035);
double r1195037 = a;
double r1195038 = cos(r1195037);
double r1195039 = cos(r1195035);
double r1195040 = -r1195036;
double r1195041 = sin(r1195037);
double r1195042 = r1195040 * r1195041;
double r1195043 = fma(r1195038, r1195039, r1195042);
double r1195044 = r1195036 / r1195043;
double r1195045 = r1195034 * r1195044;
return r1195045;
}



Bits error versus r



Bits error versus a



Bits error versus b
Initial program 15.4
rmApplied cos-sum0.4
rmApplied *-un-lft-identity0.4
Applied times-frac0.3
Simplified0.3
rmApplied fma-neg0.3
Final simplification0.3
herbie shell --seed 2019124 +o rules:numerics
(FPCore (r a b)
:name "r*sin(b)/cos(a+b), A"
(/ (* r (sin b)) (cos (+ a b))))