(FPCore (x) :precision binary64 (/ (- 1.0 (* (tan x) (tan x))) (+ 1.0 (* (tan x) (tan x)))))
(FPCore (x) :precision binary64 (/ (fma (tan x) (tan x) -1.0) (- (fma (tan x) (tan x) 1.0))))
double code(double x) {
return (1.0 - (tan(x) * tan(x))) / (1.0 + (tan(x) * tan(x)));
}
double code(double x) {
return fma(tan(x), tan(x), -1.0) / -fma(tan(x), tan(x), 1.0);
}
function code(x) return Float64(Float64(1.0 - Float64(tan(x) * tan(x))) / Float64(1.0 + Float64(tan(x) * tan(x)))) end
function code(x) return Float64(fma(tan(x), tan(x), -1.0) / Float64(-fma(tan(x), tan(x), 1.0))) end
code[x_] := N[(N[(1.0 - N[(N[Tan[x], $MachinePrecision] * N[Tan[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(N[Tan[x], $MachinePrecision] * N[Tan[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
code[x_] := N[(N[(N[Tan[x], $MachinePrecision] * N[Tan[x], $MachinePrecision] + -1.0), $MachinePrecision] / (-N[(N[Tan[x], $MachinePrecision] * N[Tan[x], $MachinePrecision] + 1.0), $MachinePrecision])), $MachinePrecision]
\frac{1 - \tan x \cdot \tan x}{1 + \tan x \cdot \tan x}
\frac{\mathsf{fma}\left(\tan x, \tan x, -1\right)}{-\mathsf{fma}\left(\tan x, \tan x, 1\right)}



Bits error versus x
Initial program 0.3
Applied frac-2neg_binary640.3
Simplified0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2022131
(FPCore (x)
:name "Trigonometry B"
:precision binary64
(/ (- 1.0 (* (tan x) (tan x))) (+ 1.0 (* (tan x) (tan x)))))