double code(double t, double l, double k) {
return 2.0 / ((((pow(t, 3.0) / (l * l)) * sin(k)) * tan(k)) * ((1.0 + pow((k / t), 2.0)) + 1.0));
}
↓
double code(double t, double l, double k) {
double t_1 = pow(sin(k), 2.0);
double t_2 = 2.0 / fma(2.0, ((t / cos(k)) * (t_1 / ((l / t) * (l / t)))), (t / ((cos(k) * (l * l)) / (t_1 * pow((-1.0 / k), -2.0)))));
double t_3 = l / (k * t);
double tmp;
if (k <= -5.9e+144) {
tmp = (l / k) * ((2.0 / t) * ((l / k) / (sin(k) * tan(k))));
} else if (k <= -3.7e-56) {
tmp = t_2;
} else if (k <= 2.6e-78) {
tmp = t_3 * (t_3 / t);
} else if (k <= 2.5e+121) {
tmp = t_2;
} else {
tmp = 2.0 / (((t * sin(k)) * (tan(k) * (k / l))) / (l / k));
}
return tmp;
}
function code(t, l, k)
return Float64(2.0 / Float64(Float64(Float64(Float64((t ^ 3.0) / Float64(l * l)) * sin(k)) * tan(k)) * Float64(Float64(1.0 + (Float64(k / t) ^ 2.0)) + 1.0)))
end
↓
function code(t, l, k)
t_1 = sin(k) ^ 2.0
t_2 = Float64(2.0 / fma(2.0, Float64(Float64(t / cos(k)) * Float64(t_1 / Float64(Float64(l / t) * Float64(l / t)))), Float64(t / Float64(Float64(cos(k) * Float64(l * l)) / Float64(t_1 * (Float64(-1.0 / k) ^ -2.0))))))
t_3 = Float64(l / Float64(k * t))
tmp = 0.0
if (k <= -5.9e+144)
tmp = Float64(Float64(l / k) * Float64(Float64(2.0 / t) * Float64(Float64(l / k) / Float64(sin(k) * tan(k)))));
elseif (k <= -3.7e-56)
tmp = t_2;
elseif (k <= 2.6e-78)
tmp = Float64(t_3 * Float64(t_3 / t));
elseif (k <= 2.5e+121)
tmp = t_2;
else
tmp = Float64(2.0 / Float64(Float64(Float64(t * sin(k)) * Float64(tan(k) * Float64(k / l))) / Float64(l / k)));
end
return tmp
end
code[t_, l_, k_] := N[(2.0 / N[(N[(N[(N[(N[Power[t, 3.0], $MachinePrecision] / N[(l * l), $MachinePrecision]), $MachinePrecision] * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision] * N[(N[(1.0 + N[Power[N[(k / t), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
↓
code[t_, l_, k_] := Block[{t$95$1 = N[Power[N[Sin[k], $MachinePrecision], 2.0], $MachinePrecision]}, Block[{t$95$2 = N[(2.0 / N[(2.0 * N[(N[(t / N[Cos[k], $MachinePrecision]), $MachinePrecision] * N[(t$95$1 / N[(N[(l / t), $MachinePrecision] * N[(l / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t / N[(N[(N[Cos[k], $MachinePrecision] * N[(l * l), $MachinePrecision]), $MachinePrecision] / N[(t$95$1 * N[Power[N[(-1.0 / k), $MachinePrecision], -2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(l / N[(k * t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[k, -5.9e+144], N[(N[(l / k), $MachinePrecision] * N[(N[(2.0 / t), $MachinePrecision] * N[(N[(l / k), $MachinePrecision] / N[(N[Sin[k], $MachinePrecision] * N[Tan[k], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, -3.7e-56], t$95$2, If[LessEqual[k, 2.6e-78], N[(t$95$3 * N[(t$95$3 / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[k, 2.5e+121], t$95$2, N[(2.0 / N[(N[(N[(t * N[Sin[k], $MachinePrecision]), $MachinePrecision] * N[(N[Tan[k], $MachinePrecision] * N[(k / l), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(l / k), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]