double code(double A, double B, double C, double F) {
return -sqrt(((2.0 * ((pow(B, 2.0) - ((4.0 * A) * C)) * F)) * ((A + C) + sqrt((pow((A - C), 2.0) + pow(B, 2.0)))))) / (pow(B, 2.0) - ((4.0 * A) * C));
}
↓
double code(double A, double B, double C, double F) {
double t_0 = C * (A * -4.0);
double t_1 = -4.0 * (A * C);
double t_2 = hypot(B, (A - C));
double t_3 = fma(B, B, t_0);
double t_4 = pow(B, 2.0) + t_0;
double t_5 = -sqrt(((2.0 * (t_4 * F)) * ((A + C) + sqrt((pow(B, 2.0) + pow((A - C), 2.0)))))) / t_4;
double tmp;
if (t_5 <= 0.0) {
tmp = ((sqrt(((A + C) + t_2)) * sqrt(F)) * -sqrt((2.0 * t_3))) / t_3;
} else if (t_5 <= ((double) INFINITY)) {
tmp = (sqrt((2.0 * (F * fma(B, B, t_1)))) * -sqrt((A + (C + t_2)))) / (t_1 + (B * B));
} else {
tmp = (sqrt(2.0) / B) * (sqrt((C + fma(0.5, (C * (C / B)), B))) * -sqrt(F));
}
return tmp;
}
function code(A, B, C, F)
return Float64(Float64(-sqrt(Float64(Float64(2.0 * Float64(Float64((B ^ 2.0) - Float64(Float64(4.0 * A) * C)) * F)) * Float64(Float64(A + C) + sqrt(Float64((Float64(A - C) ^ 2.0) + (B ^ 2.0))))))) / Float64((B ^ 2.0) - Float64(Float64(4.0 * A) * C)))
end
↓
function code(A, B, C, F)
t_0 = Float64(C * Float64(A * -4.0))
t_1 = Float64(-4.0 * Float64(A * C))
t_2 = hypot(B, Float64(A - C))
t_3 = fma(B, B, t_0)
t_4 = Float64((B ^ 2.0) + t_0)
t_5 = Float64(Float64(-sqrt(Float64(Float64(2.0 * Float64(t_4 * F)) * Float64(Float64(A + C) + sqrt(Float64((B ^ 2.0) + (Float64(A - C) ^ 2.0))))))) / t_4)
tmp = 0.0
if (t_5 <= 0.0)
tmp = Float64(Float64(Float64(sqrt(Float64(Float64(A + C) + t_2)) * sqrt(F)) * Float64(-sqrt(Float64(2.0 * t_3)))) / t_3);
elseif (t_5 <= Inf)
tmp = Float64(Float64(sqrt(Float64(2.0 * Float64(F * fma(B, B, t_1)))) * Float64(-sqrt(Float64(A + Float64(C + t_2))))) / Float64(t_1 + Float64(B * B)));
else
tmp = Float64(Float64(sqrt(2.0) / B) * Float64(sqrt(Float64(C + fma(0.5, Float64(C * Float64(C / B)), B))) * Float64(-sqrt(F))));
end
return tmp
end
code[A_, B_, C_, F_] := N[((-N[Sqrt[N[(N[(2.0 * N[(N[(N[Power[B, 2.0], $MachinePrecision] - N[(N[(4.0 * A), $MachinePrecision] * C), $MachinePrecision]), $MachinePrecision] * F), $MachinePrecision]), $MachinePrecision] * N[(N[(A + C), $MachinePrecision] + N[Sqrt[N[(N[Power[N[(A - C), $MachinePrecision], 2.0], $MachinePrecision] + N[Power[B, 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / N[(N[Power[B, 2.0], $MachinePrecision] - N[(N[(4.0 * A), $MachinePrecision] * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
↓
code[A_, B_, C_, F_] := Block[{t$95$0 = N[(C * N[(A * -4.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(-4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[Sqrt[B ^ 2 + N[(A - C), $MachinePrecision] ^ 2], $MachinePrecision]}, Block[{t$95$3 = N[(B * B + t$95$0), $MachinePrecision]}, Block[{t$95$4 = N[(N[Power[B, 2.0], $MachinePrecision] + t$95$0), $MachinePrecision]}, Block[{t$95$5 = N[((-N[Sqrt[N[(N[(2.0 * N[(t$95$4 * F), $MachinePrecision]), $MachinePrecision] * N[(N[(A + C), $MachinePrecision] + N[Sqrt[N[(N[Power[B, 2.0], $MachinePrecision] + N[Power[N[(A - C), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / t$95$4), $MachinePrecision]}, If[LessEqual[t$95$5, 0.0], N[(N[(N[(N[Sqrt[N[(N[(A + C), $MachinePrecision] + t$95$2), $MachinePrecision]], $MachinePrecision] * N[Sqrt[F], $MachinePrecision]), $MachinePrecision] * (-N[Sqrt[N[(2.0 * t$95$3), $MachinePrecision]], $MachinePrecision])), $MachinePrecision] / t$95$3), $MachinePrecision], If[LessEqual[t$95$5, Infinity], N[(N[(N[Sqrt[N[(2.0 * N[(F * N[(B * B + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * (-N[Sqrt[N[(A + N[(C + t$95$2), $MachinePrecision]), $MachinePrecision]], $MachinePrecision])), $MachinePrecision] / N[(t$95$1 + N[(B * B), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[Sqrt[2.0], $MachinePrecision] / B), $MachinePrecision] * N[(N[Sqrt[N[(C + N[(0.5 * N[(C * N[(C / B), $MachinePrecision]), $MachinePrecision] + B), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * (-N[Sqrt[F], $MachinePrecision])), $MachinePrecision]), $MachinePrecision]]]]]]]]]