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) (((double) (((double) cos(b)) / ((double) (((double) sin(b)) / ((double) cos(a)))))) - ((double) sin(a))))));
}



Bits error versus r



Bits error versus a



Bits error versus b
Results
Initial program 14.8
Simplified14.8
rmApplied cos-sum0.3
rmApplied associate-*r/0.3
Simplified0.3
Taylor expanded around inf 0.3
Simplified0.4
Final simplification0.4
herbie shell --seed 2020192
(FPCore (r a b)
:name "rsin A"
:precision binary64
(/ (* r (sin b)) (cos (+ a b))))