\frac{\left(-b_2\right) - \sqrt{b_2 \cdot b_2 - a \cdot c}}{a}\begin{array}{l}
\mathbf{if}\;b_2 \le -1.8746290509448952 \cdot 10^{74}:\\
\;\;\;\;1 \cdot \left(\frac{-1}{2} \cdot \frac{c}{b_2}\right)\\
\mathbf{elif}\;b_2 \le -6.57941275339975712 \cdot 10^{-183}:\\
\;\;\;\;\frac{\frac{1 \cdot \left(\left({b_2}^{2} - {b_2}^{2}\right) + a \cdot c\right)}{a}}{\left(-b_2\right) + \sqrt{b_2 \cdot b_2 - a \cdot c}}\\
\mathbf{elif}\;b_2 \le 2.89307072747620556 \cdot 10^{122}:\\
\;\;\;\;1 \cdot \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 r14499 = b_2;
double r14500 = -r14499;
double r14501 = r14499 * r14499;
double r14502 = a;
double r14503 = c;
double r14504 = r14502 * r14503;
double r14505 = r14501 - r14504;
double r14506 = sqrt(r14505);
double r14507 = r14500 - r14506;
double r14508 = r14507 / r14502;
return r14508;
}
double f(double a, double b_2, double c) {
double r14509 = b_2;
double r14510 = -1.8746290509448952e+74;
bool r14511 = r14509 <= r14510;
double r14512 = 1.0;
double r14513 = -0.5;
double r14514 = c;
double r14515 = r14514 / r14509;
double r14516 = r14513 * r14515;
double r14517 = r14512 * r14516;
double r14518 = -6.579412753399757e-183;
bool r14519 = r14509 <= r14518;
double r14520 = 2.0;
double r14521 = pow(r14509, r14520);
double r14522 = r14521 - r14521;
double r14523 = a;
double r14524 = r14523 * r14514;
double r14525 = r14522 + r14524;
double r14526 = r14512 * r14525;
double r14527 = r14526 / r14523;
double r14528 = -r14509;
double r14529 = r14509 * r14509;
double r14530 = r14529 - r14524;
double r14531 = sqrt(r14530);
double r14532 = r14528 + r14531;
double r14533 = r14527 / r14532;
double r14534 = 2.8930707274762056e+122;
bool r14535 = r14509 <= r14534;
double r14536 = r14528 - r14531;
double r14537 = r14536 / r14523;
double r14538 = r14512 * r14537;
double r14539 = -2.0;
double r14540 = r14509 / r14523;
double r14541 = r14539 * r14540;
double r14542 = r14535 ? r14538 : r14541;
double r14543 = r14519 ? r14533 : r14542;
double r14544 = r14511 ? r14517 : r14543;
return r14544;
}



Bits error versus a



Bits error versus b_2



Bits error versus c
Results
if b_2 < -1.8746290509448952e+74Initial program 58.6
rmApplied clear-num58.6
rmApplied *-un-lft-identity58.6
Applied *-un-lft-identity58.6
Applied times-frac58.6
Applied add-cube-cbrt58.6
Applied times-frac58.6
Simplified58.6
Simplified58.6
Taylor expanded around -inf 3.3
if -1.8746290509448952e+74 < b_2 < -6.579412753399757e-183Initial program 36.6
rmApplied clear-num36.7
rmApplied flip--36.7
Applied associate-/r/36.7
Applied associate-/r*36.7
Simplified15.8
if -6.579412753399757e-183 < b_2 < 2.8930707274762056e+122Initial program 10.2
rmApplied clear-num10.4
rmApplied *-un-lft-identity10.4
Applied *-un-lft-identity10.4
Applied times-frac10.4
Applied add-cube-cbrt10.4
Applied times-frac10.4
Simplified10.4
Simplified10.2
if 2.8930707274762056e+122 < b_2 Initial program 52.8
rmApplied clear-num52.9
Taylor expanded around 0 3.3
Final simplification8.6
herbie shell --seed 2020036
(FPCore (a b_2 c)
:name "quad2m (problem 3.2.1, negative)"
:precision binary64
(/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a))