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