\frac{\left(-b\right) + \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c}}{3 \cdot a}\frac{\frac{-3 \cdot a}{-1} \cdot \frac{c}{\left(-b\right) - \sqrt{b \cdot b - \left(3 \cdot a\right) \cdot c}}}{3 \cdot a}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 (((-(3.0 * a) / -1.0) * (c / (-b - sqrt(((b * b) - ((3.0 * a) * c)))))) / (3.0 * a));
}



Bits error versus a



Bits error versus b



Bits error versus c
Results
Initial program 28.4
rmApplied flip-+28.3
Simplified0.6
rmApplied frac-2neg0.6
Simplified0.5
rmApplied *-un-lft-identity0.5
Applied distribute-lft-neg-in0.5
Applied distribute-lft-neg-in0.5
Applied times-frac0.3
Final simplification0.3
herbie shell --seed 2020060
(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)))