\frac{1 - \tan x \cdot \tan x}{1 + \tan x \cdot \tan x}\mathsf{expm1}\left(\mathsf{log1p}\left(\frac{\sqrt{1} + \tan x}{\frac{1 + \tan x \cdot \tan x}{\sqrt{1} - \tan x}}\right)\right)double code(double x) {
return ((1.0 - (tan(x) * tan(x))) / (1.0 + (tan(x) * tan(x))));
}
double code(double x) {
return expm1(log1p(((sqrt(1.0) + tan(x)) / ((1.0 + (tan(x) * tan(x))) / (sqrt(1.0) - tan(x))))));
}



Bits error versus x
Results
Initial program 0.3
rmApplied expm1-log1p-u0.4
rmApplied add-sqr-sqrt0.4
Applied difference-of-squares0.4
Applied associate-/l*0.4
Final simplification0.4
herbie shell --seed 2020053 +o rules:numerics
(FPCore (x)
:name "Trigonometry B"
:precision binary64
(/ (- 1 (* (tan x) (tan x))) (+ 1 (* (tan x) (tan x)))))