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 = j * ((c * a) - (y * i));
double t_2 = (t_1 + (j * fma(-i, y, (y * i)))) - ((c * (z * b)) + (x * ((t * a) - (z * y))));
double t_3 = z * ((y * x) - (c * b));
double t_4 = c * (a * j);
double t_5 = ((y * (z * x)) + (b * ((i * t) - (z * c)))) + (t_4 - (i * (y * j)));
double tmp;
if (z <= -1e+145) {
tmp = ((i * (b * t)) - (a * (x * t))) + t_3;
} else if (z <= -1.462298553028206e+36) {
tmp = t_5;
} else if (z <= -1.3360345713096714e+24) {
tmp = a * ((c * j) - (x * t));
} else if (z <= -5.706431800087509e-60) {
tmp = t_4 + ((x * ((z * y) - (t * a))) - (b * ((z * c) - (i * t))));
} else if (z <= -7.269775795587896e-114) {
tmp = t_2;
} else if (z <= 3.3206315369560735e-219) {
tmp = t_5;
} else if (z <= 1e+110) {
tmp = t_2;
} else if (z <= 1e+154) {
tmp = t_5;
} else {
tmp = t_1 + 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(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(j * Float64(Float64(c * a) - Float64(y * i)))
t_2 = Float64(Float64(t_1 + Float64(j * fma(Float64(-i), y, Float64(y * i)))) - Float64(Float64(c * Float64(z * b)) + Float64(x * Float64(Float64(t * a) - Float64(z * y)))))
t_3 = Float64(z * Float64(Float64(y * x) - Float64(c * b)))
t_4 = Float64(c * Float64(a * j))
t_5 = Float64(Float64(Float64(y * Float64(z * x)) + Float64(b * Float64(Float64(i * t) - Float64(z * c)))) + Float64(t_4 - Float64(i * Float64(y * j))))
tmp = 0.0
if (z <= -1e+145)
tmp = Float64(Float64(Float64(i * Float64(b * t)) - Float64(a * Float64(x * t))) + t_3);
elseif (z <= -1.462298553028206e+36)
tmp = t_5;
elseif (z <= -1.3360345713096714e+24)
tmp = Float64(a * Float64(Float64(c * j) - Float64(x * t)));
elseif (z <= -5.706431800087509e-60)
tmp = Float64(t_4 + Float64(Float64(x * Float64(Float64(z * y) - Float64(t * a))) - Float64(b * Float64(Float64(z * c) - Float64(i * t)))));
elseif (z <= -7.269775795587896e-114)
tmp = t_2;
elseif (z <= 3.3206315369560735e-219)
tmp = t_5;
elseif (z <= 1e+110)
tmp = t_2;
elseif (z <= 1e+154)
tmp = t_5;
else
tmp = Float64(t_1 + 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[(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[(j * N[(N[(c * a), $MachinePrecision] - N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(t$95$1 + N[(j * N[((-i) * y + N[(y * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(c * N[(z * b), $MachinePrecision]), $MachinePrecision] + N[(x * N[(N[(t * a), $MachinePrecision] - N[(z * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(z * N[(N[(y * x), $MachinePrecision] - N[(c * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(c * N[(a * j), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[(N[(N[(y * N[(z * x), $MachinePrecision]), $MachinePrecision] + N[(b * N[(N[(i * t), $MachinePrecision] - N[(z * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$4 - N[(i * N[(y * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1e+145], N[(N[(N[(i * N[(b * t), $MachinePrecision]), $MachinePrecision] - N[(a * N[(x * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t$95$3), $MachinePrecision], If[LessEqual[z, -1.462298553028206e+36], t$95$5, If[LessEqual[z, -1.3360345713096714e+24], N[(a * N[(N[(c * j), $MachinePrecision] - N[(x * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -5.706431800087509e-60], N[(t$95$4 + N[(N[(x * N[(N[(z * y), $MachinePrecision] - N[(t * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(b * N[(N[(z * c), $MachinePrecision] - N[(i * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -7.269775795587896e-114], t$95$2, If[LessEqual[z, 3.3206315369560735e-219], t$95$5, If[LessEqual[z, 1e+110], t$95$2, If[LessEqual[z, 1e+154], t$95$5, N[(t$95$1 + t$95$3), $MachinePrecision]]]]]]]]]]]]]]