\[\tan \left(x + \varepsilon\right) - \tan x
\]
↓
\[\begin{array}{l}
t_0 := {\cos x}^{3}\\
t_1 := {\sin x}^{3}\\
t_2 := \tan x + \tan \varepsilon\\
t_3 := \frac{{\sin x}^{4}}{{\cos x}^{4}}\\
t_4 := {\cos x}^{2}\\
t_5 := \mathsf{fma}\left(\tan x, \tan \varepsilon, -1\right)\\
t_6 := -\tan x\\
t_7 := {\sin x}^{2}\\
t_8 := \frac{-0.3333333333333333}{\frac{t_4}{t_7}}\\
\mathbf{if}\;\varepsilon \leq -0.046:\\
\;\;\;\;\mathsf{fma}\left(t_2, \frac{-1}{t_5}, t_6\right)\\
\mathbf{elif}\;\varepsilon \leq 0.000235:\\
\;\;\;\;\frac{\tan \varepsilon}{1 - \tan x \cdot \tan \varepsilon} + \left(\left(\frac{\varepsilon}{t_4} \cdot t_7 + \left(\frac{\varepsilon \cdot \varepsilon}{\frac{t_0}{t_1}} + {\varepsilon}^{3} \cdot \left(t_3 - t_8\right)\right)\right) - \mathsf{fma}\left(-0.3333333333333333, \frac{t_1}{t_0}, \left(t_8 - t_3\right) \cdot \frac{\sin x}{\cos x}\right) \cdot {\varepsilon}^{4}\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t_2, \frac{-1}{\log \left(e^{t_5}\right)}, t_6\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(cos(x), 3.0);
double t_1 = pow(sin(x), 3.0);
double t_2 = tan(x) + tan(eps);
double t_3 = pow(sin(x), 4.0) / pow(cos(x), 4.0);
double t_4 = pow(cos(x), 2.0);
double t_5 = fma(tan(x), tan(eps), -1.0);
double t_6 = -tan(x);
double t_7 = pow(sin(x), 2.0);
double t_8 = -0.3333333333333333 / (t_4 / t_7);
double tmp;
if (eps <= -0.046) {
tmp = fma(t_2, (-1.0 / t_5), t_6);
} else if (eps <= 0.000235) {
tmp = (tan(eps) / (1.0 - (tan(x) * tan(eps)))) + ((((eps / t_4) * t_7) + (((eps * eps) / (t_0 / t_1)) + (pow(eps, 3.0) * (t_3 - t_8)))) - (fma(-0.3333333333333333, (t_1 / t_0), ((t_8 - t_3) * (sin(x) / cos(x)))) * pow(eps, 4.0)));
} else {
tmp = fma(t_2, (-1.0 / log(exp(t_5))), t_6);
}
return tmp;
}
function code(x, eps)
return Float64(tan(Float64(x + eps)) - tan(x))
end
↓
function code(x, eps)
t_0 = cos(x) ^ 3.0
t_1 = sin(x) ^ 3.0
t_2 = Float64(tan(x) + tan(eps))
t_3 = Float64((sin(x) ^ 4.0) / (cos(x) ^ 4.0))
t_4 = cos(x) ^ 2.0
t_5 = fma(tan(x), tan(eps), -1.0)
t_6 = Float64(-tan(x))
t_7 = sin(x) ^ 2.0
t_8 = Float64(-0.3333333333333333 / Float64(t_4 / t_7))
tmp = 0.0
if (eps <= -0.046)
tmp = fma(t_2, Float64(-1.0 / t_5), t_6);
elseif (eps <= 0.000235)
tmp = Float64(Float64(tan(eps) / Float64(1.0 - Float64(tan(x) * tan(eps)))) + Float64(Float64(Float64(Float64(eps / t_4) * t_7) + Float64(Float64(Float64(eps * eps) / Float64(t_0 / t_1)) + Float64((eps ^ 3.0) * Float64(t_3 - t_8)))) - Float64(fma(-0.3333333333333333, Float64(t_1 / t_0), Float64(Float64(t_8 - t_3) * Float64(sin(x) / cos(x)))) * (eps ^ 4.0))));
else
tmp = fma(t_2, Float64(-1.0 / log(exp(t_5))), t_6);
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[Power[N[Cos[x], $MachinePrecision], 3.0], $MachinePrecision]}, Block[{t$95$1 = N[Power[N[Sin[x], $MachinePrecision], 3.0], $MachinePrecision]}, Block[{t$95$2 = N[(N[Tan[x], $MachinePrecision] + N[Tan[eps], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[Power[N[Sin[x], $MachinePrecision], 4.0], $MachinePrecision] / N[Power[N[Cos[x], $MachinePrecision], 4.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[Power[N[Cos[x], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$5 = N[(N[Tan[x], $MachinePrecision] * N[Tan[eps], $MachinePrecision] + -1.0), $MachinePrecision]}, Block[{t$95$6 = (-N[Tan[x], $MachinePrecision])}, Block[{t$95$7 = N[Power[N[Sin[x], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$8 = N[(-0.3333333333333333 / N[(t$95$4 / t$95$7), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[eps, -0.046], N[(t$95$2 * N[(-1.0 / t$95$5), $MachinePrecision] + t$95$6), $MachinePrecision], If[LessEqual[eps, 0.000235], N[(N[(N[Tan[eps], $MachinePrecision] / N[(1.0 - N[(N[Tan[x], $MachinePrecision] * N[Tan[eps], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(N[(eps / t$95$4), $MachinePrecision] * t$95$7), $MachinePrecision] + N[(N[(N[(eps * eps), $MachinePrecision] / N[(t$95$0 / t$95$1), $MachinePrecision]), $MachinePrecision] + N[(N[Power[eps, 3.0], $MachinePrecision] * N[(t$95$3 - t$95$8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(-0.3333333333333333 * N[(t$95$1 / t$95$0), $MachinePrecision] + N[(N[(t$95$8 - t$95$3), $MachinePrecision] * N[(N[Sin[x], $MachinePrecision] / N[Cos[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[Power[eps, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$2 * N[(-1.0 / N[Log[N[Exp[t$95$5], $MachinePrecision]], $MachinePrecision]), $MachinePrecision] + t$95$6), $MachinePrecision]]]]]]]]]]]]
\tan \left(x + \varepsilon\right) - \tan x
↓
\begin{array}{l}
t_0 := {\cos x}^{3}\\
t_1 := {\sin x}^{3}\\
t_2 := \tan x + \tan \varepsilon\\
t_3 := \frac{{\sin x}^{4}}{{\cos x}^{4}}\\
t_4 := {\cos x}^{2}\\
t_5 := \mathsf{fma}\left(\tan x, \tan \varepsilon, -1\right)\\
t_6 := -\tan x\\
t_7 := {\sin x}^{2}\\
t_8 := \frac{-0.3333333333333333}{\frac{t_4}{t_7}}\\
\mathbf{if}\;\varepsilon \leq -0.046:\\
\;\;\;\;\mathsf{fma}\left(t_2, \frac{-1}{t_5}, t_6\right)\\
\mathbf{elif}\;\varepsilon \leq 0.000235:\\
\;\;\;\;\frac{\tan \varepsilon}{1 - \tan x \cdot \tan \varepsilon} + \left(\left(\frac{\varepsilon}{t_4} \cdot t_7 + \left(\frac{\varepsilon \cdot \varepsilon}{\frac{t_0}{t_1}} + {\varepsilon}^{3} \cdot \left(t_3 - t_8\right)\right)\right) - \mathsf{fma}\left(-0.3333333333333333, \frac{t_1}{t_0}, \left(t_8 - t_3\right) \cdot \frac{\sin x}{\cos x}\right) \cdot {\varepsilon}^{4}\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t_2, \frac{-1}{\log \left(e^{t_5}\right)}, t_6\right)\\
\end{array}