\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.6239127264630285 \cdot 10^{-63}:\\
\;\;\;\;\frac{-2 \cdot \frac{c}{b}}{2}\\
\mathbf{elif}\;b \le 7.052614559736995 \cdot 10^{+62}:\\
\;\;\;\;\frac{\frac{1}{a} \cdot \left(\left(-b\right) - \sqrt{\mathsf{fma}\left(-4 \cdot a, c, b \cdot b\right)}\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\frac{c}{b} - \frac{b}{a}\right) \cdot 2}{2}\\
\end{array}double f(double a, double b, double c) {
double r3101768 = b;
double r3101769 = -r3101768;
double r3101770 = r3101768 * r3101768;
double r3101771 = 4.0;
double r3101772 = a;
double r3101773 = c;
double r3101774 = r3101772 * r3101773;
double r3101775 = r3101771 * r3101774;
double r3101776 = r3101770 - r3101775;
double r3101777 = sqrt(r3101776);
double r3101778 = r3101769 - r3101777;
double r3101779 = 2.0;
double r3101780 = r3101779 * r3101772;
double r3101781 = r3101778 / r3101780;
return r3101781;
}
double f(double a, double b, double c) {
double r3101782 = b;
double r3101783 = -1.6239127264630285e-63;
bool r3101784 = r3101782 <= r3101783;
double r3101785 = -2.0;
double r3101786 = c;
double r3101787 = r3101786 / r3101782;
double r3101788 = r3101785 * r3101787;
double r3101789 = 2.0;
double r3101790 = r3101788 / r3101789;
double r3101791 = 7.052614559736995e+62;
bool r3101792 = r3101782 <= r3101791;
double r3101793 = 1.0;
double r3101794 = a;
double r3101795 = r3101793 / r3101794;
double r3101796 = -r3101782;
double r3101797 = -4.0;
double r3101798 = r3101797 * r3101794;
double r3101799 = r3101782 * r3101782;
double r3101800 = fma(r3101798, r3101786, r3101799);
double r3101801 = sqrt(r3101800);
double r3101802 = r3101796 - r3101801;
double r3101803 = r3101795 * r3101802;
double r3101804 = r3101803 / r3101789;
double r3101805 = r3101782 / r3101794;
double r3101806 = r3101787 - r3101805;
double r3101807 = r3101806 * r3101789;
double r3101808 = r3101807 / r3101789;
double r3101809 = r3101792 ? r3101804 : r3101808;
double r3101810 = r3101784 ? r3101790 : r3101809;
return r3101810;
}




Bits error versus a




Bits error versus b




Bits error versus c
| Original | 32.8 |
|---|---|
| Target | 20.1 |
| Herbie | 10.3 |
if b < -1.6239127264630285e-63Initial program 52.2
Simplified52.3
Taylor expanded around -inf 8.6
if -1.6239127264630285e-63 < b < 7.052614559736995e+62Initial program 13.9
Simplified14.0
rmApplied clear-num14.1
rmApplied div-inv14.1
rmApplied un-div-inv14.1
Applied associate-/r/14.1
if 7.052614559736995e+62 < b Initial program 38.1
Simplified38.0
Taylor expanded around inf 4.7
Simplified4.7
Final simplification10.3
herbie shell --seed 2019158 +o rules:numerics
(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)))