\frac{\left(-b_2\right) - \sqrt{b_2 \cdot b_2 - a \cdot c}}{a}\begin{array}{l}
\mathbf{if}\;b_2 \leq -8.055509887818825 \cdot 10^{+155}:\\
\;\;\;\;\frac{c}{\left(a \cdot \frac{0.5}{\frac{b_2}{c}} - b_2\right) - b_2}\\
\mathbf{elif}\;b_2 \leq 1.2144381520308056 \cdot 10^{-219}:\\
\;\;\;\;\frac{c}{\sqrt{b_2 \cdot b_2 - c \cdot a} - b_2}\\
\mathbf{elif}\;b_2 \leq 2.706784727683868 \cdot 10^{+149}:\\
\;\;\;\;\left(b_2 + \sqrt{b_2 \cdot b_2 - c \cdot a}\right) \cdot \frac{-1}{a}\\
\mathbf{else}:\\
\;\;\;\;\frac{b_2 \cdot -2}{a}\\
\end{array}double code(double a, double b_2, double c) {
return (((double) (((double) -(b_2)) - ((double) sqrt(((double) (((double) (b_2 * b_2)) - ((double) (a * c)))))))) / a);
}
double code(double a, double b_2, double c) {
double VAR;
if ((b_2 <= -8.055509887818825e+155)) {
VAR = (c / ((double) (((double) (((double) (a * (0.5 / (b_2 / c)))) - b_2)) - b_2)));
} else {
double VAR_1;
if ((b_2 <= 1.2144381520308056e-219)) {
VAR_1 = (c / ((double) (((double) sqrt(((double) (((double) (b_2 * b_2)) - ((double) (c * a)))))) - b_2)));
} else {
double VAR_2;
if ((b_2 <= 2.706784727683868e+149)) {
VAR_2 = ((double) (((double) (b_2 + ((double) sqrt(((double) (((double) (b_2 * b_2)) - ((double) (c * a)))))))) * (-1.0 / a)));
} else {
VAR_2 = (((double) (b_2 * -2.0)) / a);
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus a



Bits error versus b_2



Bits error versus c
Results
if b_2 < -8.0555098878188251e155Initial program Error: 64.0 bits
rmApplied flip--Error: 64.0 bits
SimplifiedError: 38.4 bits
SimplifiedError: 38.4 bits
rmApplied *-un-lft-identityError: 38.4 bits
Applied *-un-lft-identityError: 38.4 bits
Applied times-fracError: 38.4 bits
SimplifiedError: 38.4 bits
SimplifiedError: 38.2 bits
Taylor expanded around -inf Error: 6.9 bits
SimplifiedError: 1.2 bits
if -8.0555098878188251e155 < b_2 < 1.21443815203080564e-219Initial program Error: 32.7 bits
rmApplied flip--Error: 32.7 bits
SimplifiedError: 15.8 bits
SimplifiedError: 15.8 bits
rmApplied *-un-lft-identityError: 15.8 bits
Applied *-un-lft-identityError: 15.8 bits
Applied times-fracError: 15.8 bits
SimplifiedError: 15.8 bits
SimplifiedError: 9.0 bits
if 1.21443815203080564e-219 < b_2 < 2.70678472768386803e149Initial program Error: 7.2 bits
rmApplied div-invError: 7.4 bits
if 2.70678472768386803e149 < b_2 Initial program Error: 62.0 bits
rmApplied flip--Error: 64.0 bits
SimplifiedError: 62.8 bits
SimplifiedError: 62.8 bits
Taylor expanded around 0 Error: 1.9 bits
SimplifiedError: 1.9 bits
Final simplificationError: 6.4 bits
herbie shell --seed 2020200
(FPCore (a b_2 c)
:name "quad2m (problem 3.2.1, negative)"
:precision binary64
(/ (- (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a))