\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 -1.970010565552108757188050455448622102575 \cdot 10^{58}:\\
\;\;\;\;-1 \cdot \frac{c}{b}\\
\mathbf{elif}\;b \le -1.341090161490938310878248498686683235407 \cdot 10^{-308}:\\
\;\;\;\;\frac{\frac{\left(a \cdot c\right) \cdot 4}{\sqrt{b \cdot b - \left(a \cdot c\right) \cdot 4} - b}}{a \cdot 2}\\
\mathbf{elif}\;b \le 3.628799960716311990444092539387346352569 \cdot 10^{50}:\\
\;\;\;\;\frac{\left(-b\right) - \sqrt{b \cdot b - \left(a \cdot c\right) \cdot 4}}{a \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{b}{a} \cdot -1\\
\end{array}double f(double a, double b, double c) {
double r4704593 = b;
double r4704594 = -r4704593;
double r4704595 = r4704593 * r4704593;
double r4704596 = 4.0;
double r4704597 = a;
double r4704598 = c;
double r4704599 = r4704597 * r4704598;
double r4704600 = r4704596 * r4704599;
double r4704601 = r4704595 - r4704600;
double r4704602 = sqrt(r4704601);
double r4704603 = r4704594 - r4704602;
double r4704604 = 2.0;
double r4704605 = r4704604 * r4704597;
double r4704606 = r4704603 / r4704605;
return r4704606;
}
double f(double a, double b, double c) {
double r4704607 = b;
double r4704608 = -1.9700105655521088e+58;
bool r4704609 = r4704607 <= r4704608;
double r4704610 = -1.0;
double r4704611 = c;
double r4704612 = r4704611 / r4704607;
double r4704613 = r4704610 * r4704612;
double r4704614 = -1.3410901614909383e-308;
bool r4704615 = r4704607 <= r4704614;
double r4704616 = a;
double r4704617 = r4704616 * r4704611;
double r4704618 = 4.0;
double r4704619 = r4704617 * r4704618;
double r4704620 = r4704607 * r4704607;
double r4704621 = r4704620 - r4704619;
double r4704622 = sqrt(r4704621);
double r4704623 = r4704622 - r4704607;
double r4704624 = r4704619 / r4704623;
double r4704625 = 2.0;
double r4704626 = r4704616 * r4704625;
double r4704627 = r4704624 / r4704626;
double r4704628 = 3.628799960716312e+50;
bool r4704629 = r4704607 <= r4704628;
double r4704630 = -r4704607;
double r4704631 = r4704630 - r4704622;
double r4704632 = r4704631 / r4704626;
double r4704633 = r4704607 / r4704616;
double r4704634 = r4704633 * r4704610;
double r4704635 = r4704629 ? r4704632 : r4704634;
double r4704636 = r4704615 ? r4704627 : r4704635;
double r4704637 = r4704609 ? r4704613 : r4704636;
return r4704637;
}




Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.8 |
|---|---|
| Target | 20.7 |
| Herbie | 9.1 |
if b < -1.9700105655521088e+58Initial program 57.5
Taylor expanded around -inf 3.4
if -1.9700105655521088e+58 < b < -1.3410901614909383e-308Initial program 29.3
rmApplied flip--29.4
Simplified16.6
Simplified16.6
if -1.3410901614909383e-308 < b < 3.628799960716312e+50Initial program 9.7
if 3.628799960716312e+50 < b Initial program 38.2
rmApplied clear-num38.3
Taylor expanded around 0 6.3
Final simplification9.1
herbie shell --seed 2019179
(FPCore (a b c)
:name "The quadratic formula (r2)"
: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)))