\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 -2.569494919068124572690421335939486791404 \cdot 10^{-64}:\\
\;\;\;\;-1 \cdot \frac{c}{b}\\
\mathbf{elif}\;b \le 2.8653816703769607550753035783606354728 \cdot 10^{117}:\\
\;\;\;\;\frac{\left(-b\right) - \sqrt{b \cdot b - \left(4 \cdot c\right) \cdot a}}{a} \cdot \frac{1}{2}\\
\mathbf{else}:\\
\;\;\;\;-1 \cdot \frac{b}{a}\\
\end{array}double f(double a, double b, double c) {
double r4070563 = b;
double r4070564 = -r4070563;
double r4070565 = r4070563 * r4070563;
double r4070566 = 4.0;
double r4070567 = a;
double r4070568 = c;
double r4070569 = r4070567 * r4070568;
double r4070570 = r4070566 * r4070569;
double r4070571 = r4070565 - r4070570;
double r4070572 = sqrt(r4070571);
double r4070573 = r4070564 - r4070572;
double r4070574 = 2.0;
double r4070575 = r4070574 * r4070567;
double r4070576 = r4070573 / r4070575;
return r4070576;
}
double f(double a, double b, double c) {
double r4070577 = b;
double r4070578 = -2.5694949190681246e-64;
bool r4070579 = r4070577 <= r4070578;
double r4070580 = -1.0;
double r4070581 = c;
double r4070582 = r4070581 / r4070577;
double r4070583 = r4070580 * r4070582;
double r4070584 = 2.865381670376961e+117;
bool r4070585 = r4070577 <= r4070584;
double r4070586 = -r4070577;
double r4070587 = r4070577 * r4070577;
double r4070588 = 4.0;
double r4070589 = r4070588 * r4070581;
double r4070590 = a;
double r4070591 = r4070589 * r4070590;
double r4070592 = r4070587 - r4070591;
double r4070593 = sqrt(r4070592);
double r4070594 = r4070586 - r4070593;
double r4070595 = r4070594 / r4070590;
double r4070596 = 1.0;
double r4070597 = 2.0;
double r4070598 = r4070596 / r4070597;
double r4070599 = r4070595 * r4070598;
double r4070600 = r4070577 / r4070590;
double r4070601 = r4070580 * r4070600;
double r4070602 = r4070585 ? r4070599 : r4070601;
double r4070603 = r4070579 ? r4070583 : r4070602;
return r4070603;
}




Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.8 |
|---|---|
| Target | 20.9 |
| Herbie | 10.2 |
if b < -2.5694949190681246e-64Initial program 53.0
Taylor expanded around -inf 9.1
if -2.5694949190681246e-64 < b < 2.865381670376961e+117Initial program 13.1
rmApplied clear-num13.2
rmApplied *-un-lft-identity13.2
Applied times-frac13.2
Applied add-cube-cbrt13.2
Applied times-frac13.2
Simplified13.2
Simplified13.1
if 2.865381670376961e+117 < b Initial program 52.1
rmApplied clear-num52.2
Taylor expanded around 0 3.1
Final simplification10.2
herbie shell --seed 2019169
(FPCore (a b c)
:name "quadm (p42, negative)"
:herbie-target
(if (< b 0.0) (/ c (* a (/ (+ (- b) (sqrt (- (* b b) (* 4.0 (* a c))))) (* 2.0 a)))) (/ (- (- b) (sqrt (- (* b b) (* 4.0 (* a c))))) (* 2.0 a)))
(/ (- (- b) (sqrt (- (* b b) (* 4.0 (* a c))))) (* 2.0 a)))