\frac{\left(-b\right) + \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}}{2 \cdot a}\frac{\frac{\frac{\frac{\sqrt{{\left(b \cdot b - \left(4 \cdot a\right) \cdot c\right)}^{6}} \cdot \sqrt{{\left(b \cdot b - \left(4 \cdot a\right) \cdot c\right)}^{6}} - {b}^{12}}{{\left(b \cdot b - \left(4 \cdot a\right) \cdot c\right)}^{3} + {b}^{6}}}{\left(b \cdot b - \left(4 \cdot a\right) \cdot c\right) \cdot \left(b \cdot b - \left(4 \cdot a\right) \cdot c\right) + \left({b}^{4} + \left(b \cdot b\right) \cdot \left(b \cdot b - \left(4 \cdot a\right) \cdot c\right)\right)}}{b + \sqrt{b \cdot b - \left(4 \cdot a\right) \cdot c}}}{a \cdot 2}(FPCore (a b c) :precision binary64 (/ (+ (- b) (sqrt (- (* b b) (* (* 4.0 a) c)))) (* 2.0 a)))
(FPCore (a b c)
:precision binary64
(/
(/
(/
(/
(-
(*
(sqrt (pow (- (* b b) (* (* 4.0 a) c)) 6.0))
(sqrt (pow (- (* b b) (* (* 4.0 a) c)) 6.0)))
(pow b 12.0))
(+ (pow (- (* b b) (* (* 4.0 a) c)) 3.0) (pow b 6.0)))
(+
(* (- (* b b) (* (* 4.0 a) c)) (- (* b b) (* (* 4.0 a) c)))
(+ (pow b 4.0) (* (* b b) (- (* b b) (* (* 4.0 a) c))))))
(+ b (sqrt (- (* b b) (* (* 4.0 a) c)))))
(* a 2.0)))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 (((((sqrt(pow(((b * b) - ((4.0 * a) * c)), 6.0)) * sqrt(pow(((b * b) - ((4.0 * a) * c)), 6.0))) - pow(b, 12.0)) / (pow(((b * b) - ((4.0 * a) * c)), 3.0) + pow(b, 6.0))) / ((((b * b) - ((4.0 * a) * c)) * ((b * b) - ((4.0 * a) * c))) + (pow(b, 4.0) + ((b * b) * ((b * b) - ((4.0 * a) * c)))))) / (b + sqrt((b * b) - ((4.0 * a) * c)))) / (a * 2.0);
}



Bits error versus a



Bits error versus b



Bits error versus c
Results
Initial program 43.8
Simplified43.8
rmApplied flip--_binary6443.9
Simplified43.3
Simplified43.3
rmApplied flip3--_binary6443.3
Simplified42.8
Simplified42.8
rmApplied flip--_binary6442.8
Simplified42.6
rmApplied add-sqr-sqrt_binary6442.6
Final simplification42.6
herbie shell --seed 2020299
(FPCore (a b c)
:name "Quadratic roots, 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) (* (* 4.0 a) c)))) (* 2.0 a)))