\[\frac{-\sqrt{\left(2 \cdot \left(\left({B}^{2} - \left(4 \cdot A\right) \cdot C\right) \cdot F\right)\right) \cdot \left(\left(A + C\right) - \sqrt{{\left(A - C\right)}^{2} + {B}^{2}}\right)}}{{B}^{2} - \left(4 \cdot A\right) \cdot C}
\]
↓
\[\begin{array}{l}
t_0 := \left(B \cdot B + C \cdot C\right) - {\left(-C\right)}^{2}\\
t_1 := \left(C - \left(-C\right)\right) \cdot F\\
\mathbf{if}\;A \leq -5.5 \cdot 10^{+54}:\\
\;\;\;\;\left(-\sqrt{-16 \cdot \left(C \cdot F\right)}\right) \cdot \frac{0.25}{C}\\
\mathbf{elif}\;A \leq 2.6 \cdot 10^{+48}:\\
\;\;\;\;\frac{-\sqrt{\left(\left(F + F\right) \cdot \left(C \cdot \left(-4 \cdot A\right)\right) + \left(F + F\right) \cdot \left(B \cdot B\right)\right) \cdot \left(C - \left(\mathsf{hypot}\left(A - C, B\right) - A\right)\right)}}{B \cdot B - 4 \cdot \left(A \cdot C\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{\mathsf{fma}\left(2, \mathsf{fma}\left(2, t_0 \cdot \left(C \cdot F\right), \left(B \cdot B\right) \cdot t_1\right), \mathsf{fma}\left(2, \frac{\mathsf{fma}\left(-0.5, \left(B \cdot B\right) \cdot \left(t_0 \cdot F\right), 2 \cdot \left(t_0 \cdot \left(\left(C \cdot C\right) \cdot F\right)\right)\right)}{A}, \left(-8 \cdot A\right) \cdot \left(t_1 \cdot C\right)\right)\right)}}{\mathsf{fma}\left(4, A \cdot C, -B \cdot B\right)}\\
\end{array}
\]
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 = ((B * B) + (C * C)) - pow(-C, 2.0);
double t_1 = (C - -C) * F;
double tmp;
if (A <= -5.5e+54) {
tmp = -sqrt((-16.0 * (C * F))) * (0.25 / C);
} else if (A <= 2.6e+48) {
tmp = -sqrt(((((F + F) * (C * (-4.0 * A))) + ((F + F) * (B * B))) * (C - (hypot((A - C), B) - A)))) / ((B * B) - (4.0 * (A * C)));
} else {
tmp = sqrt(fma(2.0, fma(2.0, (t_0 * (C * F)), ((B * B) * t_1)), fma(2.0, (fma(-0.5, ((B * B) * (t_0 * F)), (2.0 * (t_0 * ((C * C) * F)))) / A), ((-8.0 * A) * (t_1 * C))))) / fma(4.0, (A * C), -(B * B));
}
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(Float64(Float64(B * B) + Float64(C * C)) - (Float64(-C) ^ 2.0))
t_1 = Float64(Float64(C - Float64(-C)) * F)
tmp = 0.0
if (A <= -5.5e+54)
tmp = Float64(Float64(-sqrt(Float64(-16.0 * Float64(C * F)))) * Float64(0.25 / C));
elseif (A <= 2.6e+48)
tmp = Float64(Float64(-sqrt(Float64(Float64(Float64(Float64(F + F) * Float64(C * Float64(-4.0 * A))) + Float64(Float64(F + F) * Float64(B * B))) * Float64(C - Float64(hypot(Float64(A - C), B) - A))))) / Float64(Float64(B * B) - Float64(4.0 * Float64(A * C))));
else
tmp = Float64(sqrt(fma(2.0, fma(2.0, Float64(t_0 * Float64(C * F)), Float64(Float64(B * B) * t_1)), fma(2.0, Float64(fma(-0.5, Float64(Float64(B * B) * Float64(t_0 * F)), Float64(2.0 * Float64(t_0 * Float64(Float64(C * C) * F)))) / A), Float64(Float64(-8.0 * A) * Float64(t_1 * C))))) / fma(4.0, Float64(A * C), Float64(-Float64(B * B))));
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[(N[(N[(B * B), $MachinePrecision] + N[(C * C), $MachinePrecision]), $MachinePrecision] - N[Power[(-C), 2.0], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(C - (-C)), $MachinePrecision] * F), $MachinePrecision]}, If[LessEqual[A, -5.5e+54], N[((-N[Sqrt[N[(-16.0 * N[(C * F), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) * N[(0.25 / C), $MachinePrecision]), $MachinePrecision], If[LessEqual[A, 2.6e+48], N[((-N[Sqrt[N[(N[(N[(N[(F + F), $MachinePrecision] * N[(C * N[(-4.0 * A), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(F + F), $MachinePrecision] * N[(B * B), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(C - N[(N[Sqrt[N[(A - C), $MachinePrecision] ^ 2 + B ^ 2], $MachinePrecision] - A), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]) / N[(N[(B * B), $MachinePrecision] - N[(4.0 * N[(A * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Sqrt[N[(2.0 * N[(2.0 * N[(t$95$0 * N[(C * F), $MachinePrecision]), $MachinePrecision] + N[(N[(B * B), $MachinePrecision] * t$95$1), $MachinePrecision]), $MachinePrecision] + N[(2.0 * N[(N[(-0.5 * N[(N[(B * B), $MachinePrecision] * N[(t$95$0 * F), $MachinePrecision]), $MachinePrecision] + N[(2.0 * N[(t$95$0 * N[(N[(C * C), $MachinePrecision] * F), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / A), $MachinePrecision] + N[(N[(-8.0 * A), $MachinePrecision] * N[(t$95$1 * C), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / N[(4.0 * N[(A * C), $MachinePrecision] + (-N[(B * B), $MachinePrecision])), $MachinePrecision]), $MachinePrecision]]]]]
\frac{-\sqrt{\left(2 \cdot \left(\left({B}^{2} - \left(4 \cdot A\right) \cdot C\right) \cdot F\right)\right) \cdot \left(\left(A + C\right) - \sqrt{{\left(A - C\right)}^{2} + {B}^{2}}\right)}}{{B}^{2} - \left(4 \cdot A\right) \cdot C}
↓
\begin{array}{l}
t_0 := \left(B \cdot B + C \cdot C\right) - {\left(-C\right)}^{2}\\
t_1 := \left(C - \left(-C\right)\right) \cdot F\\
\mathbf{if}\;A \leq -5.5 \cdot 10^{+54}:\\
\;\;\;\;\left(-\sqrt{-16 \cdot \left(C \cdot F\right)}\right) \cdot \frac{0.25}{C}\\
\mathbf{elif}\;A \leq 2.6 \cdot 10^{+48}:\\
\;\;\;\;\frac{-\sqrt{\left(\left(F + F\right) \cdot \left(C \cdot \left(-4 \cdot A\right)\right) + \left(F + F\right) \cdot \left(B \cdot B\right)\right) \cdot \left(C - \left(\mathsf{hypot}\left(A - C, B\right) - A\right)\right)}}{B \cdot B - 4 \cdot \left(A \cdot C\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{\mathsf{fma}\left(2, \mathsf{fma}\left(2, t_0 \cdot \left(C \cdot F\right), \left(B \cdot B\right) \cdot t_1\right), \mathsf{fma}\left(2, \frac{\mathsf{fma}\left(-0.5, \left(B \cdot B\right) \cdot \left(t_0 \cdot F\right), 2 \cdot \left(t_0 \cdot \left(\left(C \cdot C\right) \cdot F\right)\right)\right)}{A}, \left(-8 \cdot A\right) \cdot \left(t_1 \cdot C\right)\right)\right)}}{\mathsf{fma}\left(4, A \cdot C, -B \cdot B\right)}\\
\end{array}