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



Bits error versus a



Bits error versus b



Bits error versus c
Results
Initial program 28.4
rmApplied flip-+28.4
Simplified0.4
rmApplied *-commutative0.4
Applied associate-*r*0.4
Applied associate-/l*0.4
rmApplied frac-2neg0.4
Applied associate-/l/0.3
Simplified0.3
rmApplied flip3--0.3
Applied sqrt-div0.3
Final simplification0.3
herbie shell --seed 2020078
(FPCore (a b c)
:name "Quadratic roots, narrow range"
:precision binary64
:pre (and (< 1.0536712127723509e-08 a 94906265.62425156) (< 1.0536712127723509e-08 b 94906265.62425156) (< 1.0536712127723509e-08 c 94906265.62425156))
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))