\[\tan \left(x + \varepsilon\right) - \tan x
\]
↓
\[\begin{array}{l}
t_0 := \frac{{\sin x}^{2}}{{\cos x}^{2}}\\
t_1 := \frac{\sin x}{\cos x}\\
t_2 := 1 - \frac{\tan x}{\frac{1}{\tan \varepsilon}}\\
t_3 := \tan x + \tan \varepsilon\\
t_4 := \frac{{\sin x}^{3}}{{\cos x}^{3}}\\
\mathbf{if}\;\varepsilon \leq -5.8 \cdot 10^{-5}:\\
\;\;\;\;\frac{t_3}{t_2} - \tan x\\
\mathbf{elif}\;\varepsilon \leq 1.1 \cdot 10^{-6}:\\
\;\;\;\;\mathsf{fma}\left(t_1 + t_4, \varepsilon \cdot \varepsilon, \mathsf{fma}\left({\varepsilon}^{3}, t_0 + \left(0.3333333333333333 - \frac{\sin x}{\frac{\cos x}{t_1 \cdot -0.3333333333333333 - t_4}}\right), \varepsilon + \varepsilon \cdot t_0\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t_3, \frac{1}{t_2}, -\tan x\right)\\
\end{array}
\]
double code(double x, double eps) {
return tan((x + eps)) - tan(x);
}
↓
double code(double x, double eps) {
double t_0 = pow(sin(x), 2.0) / pow(cos(x), 2.0);
double t_1 = sin(x) / cos(x);
double t_2 = 1.0 - (tan(x) / (1.0 / tan(eps)));
double t_3 = tan(x) + tan(eps);
double t_4 = pow(sin(x), 3.0) / pow(cos(x), 3.0);
double tmp;
if (eps <= -5.8e-5) {
tmp = (t_3 / t_2) - tan(x);
} else if (eps <= 1.1e-6) {
tmp = fma((t_1 + t_4), (eps * eps), fma(pow(eps, 3.0), (t_0 + (0.3333333333333333 - (sin(x) / (cos(x) / ((t_1 * -0.3333333333333333) - t_4))))), (eps + (eps * t_0))));
} else {
tmp = fma(t_3, (1.0 / t_2), -tan(x));
}
return tmp;
}
function code(x, eps)
return Float64(tan(Float64(x + eps)) - tan(x))
end
↓
function code(x, eps)
t_0 = Float64((sin(x) ^ 2.0) / (cos(x) ^ 2.0))
t_1 = Float64(sin(x) / cos(x))
t_2 = Float64(1.0 - Float64(tan(x) / Float64(1.0 / tan(eps))))
t_3 = Float64(tan(x) + tan(eps))
t_4 = Float64((sin(x) ^ 3.0) / (cos(x) ^ 3.0))
tmp = 0.0
if (eps <= -5.8e-5)
tmp = Float64(Float64(t_3 / t_2) - tan(x));
elseif (eps <= 1.1e-6)
tmp = fma(Float64(t_1 + t_4), Float64(eps * eps), fma((eps ^ 3.0), Float64(t_0 + Float64(0.3333333333333333 - Float64(sin(x) / Float64(cos(x) / Float64(Float64(t_1 * -0.3333333333333333) - t_4))))), Float64(eps + Float64(eps * t_0))));
else
tmp = fma(t_3, Float64(1.0 / t_2), Float64(-tan(x)));
end
return tmp
end
code[x_, eps_] := N[(N[Tan[N[(x + eps), $MachinePrecision]], $MachinePrecision] - N[Tan[x], $MachinePrecision]), $MachinePrecision]
↓
code[x_, eps_] := Block[{t$95$0 = N[(N[Power[N[Sin[x], $MachinePrecision], 2.0], $MachinePrecision] / N[Power[N[Cos[x], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Sin[x], $MachinePrecision] / N[Cos[x], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(1.0 - N[(N[Tan[x], $MachinePrecision] / N[(1.0 / N[Tan[eps], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[Tan[x], $MachinePrecision] + N[Tan[eps], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(N[Power[N[Sin[x], $MachinePrecision], 3.0], $MachinePrecision] / N[Power[N[Cos[x], $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[eps, -5.8e-5], N[(N[(t$95$3 / t$95$2), $MachinePrecision] - N[Tan[x], $MachinePrecision]), $MachinePrecision], If[LessEqual[eps, 1.1e-6], N[(N[(t$95$1 + t$95$4), $MachinePrecision] * N[(eps * eps), $MachinePrecision] + N[(N[Power[eps, 3.0], $MachinePrecision] * N[(t$95$0 + N[(0.3333333333333333 - N[(N[Sin[x], $MachinePrecision] / N[(N[Cos[x], $MachinePrecision] / N[(N[(t$95$1 * -0.3333333333333333), $MachinePrecision] - t$95$4), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(eps + N[(eps * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$3 * N[(1.0 / t$95$2), $MachinePrecision] + (-N[Tan[x], $MachinePrecision])), $MachinePrecision]]]]]]]]
\tan \left(x + \varepsilon\right) - \tan x
↓
\begin{array}{l}
t_0 := \frac{{\sin x}^{2}}{{\cos x}^{2}}\\
t_1 := \frac{\sin x}{\cos x}\\
t_2 := 1 - \frac{\tan x}{\frac{1}{\tan \varepsilon}}\\
t_3 := \tan x + \tan \varepsilon\\
t_4 := \frac{{\sin x}^{3}}{{\cos x}^{3}}\\
\mathbf{if}\;\varepsilon \leq -5.8 \cdot 10^{-5}:\\
\;\;\;\;\frac{t_3}{t_2} - \tan x\\
\mathbf{elif}\;\varepsilon \leq 1.1 \cdot 10^{-6}:\\
\;\;\;\;\mathsf{fma}\left(t_1 + t_4, \varepsilon \cdot \varepsilon, \mathsf{fma}\left({\varepsilon}^{3}, t_0 + \left(0.3333333333333333 - \frac{\sin x}{\frac{\cos x}{t_1 \cdot -0.3333333333333333 - t_4}}\right), \varepsilon + \varepsilon \cdot t_0\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t_3, \frac{1}{t_2}, -\tan x\right)\\
\end{array}