\frac{\left(-b\right) + \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}}{2 \cdot a}\frac{\frac{\left({b}^{2} - {b}^{2}\right) + 4 \cdot \left(a \cdot c\right)}{\left(-b\right) - \sqrt{\left(\left(\sqrt{4} \cdot \sqrt{a}\right) \cdot \sqrt{c} + b\right) \cdot \left(b - \left(\sqrt{4} \cdot \sqrt{a}\right) \cdot \sqrt{c}\right)}}}{2 \cdot a}double code(double a, double b, double c) {
return ((double) (((double) (((double) -(b)) + ((double) sqrt(((double) (((double) (b * b)) - ((double) (((double) (4.0 * a)) * c)))))))) / ((double) (2.0 * a))));
}
double code(double a, double b, double c) {
return ((double) (((double) (((double) (((double) (((double) pow(b, 2.0)) - ((double) pow(b, 2.0)))) + ((double) (4.0 * ((double) (a * c)))))) / ((double) (((double) -(b)) - ((double) sqrt(((double) (((double) (((double) (((double) (((double) sqrt(4.0)) * ((double) sqrt(a)))) * ((double) sqrt(c)))) + b)) * ((double) (b - ((double) (((double) (((double) sqrt(4.0)) * ((double) sqrt(a)))) * ((double) sqrt(c)))))))))))))) / ((double) (2.0 * a))));
}



Bits error versus a



Bits error versus b



Bits error versus c
Results
Initial program 44.0
rmApplied flip-+44.0
Simplified0.4
rmApplied add-sqr-sqrt0.4
Applied add-sqr-sqrt0.4
Applied add-sqr-sqrt0.4
Applied unswap-sqr0.4
Applied unswap-sqr0.4
Applied difference-of-squares0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020114
(FPCore (a b c)
:name "Quadratic roots, medium range"
:precision binary64
:pre (and (< 1.11022e-16 a 9.0072e+15) (< 1.11022e-16 b 9.0072e+15) (< 1.11022e-16 c 9.0072e+15))
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))