


Bits error versus a



Bits error versus b



Bits error versus c
Initial program 43.6
Simplified43.6
rmApplied flip--43.6
Applied associate-/l/43.6
Simplified0.4
rmApplied times-frac0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019002 +o rules:numerics
(FPCore (a b c)
:name "Quadratic roots, medium range"
: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 a) c)))) (* 2 a)))
Time bar (total: 17.6s)Debug log
1 calls. Slowest were:
| 92.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 43.9b
Found 4 expressions with local error:
| 39.7b | (- (sqrt (fma (* a -4) c (* b b))) b) |
| 0.2b | (/ (- (sqrt (fma (* a -4) c (* b b))) b) (* 2 a)) |
| 0.1b | (fma (* a -4) c (* b b)) |
| 0.1b | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 10.0ms | (/ (- (sqrt (fma (* a -4) c (* b b))) b) (* 2 a)) |
| 5.0ms | (- (sqrt (fma (* a -4) c (* b b))) b) |
| 1.0ms | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 106.0ms | (/ (- (sqrt (fma (* a -4) c (* b b))) b) (* 2 a)) |
| 68.0ms | (- (sqrt (fma (* a -4) c (* b b))) b) |
| 48.0ms | (sqrt (fma (* a -4) c (* b b))) |
| 29.0ms | (fma (* a -4) c (* b b)) |
34 calls. Slowest were:
| 780.0ms | (- (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b)))) (* b b)) |
| 125.0ms | (- (pow (sqrt (fma (* a -4) c (* b b))) 3) (pow b 3)) |
| 62.0ms | (/ (* (* (- (sqrt (fma (* a -4) c (* b b))) b) (- (sqrt (fma (* a -4) c (* b b))) b)) (- (sqrt (fma (* a -4) c (* b b))) b)) (* (* (* 2 a) (* 2 a)) (* 2 a))) |
3 alts after pruning (2 fresh and 1 done)
Merged error: 0.4b
Found 4 expressions with local error:
| 0.3b | (/ (* (* a -4) c) (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b))) |
| 0.2b | (+ (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)) |
4 calls. Slowest were:
| 11.0ms | (/ (* (* a -4) c) (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b))) |
| 8.0ms | (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 4.0ms | (+ (sqrt (fma (* a -4) c (* b b))) b) |
4 calls. Slowest were:
| 144.0ms | (/ (* (* a -4) c) (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b))) |
| 88.0ms | (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 67.0ms | (+ (sqrt (fma (* a -4) c (* b b))) b) |
| 32.0ms | (fma (* a -4) c (* b b)) |
37 calls. Slowest were:
| 683.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)))) |
| 682.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)))) |
| 648.0ms | (- (log (* (* a -4) c)) (+ (log (* 2 a)) (log (+ (sqrt (fma (* a -4) c (* b b))) b)))) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.2b | (+ (sqrt (fma (* a -4) c (* b b))) b) |
| 0.2b | (/ c (+ (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))) |
4 calls. Slowest were:
| 9.0ms | (+ (sqrt (fma (* a -4) c (* b b))) b) |
| 5.0ms | (/ c (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 2.0ms | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 125.0ms | (/ c (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 56.0ms | (sqrt (fma (* a -4) c (* b b))) |
| 54.0ms | (+ (sqrt (fma (* a -4) c (* b b))) b) |
| 41.0ms | (fma (* a -4) c (* b b)) |
26 calls. Slowest were:
| 219.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))) |
| 53.0ms | (- (pow b 2) (* 4 (* a c))) |
| 50.0ms | (- (pow b 2) (* 4 (* a c))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.2b | (fma (sqrt (sqrt (fma (* a -4) c (* b b)))) (sqrt (sqrt (fma (* a -4) c (* b b)))) b) |
| 0.2b | (/ c (fma (sqrt (sqrt (fma (* a -4) c (* b b)))) (sqrt (sqrt (fma (* a -4) c (* b b)))) b)) |
| 0.1b | (fma (* a -4) c (* b b)) |
| 0.1b | (fma (* a -4) c (* b b)) |
4 calls. Slowest were:
| 1.0ms | (/ c (fma (sqrt (sqrt (fma (* a -4) c (* b b)))) (sqrt (sqrt (fma (* a -4) c (* b b)))) b)) |
| 0.0ms | (fma (sqrt (sqrt (fma (* a -4) c (* b b)))) (sqrt (sqrt (fma (* a -4) c (* b b)))) b) |
| 0.0ms | (fma (* a -4) c (* b b)) |
4 calls. Slowest were:
| 332.0ms | (fma (sqrt (sqrt (fma (* a -4) c (* b b)))) (sqrt (sqrt (fma (* a -4) c (* b b)))) b) |
| 297.0ms | (/ c (fma (sqrt (sqrt (fma (* a -4) c (* b b)))) (sqrt (sqrt (fma (* a -4) c (* b b)))) b)) |
| 50.0ms | (fma (* a -4) c (* b b)) |
| 34.0ms | (fma (* a -4) c (* b b)) |
15 calls. Slowest were:
| 321.0ms | (+ (pow (exp (* 1/4 (- (log -4) (+ (log (/ 1 a)) (log (/ 1 c)))))) 2) b) |
| 278.0ms | (- (/ c b) (+ (* 1/4 (/ (pow (exp (* 1/4 (- (log -4) (+ (log (/ -1 c)) (log (/ -1 a)))))) 4) (* a b))) (/ (* c (pow (exp (* 1/4 (- (log -4) (+ (log (/ -1 c)) (log (/ -1 a)))))) 2)) (pow b 2)))) |
| 156.0ms | (- (/ c b) (+ (/ (* (pow (exp (* 1/4 (- (log -4) (+ (log (/ 1 a)) (log (/ 1 c)))))) 2) c) (pow b 2)) (* 1/4 (/ (pow (exp (* 1/4 (- (log -4) (+ (log (/ 1 a)) (log (/ 1 c)))))) 4) (* a b))))) |
7 alts after pruning (5 fresh and 2 done)
Merged error: 0b
0% (0.1b remaining)
Error of 0.2b against oracle of 0.0b and baseline of 0.2b