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) - (i * a)))) + (j * ((c * t) - (i * y)));
}
↓
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double t_1 = b * ((a * i) - (z * c));
double t_2 = x * ((y * z) - (t * a));
double t_3 = j * ((t * c) - (y * i));
double t_4 = (t_2 + t_1) + t_3;
double tmp;
if ((t_4 <= -((double) INFINITY)) || !(t_4 <= ((double) INFINITY))) {
tmp = ((y * ((x * z) - (i * j))) - (t * ((x * a) - (c * j)))) + (i * (a * b));
} else {
tmp = (t_2 + (t_1 - (b * fma(-a, i, (a * i))))) + t_3;
}
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(i * a)))) + Float64(j * Float64(Float64(c * t) - Float64(i * y))))
end
↓
function code(x, y, z, t, a, b, c, i, j)
t_1 = Float64(b * Float64(Float64(a * i) - Float64(z * c)))
t_2 = Float64(x * Float64(Float64(y * z) - Float64(t * a)))
t_3 = Float64(j * Float64(Float64(t * c) - Float64(y * i)))
t_4 = Float64(Float64(t_2 + t_1) + t_3)
tmp = 0.0
if ((t_4 <= Float64(-Inf)) || !(t_4 <= Inf))
tmp = Float64(Float64(Float64(y * Float64(Float64(x * z) - Float64(i * j))) - Float64(t * Float64(Float64(x * a) - Float64(c * j)))) + Float64(i * Float64(a * b)));
else
tmp = Float64(Float64(t_2 + Float64(t_1 - Float64(b * fma(Float64(-a), i, Float64(a * i))))) + t_3);
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[(i * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(j * N[(N[(c * t), $MachinePrecision] - N[(i * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
↓
code[x_, y_, z_, t_, a_, b_, c_, i_, j_] := Block[{t$95$1 = N[(b * N[(N[(a * i), $MachinePrecision] - N[(z * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x * N[(N[(y * z), $MachinePrecision] - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(j * N[(N[(t * c), $MachinePrecision] - N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(N[(t$95$2 + t$95$1), $MachinePrecision] + t$95$3), $MachinePrecision]}, If[Or[LessEqual[t$95$4, (-Infinity)], N[Not[LessEqual[t$95$4, Infinity]], $MachinePrecision]], N[(N[(N[(y * N[(N[(x * z), $MachinePrecision] - N[(i * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(t * N[(N[(x * a), $MachinePrecision] - N[(c * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(i * N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(t$95$2 + N[(t$95$1 - N[(b * N[((-a) * i + N[(a * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$3), $MachinePrecision]]]]]]