\frac{\left(-b_2\right) - \sqrt{b_2 \cdot b_2 - a \cdot c}}{a}\begin{array}{l}
\mathbf{if}\;b_2 \le -4.2887136042886476 \cdot 10^{+71}:\\
\;\;\;\;\frac{-1}{2} \cdot \frac{c}{b_2}\\
\mathbf{elif}\;b_2 \le -3.407079315314288 \cdot 10^{-176}:\\
\;\;\;\;\frac{a \cdot \frac{c}{{\left(\sqrt{e}\right)}^{\left(\log \left(b_2 \cdot b_2 - c \cdot a\right)\right)} - b_2}}{a}\\
\mathbf{elif}\;b_2 \le 8.016779424032652 \cdot 10^{+82}:\\
\;\;\;\;\frac{1}{\frac{a}{\left(-b_2\right) - \sqrt{b_2 \cdot b_2 - c \cdot a}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{c}{b_2} \cdot \frac{1}{2} - 2 \cdot \frac{b_2}{a}\\
\end{array}double f(double a, double b_2, double c) {
double r1183852 = b_2;
double r1183853 = -r1183852;
double r1183854 = r1183852 * r1183852;
double r1183855 = a;
double r1183856 = c;
double r1183857 = r1183855 * r1183856;
double r1183858 = r1183854 - r1183857;
double r1183859 = sqrt(r1183858);
double r1183860 = r1183853 - r1183859;
double r1183861 = r1183860 / r1183855;
return r1183861;
}
double f(double a, double b_2, double c) {
double r1183862 = b_2;
double r1183863 = -4.2887136042886476e+71;
bool r1183864 = r1183862 <= r1183863;
double r1183865 = -0.5;
double r1183866 = c;
double r1183867 = r1183866 / r1183862;
double r1183868 = r1183865 * r1183867;
double r1183869 = -3.407079315314288e-176;
bool r1183870 = r1183862 <= r1183869;
double r1183871 = a;
double r1183872 = exp(1.0);
double r1183873 = sqrt(r1183872);
double r1183874 = r1183862 * r1183862;
double r1183875 = r1183866 * r1183871;
double r1183876 = r1183874 - r1183875;
double r1183877 = log(r1183876);
double r1183878 = pow(r1183873, r1183877);
double r1183879 = r1183878 - r1183862;
double r1183880 = r1183866 / r1183879;
double r1183881 = r1183871 * r1183880;
double r1183882 = r1183881 / r1183871;
double r1183883 = 8.016779424032652e+82;
bool r1183884 = r1183862 <= r1183883;
double r1183885 = 1.0;
double r1183886 = -r1183862;
double r1183887 = sqrt(r1183876);
double r1183888 = r1183886 - r1183887;
double r1183889 = r1183871 / r1183888;
double r1183890 = r1183885 / r1183889;
double r1183891 = 0.5;
double r1183892 = r1183867 * r1183891;
double r1183893 = 2.0;
double r1183894 = r1183862 / r1183871;
double r1183895 = r1183893 * r1183894;
double r1183896 = r1183892 - r1183895;
double r1183897 = r1183884 ? r1183890 : r1183896;
double r1183898 = r1183870 ? r1183882 : r1183897;
double r1183899 = r1183864 ? r1183868 : r1183898;
return r1183899;
}



Bits error versus a



Bits error versus b_2



Bits error versus c
Results
if b_2 < -4.2887136042886476e+71Initial program 57.3
Taylor expanded around -inf 3.3
if -4.2887136042886476e+71 < b_2 < -3.407079315314288e-176Initial program 36.0
rmApplied flip--36.1
Simplified15.5
Simplified15.5
rmApplied *-un-lft-identity15.5
Applied times-frac12.6
Simplified12.6
rmApplied add-exp-log15.3
rmApplied pow1/215.3
Applied log-pow15.3
Applied exp-prod15.7
Simplified15.7
if -3.407079315314288e-176 < b_2 < 8.016779424032652e+82Initial program 12.0
rmApplied clear-num12.1
if 8.016779424032652e+82 < b_2 Initial program 42.4
Taylor expanded around inf 3.6
Final simplification9.1
herbie shell --seed 2019163
(FPCore (a b_2 c)
:name "quad2m (problem 3.2.1, negative)"
(/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a))