\frac{\left(-b\right) + \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c}}{3 \cdot a}\frac{1}{-\frac{b + \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c}}{c}}double code(double a, double b, double c) {
return ((-b + sqrt(((b * b) - ((3.0 * a) * c)))) / (3.0 * a));
}
double code(double a, double b, double c) {
return (1.0 / -((b + sqrt(((b * b) - ((3.0 * a) * c)))) / c));
}



Bits error versus a



Bits error versus b



Bits error versus c
Results
Initial program 28.6
rmApplied flip-+28.6
Simplified0.6
rmApplied associate-*r*0.4
rmApplied *-commutative0.4
Applied associate-/l*0.5
Applied associate-/l/0.3
Simplified0.3
rmApplied clear-num0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2020078
(FPCore (a b c)
:name "Cubic critical, 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) (* (* 3 a) c)))) (* 3 a)))