\frac{\left(-b_2\right) - \sqrt{b_2 \cdot b_2 - a \cdot c}}{a}\begin{array}{l}
\mathbf{if}\;b_2 \le -5.674469085146396739103610609439188639717 \cdot 10^{110}:\\
\;\;\;\;\frac{\frac{c}{b_2}}{-2}\\
\mathbf{elif}\;b_2 \le \frac{2563685591414627}{3.630412374213337555659332270873780664383 \cdot 10^{280}}:\\
\;\;\;\;1 \cdot \frac{c}{\sqrt{b_2 \cdot b_2 - a \cdot c} - b_2}\\
\mathbf{elif}\;b_2 \le 1.715181108188238274259588142060201574853 \cdot 10^{78}:\\
\;\;\;\;\frac{\left(-b_2\right) - \sqrt{b_2 \cdot b_2 - a \cdot c}}{a}\\
\mathbf{else}:\\
\;\;\;\;-2 \cdot \frac{b_2}{a}\\
\end{array}double f(double a, double b_2, double c) {
double r30883 = b_2;
double r30884 = -r30883;
double r30885 = r30883 * r30883;
double r30886 = a;
double r30887 = c;
double r30888 = r30886 * r30887;
double r30889 = r30885 - r30888;
double r30890 = sqrt(r30889);
double r30891 = r30884 - r30890;
double r30892 = r30891 / r30886;
return r30892;
}
double f(double a, double b_2, double c) {
double r30893 = b_2;
double r30894 = -5.674469085146397e+110;
bool r30895 = r30893 <= r30894;
double r30896 = c;
double r30897 = r30896 / r30893;
double r30898 = -2.0;
double r30899 = r30897 / r30898;
double r30900 = 2563685591414627.0;
double r30901 = 3.6304123742133376e+280;
double r30902 = r30900 / r30901;
bool r30903 = r30893 <= r30902;
double r30904 = 1.0;
double r30905 = r30893 * r30893;
double r30906 = a;
double r30907 = r30906 * r30896;
double r30908 = r30905 - r30907;
double r30909 = sqrt(r30908);
double r30910 = r30909 - r30893;
double r30911 = r30896 / r30910;
double r30912 = r30904 * r30911;
double r30913 = 1.7151811081882383e+78;
bool r30914 = r30893 <= r30913;
double r30915 = -r30893;
double r30916 = r30915 - r30909;
double r30917 = r30916 / r30906;
double r30918 = r30893 / r30906;
double r30919 = r30898 * r30918;
double r30920 = r30914 ? r30917 : r30919;
double r30921 = r30903 ? r30912 : r30920;
double r30922 = r30895 ? r30899 : r30921;
return r30922;
}



Bits error versus a



Bits error versus b_2



Bits error versus c
Results
if b_2 < -5.674469085146397e+110Initial program 59.7
Taylor expanded around -inf 2.7
Simplified2.7
if -5.674469085146397e+110 < b_2 < 7.061692521831336e-266Initial program 31.8
rmApplied flip--31.8
Simplified16.1
Simplified16.1
rmApplied *-un-lft-identity16.1
Applied associate-/r*16.1
Simplified14.1
rmApplied *-un-lft-identity14.1
Applied *-un-lft-identity14.1
Applied times-frac14.1
Applied associate-/l*14.0
Simplified9.0
rmApplied div-inv9.0
Simplified8.7
if 7.061692521831336e-266 < b_2 < 1.7151811081882383e+78Initial program 8.5
if 1.7151811081882383e+78 < b_2 Initial program 43.0
rmApplied flip--62.6
Simplified61.8
Simplified61.8
rmApplied *-un-lft-identity61.8
Applied associate-/r*61.8
Simplified61.6
Taylor expanded around 0 4.8
Final simplification6.6
herbie shell --seed 2019303
(FPCore (a b_2 c)
:name "quad2m (problem 3.2.1, negative)"
:precision binary64
(/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a))