\frac{\left(-b\right) - \sqrt{b \cdot b - 4 \cdot \left(a \cdot c\right)}}{2 \cdot a}\begin{array}{l}
\mathbf{if}\;b \le -6239046376.848015:\\
\;\;\;\;\frac{-c}{b}\\
\mathbf{elif}\;b \le -2.396838732613953 \cdot 10^{-306}:\\
\;\;\;\;\frac{\frac{\frac{1}{2}}{a} \cdot \left(\left(c \cdot 4\right) \cdot a\right)}{\left(-b\right) + \sqrt{b \cdot b - 4 \cdot \left(c \cdot a\right)}}\\
\mathbf{elif}\;b \le 9.179168538250646 \cdot 10^{+63}:\\
\;\;\;\;\frac{\frac{1}{2}}{a} \cdot \left(\left(-b\right) - \sqrt{b \cdot b - 4 \cdot \left(c \cdot a\right)}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{-b}{a}\\
\end{array}double f(double a, double b, double c) {
double r2013272 = b;
double r2013273 = -r2013272;
double r2013274 = r2013272 * r2013272;
double r2013275 = 4.0;
double r2013276 = a;
double r2013277 = c;
double r2013278 = r2013276 * r2013277;
double r2013279 = r2013275 * r2013278;
double r2013280 = r2013274 - r2013279;
double r2013281 = sqrt(r2013280);
double r2013282 = r2013273 - r2013281;
double r2013283 = 2.0;
double r2013284 = r2013283 * r2013276;
double r2013285 = r2013282 / r2013284;
return r2013285;
}
double f(double a, double b, double c) {
double r2013286 = b;
double r2013287 = -6239046376.848015;
bool r2013288 = r2013286 <= r2013287;
double r2013289 = c;
double r2013290 = -r2013289;
double r2013291 = r2013290 / r2013286;
double r2013292 = -2.396838732613953e-306;
bool r2013293 = r2013286 <= r2013292;
double r2013294 = 0.5;
double r2013295 = a;
double r2013296 = r2013294 / r2013295;
double r2013297 = 4.0;
double r2013298 = r2013289 * r2013297;
double r2013299 = r2013298 * r2013295;
double r2013300 = r2013296 * r2013299;
double r2013301 = -r2013286;
double r2013302 = r2013286 * r2013286;
double r2013303 = r2013289 * r2013295;
double r2013304 = r2013297 * r2013303;
double r2013305 = r2013302 - r2013304;
double r2013306 = sqrt(r2013305);
double r2013307 = r2013301 + r2013306;
double r2013308 = r2013300 / r2013307;
double r2013309 = 9.179168538250646e+63;
bool r2013310 = r2013286 <= r2013309;
double r2013311 = r2013301 - r2013306;
double r2013312 = r2013296 * r2013311;
double r2013313 = r2013301 / r2013295;
double r2013314 = r2013310 ? r2013312 : r2013313;
double r2013315 = r2013293 ? r2013308 : r2013314;
double r2013316 = r2013288 ? r2013291 : r2013315;
return r2013316;
}




Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.7 |
|---|---|
| Target | 20.7 |
| Herbie | 9.0 |
if b < -6239046376.848015Initial program 55.7
Taylor expanded around -inf 5.1
Simplified5.1
if -6239046376.848015 < b < -2.396838732613953e-306Initial program 27.2
rmApplied clear-num27.2
rmApplied flip--27.3
Applied associate-/r/27.3
Applied associate-/r*27.3
Simplified17.2
if -2.396838732613953e-306 < b < 9.179168538250646e+63Initial program 9.1
rmApplied clear-num9.3
rmApplied associate-/r/9.3
Simplified9.3
if 9.179168538250646e+63 < b Initial program 37.8
rmApplied clear-num37.9
Taylor expanded around 0 6.1
Simplified6.1
Final simplification9.0
herbie shell --seed 2019146
(FPCore (a b c)
:name "quadm (p42, negative)"
:herbie-target
(if (< b 0) (/ c (* a (/ (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))) (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))
(/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))