


Bits error versus a



Bits error versus b



Bits error versus c
Initial program 52.6
Initial simplification52.6
rmApplied flip--52.6
Applied associate-/l/52.6
Simplified0.4
rmApplied times-frac0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2018365 +o rules:numerics
(FPCore (a b c)
:name "Quadratic roots, wide range"
:pre (and (< 4.930380657631324e-32 a 2.028240960365167e+31) (< 4.930380657631324e-32 b 2.028240960365167e+31) (< 4.930380657631324e-32 c 2.028240960365167e+31))
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
Time bar (total: 16.2s)Debug log
2 alts after pruning (2 fresh and 0 done)
Merged error: 52.4b
Found 4 expressions with local error:
| 25.3b | (- (sqrt (fma (* a -4) c (* b b))) b) |
| 0.1b | (fma (* a -4) c (* b b)) |
| 0.1b | (/ (- (sqrt (fma (* a -4) c (* b b))) b) (* 2 a)) |
| 0.0b | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 10.0ms | (/ (- (sqrt (fma (* a -4) c (* b b))) b) (* 2 a)) |
| 7.0ms | (- (sqrt (fma (* a -4) c (* b b))) b) |
| 1.0ms | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 214.0ms | (- (sqrt (fma (* a -4) c (* b b))) b) |
| 126.0ms | (/ (- (sqrt (fma (* a -4) c (* b b))) b) (* 2 a)) |
| 50.0ms | (fma (* a -4) c (* b b)) |
| 45.0ms | (sqrt (fma (* a -4) c (* b b))) |
34 calls. Slowest were:
| 643.0ms | (- (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b)))) (* b b)) |
| 439.0ms | (/ (* (cbrt (- (sqrt (fma (* a -4) c (* b b))) b)) (cbrt (- (sqrt (fma (* a -4) c (* b b))) b))) 2) |
| 109.0ms | (- (pow (sqrt (fma (* a -4) c (* b b))) 3) (pow b 3)) |
3 alts after pruning (3 fresh and 0 done)
Merged error: 0.3b
Found 4 expressions with local error:
| 0.3b | (/ (* (* a -4) c) (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b))) |
| 0.1b | (fma (* a -4) c (* b b)) |
| 0.1b | (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 0.1b | (+ (sqrt (fma (* a -4) c (* b b))) b) |
4 calls. Slowest were:
| 19.0ms | (/ (* (* a -4) c) (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b))) |
| 16.0ms | (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 7.0ms | (+ (sqrt (fma (* a -4) c (* b b))) b) |
4 calls. Slowest were:
| 135.0ms | (/ (* (* a -4) c) (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b))) |
| 98.0ms | (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 58.0ms | (+ (sqrt (fma (* a -4) c (* b b))) b) |
| 28.0ms | (fma (* a -4) c (* b b)) |
37 calls. Slowest were:
| 710.0ms | (/ (* (* (* (* a -4) c) (* (* a -4) c)) (* (* a -4) c)) (* (* (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b)) (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b))) (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b)))) |
| 709.0ms | (- (log (* (* a -4) c)) (+ (log (* 2 a)) (log (+ (sqrt (fma (* a -4) c (* b b))) b)))) |
| 645.0ms | (/ (* (* (* (* a -4) c) (* (* a -4) c)) (* (* a -4) c)) (* (* (* (* 2 a) (* 2 a)) (* 2 a)) (* (* (+ (sqrt (fma (* a -4) c (* b b))) b) (+ (sqrt (fma (* a -4) c (* b b))) b)) (+ (sqrt (fma (* a -4) c (* b b))) b)))) |
4 alts after pruning (4 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.1b | (fma (* a -4) c (* b b)) |
| 0.1b | (+ (sqrt (fma (* a -4) c (* b b))) b) |
| 0.0b | (/ c (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 0.0b | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 4.0ms | (+ (sqrt (fma (* a -4) c (* b b))) b) |
| 2.0ms | (/ c (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 1.0ms | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 146.0ms | (/ c (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 80.0ms | (+ (sqrt (fma (* a -4) c (* b b))) b) |
| 46.0ms | (sqrt (fma (* a -4) c (* b b))) |
| 42.0ms | (fma (* a -4) c (* b b)) |
26 calls. Slowest were:
| 300.0ms | (+ (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b)))) (- (* b b) (* (sqrt (fma (* a -4) c (* b b))) b))) |
| 49.0ms | (- (pow b 2) (* 4 (* a c))) |
| 37.0ms | (- (pow b 2) (* 4 (* a c))) |
4 alts after pruning (3 fresh and 1 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.3b | (/ (* (* a -4) c) (+ (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) (* (* 2 a) b))) |
| 0.1b | (fma (* a -4) c (* b b)) |
| 0.1b | (+ (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) (* (* 2 a) b)) |
| 0.0b | (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) |
4 calls. Slowest were:
| 19.0ms | (/ (* (* a -4) c) (+ (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) (* (* 2 a) b))) |
| 11.0ms | (+ (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) (* (* 2 a) b)) |
| 8.0ms | (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) |
4 calls. Slowest were:
| 220.0ms | (/ (* (* a -4) c) (+ (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) (* (* 2 a) b))) |
| 80.0ms | (+ (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) (* (* 2 a) b)) |
| 76.0ms | (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) |
| 38.0ms | (fma (* a -4) c (* b b)) |
42 calls. Slowest were:
| 778.0ms | (* (* (* (* 2 a) (* 2 a)) (* 2 a)) (* (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b)))) (sqrt (fma (* a -4) c (* b b))))) |
| 486.0ms | (/ (* (* (* (* a -4) c) (* (* a -4) c)) (* (* a -4) c)) (* (* (+ (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) (* (* 2 a) b)) (+ (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) (* (* 2 a) b))) (+ (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) (* (* 2 a) b)))) |
| 209.0ms | (- (log (* (* a -4) c)) (log (+ (* (* 2 a) (sqrt (fma (* a -4) c (* b b)))) (* (* 2 a) b)))) |
4 alts after pruning (2 fresh and 2 done)
Merged error: 0b
0% (0.1b remaining)
Error of 0.1b against oracle of 0.0b and baseline of 0.1b