\frac{\left(-b_2\right) - \sqrt{b_2 \cdot b_2 - a \cdot c}}{a}\begin{array}{l}
\mathbf{if}\;b_2 \le -3.101473652193339 \cdot 10^{+126}:\\
\;\;\;\;\frac{-1}{2} \cdot \frac{c}{b_2}\\
\mathbf{elif}\;b_2 \le 3.3026083301884445 \cdot 10^{-291}:\\
\;\;\;\;\frac{\frac{c \cdot a}{a}}{\sqrt{b_2 \cdot b_2 - c \cdot a} - b_2}\\
\mathbf{elif}\;b_2 \le 1.1638796624534952 \cdot 10^{+125}:\\
\;\;\;\;\frac{\left(-b_2\right) - \sqrt{b_2 \cdot b_2 - c \cdot a}}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{2} \cdot \frac{c}{b_2} - \frac{b_2}{a} \cdot 2\\
\end{array}double f(double a, double b_2, double c) {
double r2035237 = b_2;
double r2035238 = -r2035237;
double r2035239 = r2035237 * r2035237;
double r2035240 = a;
double r2035241 = c;
double r2035242 = r2035240 * r2035241;
double r2035243 = r2035239 - r2035242;
double r2035244 = sqrt(r2035243);
double r2035245 = r2035238 - r2035244;
double r2035246 = r2035245 / r2035240;
return r2035246;
}
double f(double a, double b_2, double c) {
double r2035247 = b_2;
double r2035248 = -3.101473652193339e+126;
bool r2035249 = r2035247 <= r2035248;
double r2035250 = -0.5;
double r2035251 = c;
double r2035252 = r2035251 / r2035247;
double r2035253 = r2035250 * r2035252;
double r2035254 = 3.3026083301884445e-291;
bool r2035255 = r2035247 <= r2035254;
double r2035256 = a;
double r2035257 = r2035251 * r2035256;
double r2035258 = r2035257 / r2035256;
double r2035259 = r2035247 * r2035247;
double r2035260 = r2035259 - r2035257;
double r2035261 = sqrt(r2035260);
double r2035262 = r2035261 - r2035247;
double r2035263 = r2035258 / r2035262;
double r2035264 = 1.1638796624534952e+125;
bool r2035265 = r2035247 <= r2035264;
double r2035266 = -r2035247;
double r2035267 = r2035266 - r2035261;
double r2035268 = r2035267 / r2035256;
double r2035269 = 0.5;
double r2035270 = r2035269 * r2035252;
double r2035271 = r2035247 / r2035256;
double r2035272 = 2.0;
double r2035273 = r2035271 * r2035272;
double r2035274 = r2035270 - r2035273;
double r2035275 = r2035265 ? r2035268 : r2035274;
double r2035276 = r2035255 ? r2035263 : r2035275;
double r2035277 = r2035249 ? r2035253 : r2035276;
return r2035277;
}



Bits error versus a



Bits error versus b_2



Bits error versus c
Results
if b_2 < -3.101473652193339e+126Initial program 60.4
rmApplied flip--60.4
Simplified56.0
Simplified56.0
Taylor expanded around -inf 2.4
if -3.101473652193339e+126 < b_2 < 3.3026083301884445e-291Initial program 32.7
rmApplied flip--32.8
Simplified15.2
Simplified15.2
rmApplied *-un-lft-identity15.2
Applied *-un-lft-identity15.2
Applied *-un-lft-identity15.2
Applied times-frac15.2
Applied times-frac15.2
Simplified15.2
Simplified13.9
if 3.3026083301884445e-291 < b_2 < 1.1638796624534952e+125Initial program 8.6
if 1.1638796624534952e+125 < b_2 Initial program 50.6
rmApplied flip--61.9
Simplified61.9
Simplified61.9
Taylor expanded around inf 3.1
Final simplification8.5
herbie shell --seed 2019129
(FPCore (a b_2 c)
:name "NMSE problem 3.2.1"
(/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a))