double code(double r, double a, double b) {
return ((double) (((double) (r * ((double) sin(b)))) / ((double) cos(((double) (a + b))))));
}
double code(double r, double a, double b) {
return ((double) (r * ((double) (1.0 / ((double) (((double) (((double) cos(a)) * ((double) (((double) cos(b)) / ((double) sin(b)))))) - ((double) sin(a))))))));
}



Bits error versus r



Bits error versus a



Bits error versus b
Results
Initial program 15.1
rmApplied cos-sum0.3
rmApplied *-un-lft-identity0.3
Applied times-frac0.3
Simplified0.3
Simplified0.3
rmApplied clear-num0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020181
(FPCore (r a b)
:name "rsin A"
:precision binary64
(/ (* r (sin b)) (cos (+ a b))))