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 = (t * c) - (y * i);
double t_2 = (j * t_1) - ((b * ((z * c) - (a * i))) + (x * ((t * a) - (y * z))));
double tmp;
if ((t_2 <= -((double) INFINITY)) || !(t_2 <= 1e+304)) {
tmp = (t * ((c * j) - (x * a))) + ((i * (a * b)) - (y * ((i * j) - (x * z))));
} else {
tmp = fma(j, t_1, (2.0 * (j * fma(-y, i, (y * i))))) + ((x * ((y * z) - (t * a))) + (b * ((a * i) - (z * c))));
}
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(Float64(t * c) - Float64(y * i))
t_2 = Float64(Float64(j * t_1) - Float64(Float64(b * Float64(Float64(z * c) - Float64(a * i))) + Float64(x * Float64(Float64(t * a) - Float64(y * z)))))
tmp = 0.0
if ((t_2 <= Float64(-Inf)) || !(t_2 <= 1e+304))
tmp = Float64(Float64(t * Float64(Float64(c * j) - Float64(x * a))) + Float64(Float64(i * Float64(a * b)) - Float64(y * Float64(Float64(i * j) - Float64(x * z)))));
else
tmp = Float64(fma(j, t_1, Float64(2.0 * Float64(j * fma(Float64(-y), i, Float64(y * i))))) + Float64(Float64(x * Float64(Float64(y * z) - Float64(t * a))) + Float64(b * Float64(Float64(a * i) - Float64(z * c)))));
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[(N[(t * c), $MachinePrecision] - N[(y * i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(j * t$95$1), $MachinePrecision] - N[(N[(b * N[(N[(z * c), $MachinePrecision] - N[(a * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x * N[(N[(t * a), $MachinePrecision] - N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[t$95$2, (-Infinity)], N[Not[LessEqual[t$95$2, 1e+304]], $MachinePrecision]], N[(N[(t * N[(N[(c * j), $MachinePrecision] - N[(x * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(i * N[(a * b), $MachinePrecision]), $MachinePrecision] - N[(y * N[(N[(i * j), $MachinePrecision] - N[(x * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(j * t$95$1 + N[(2.0 * N[(j * N[((-y) * i + N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(x * N[(N[(y * z), $MachinePrecision] - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(b * N[(N[(a * i), $MachinePrecision] - N[(z * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]