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



Bits error versus a



Bits error versus b



Bits error versus c
Results
Initial program 44.2
Simplified44.2
rmApplied flip--44.2
Simplified0.5
Simplified0.5
rmApplied *-un-lft-identity0.5
Applied times-frac0.5
Applied associate-/l*0.5
Simplified0.4
rmApplied div-inv0.4
Applied times-frac0.5
Simplified0.4
Final simplification0.4
herbie shell --seed 2020196
(FPCore (a b c)
:name "Cubic critical, medium range"
:precision binary64
:pre (and (< 1.1102230246251565e-16 a 9007199254740992.0) (< 1.1102230246251565e-16 b 9007199254740992.0) (< 1.1102230246251565e-16 c 9007199254740992.0))
(/ (+ (- b) (sqrt (- (* b b) (* (* 3.0 a) c)))) (* 3.0 a)))