\left(-x \cdot \frac{1}{\tan B}\right) + \frac{1}{\sin B}\frac{1}{\sin B} \cdot \left(1 - x \cdot \cos B\right)double code(double B, double x) {
return ((double) (((double) -(((double) (x * ((double) (1.0 / ((double) tan(B)))))))) + ((double) (1.0 / ((double) sin(B))))));
}
double code(double B, double x) {
return ((double) (((double) (1.0 / ((double) sin(B)))) * ((double) (1.0 - ((double) (x * ((double) cos(B))))))));
}



Bits error versus B



Bits error versus x
Results
Initial program 0.2
Simplified0.2
Taylor expanded around inf 0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2020121 +o rules:numerics
(FPCore (B x)
:name "VandenBroeck and Keller, Equation (24)"
:precision binary64
(+ (- (* x (/ 1 (tan B)))) (/ 1 (sin B))))