double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
return ((x * ((y * z) - (t * a))) - (b * ((c * z) - (t * i)))) + (j * ((c * a) - (y * i)));
}
↓
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = c * (z * b);
double t_2 = c * (a * j);
double t_3 = i * fma(-b, t, (t * b));
double t_4 = x * ((y * z) - (t * a));
double t_5 = t_4 - (b * ((z * c) - (t * i)));
double t_6 = j * ((a * c) - (y * i));
double t_7 = t_5 + t_6;
double tmp;
if (t_7 <= -((double) INFINITY)) {
tmp = (t_2 - ((i * ((y * j) - (t * b))) + (t_3 + t_3))) - t_1;
} else if (t_7 <= 2e+302) {
tmp = t_5 + ((j * fma(-i, y, (y * i))) + t_6);
} else {
tmp = ((i * ((t * b) - (y * j))) + (t_2 + t_4)) - t_1;
}
return tmp;
}
function code(x, y, z, t, a, b, c, i, j)
return Float64(Float64(Float64(x * Float64(Float64(y * z) - Float64(t * a))) - Float64(b * Float64(Float64(c * z) - Float64(t * i)))) + Float64(j * Float64(Float64(c * a) - Float64(y * i))))
end
↓
function code(x, y, z, t, a, b, c, i, j)
t_1 = Float64(c * Float64(z * b))
t_2 = Float64(c * Float64(a * j))
t_3 = Float64(i * fma(Float64(-b), t, Float64(t * b)))
t_4 = Float64(x * Float64(Float64(y * z) - Float64(t * a)))
t_5 = Float64(t_4 - Float64(b * Float64(Float64(z * c) - Float64(t * i))))
t_6 = Float64(j * Float64(Float64(a * c) - Float64(y * i)))
t_7 = Float64(t_5 + t_6)
tmp = 0.0
if (t_7 <= Float64(-Inf))
tmp = Float64(Float64(t_2 - Float64(Float64(i * Float64(Float64(y * j) - Float64(t * b))) + Float64(t_3 + t_3))) - t_1);
elseif (t_7 <= 2e+302)
tmp = Float64(t_5 + Float64(Float64(j * fma(Float64(-i), y, Float64(y * i))) + t_6));
else
tmp = Float64(Float64(Float64(i * Float64(Float64(t * b) - Float64(y * j))) + Float64(t_2 + t_4)) - t_1);
end
return tmp
end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := N[(N[(N[(x * N[(N[(y * z), $MachinePrecision] - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(b * N[(N[(c * z), $MachinePrecision] - N[(t * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(j * N[(N[(c * a), $MachinePrecision] - N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
↓
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(c * N[(z * b), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(c * N[(a * j), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(i * N[((-b) * t + N[(t * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(x * N[(N[(y * z), $MachinePrecision] - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[(t$95$4 - N[(b * N[(N[(z * c), $MachinePrecision] - N[(t * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$6 = N[(j * N[(N[(a * c), $MachinePrecision] - N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$7 = N[(t$95$5 + t$95$6), $MachinePrecision]}, If[LessEqual[t$95$7, (-Infinity)], N[(N[(t$95$2 - N[(N[(i * N[(N[(y * j), $MachinePrecision] - N[(t * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$3 + t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision], If[LessEqual[t$95$7, 2e+302], N[(t$95$5 + N[(N[(j * N[((-i) * y + N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$6), $MachinePrecision]), $MachinePrecision], N[(N[(N[(i * N[(N[(t * b), $MachinePrecision] - N[(y * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$2 + t$95$4), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision]]]]]]]]]]