double code(double x, double eps) {
return tan((x + eps)) - tan(x);
}
↓
double code(double x, double eps) {
double t_0 = tan(x) * tan(eps);
double t_1 = tan(x) + tan(eps);
double t_2 = pow(cos(x), 2.0);
double t_3 = pow(sin(x), 2.0);
double tmp;
if (eps <= -6.5e-5) {
tmp = (t_1 / (1.0 - (tan(x) / (1.0 / tan(eps))))) - tan(x);
} else if (eps <= 6.8e-19) {
tmp = (tan(eps) / (1.0 - t_0)) + ((eps / (t_2 / t_3)) + (((eps * eps) / (pow(cos(x), 3.0) / pow(sin(x), 3.0))) - (pow(eps, 3.0) * fma(-1.0, (pow(sin(x), 4.0) / pow(cos(x), 4.0)), (-0.3333333333333333 * (t_3 / t_2))))));
} else {
tmp = fma(t_1, (-1.0 / (t_0 + -1.0)), -tan(x));
}
return tmp;
}
function code(x, eps)
return Float64(tan(Float64(x + eps)) - tan(x))
end
↓
function code(x, eps)
t_0 = Float64(tan(x) * tan(eps))
t_1 = Float64(tan(x) + tan(eps))
t_2 = cos(x) ^ 2.0
t_3 = sin(x) ^ 2.0
tmp = 0.0
if (eps <= -6.5e-5)
tmp = Float64(Float64(t_1 / Float64(1.0 - Float64(tan(x) / Float64(1.0 / tan(eps))))) - tan(x));
elseif (eps <= 6.8e-19)
tmp = Float64(Float64(tan(eps) / Float64(1.0 - t_0)) + Float64(Float64(eps / Float64(t_2 / t_3)) + Float64(Float64(Float64(eps * eps) / Float64((cos(x) ^ 3.0) / (sin(x) ^ 3.0))) - Float64((eps ^ 3.0) * fma(-1.0, Float64((sin(x) ^ 4.0) / (cos(x) ^ 4.0)), Float64(-0.3333333333333333 * Float64(t_3 / t_2)))))));
else
tmp = fma(t_1, Float64(-1.0 / Float64(t_0 + -1.0)), 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[Tan[x], $MachinePrecision] * N[Tan[eps], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[Tan[x], $MachinePrecision] + N[Tan[eps], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[Power[N[Cos[x], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$3 = N[Power[N[Sin[x], $MachinePrecision], 2.0], $MachinePrecision]}, If[LessEqual[eps, -6.5e-5], N[(N[(t$95$1 / N[(1.0 - N[(N[Tan[x], $MachinePrecision] / N[(1.0 / N[Tan[eps], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[Tan[x], $MachinePrecision]), $MachinePrecision], If[LessEqual[eps, 6.8e-19], N[(N[(N[Tan[eps], $MachinePrecision] / N[(1.0 - t$95$0), $MachinePrecision]), $MachinePrecision] + N[(N[(eps / N[(t$95$2 / t$95$3), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(eps * eps), $MachinePrecision] / N[(N[Power[N[Cos[x], $MachinePrecision], 3.0], $MachinePrecision] / N[Power[N[Sin[x], $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[Power[eps, 3.0], $MachinePrecision] * N[(-1.0 * N[(N[Power[N[Sin[x], $MachinePrecision], 4.0], $MachinePrecision] / N[Power[N[Cos[x], $MachinePrecision], 4.0], $MachinePrecision]), $MachinePrecision] + N[(-0.3333333333333333 * N[(t$95$3 / t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 * N[(-1.0 / N[(t$95$0 + -1.0), $MachinePrecision]), $MachinePrecision] + (-N[Tan[x], $MachinePrecision])), $MachinePrecision]]]]]]]