


Bits error versus a



Bits error versus b



Bits error versus c
Initial program 28.3
Simplified28.3
rmApplied flip--28.3
Applied associate-/l/28.3
Simplified0.5
rmApplied times-frac0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2019005 +o rules:numerics
(FPCore (a b c)
:name "Quadratic roots, narrow range"
:pre (and (< 1.0536712127723509e-08 a 94906265.62425156) (< 1.0536712127723509e-08 b 94906265.62425156) (< 1.0536712127723509e-08 c 94906265.62425156))
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
Time bar (total: 21.0s)Debug log
1 calls. Slowest were:
| 66.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 28.6b
Found 4 expressions with local error:
| 28.4b | (- (sqrt (fma (* a -4) c (* b b))) b) |
| 0.3b | (/ (- (sqrt (fma (* a -4) c (* b b))) b) (* 2 a)) |
| 0.2b | (fma (* a -4) c (* b b)) |
| 0.1b | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 21.0ms | (/ (- (sqrt (fma (* a -4) c (* b b))) b) (* 2 a)) |
| 11.0ms | (- (sqrt (fma (* a -4) c (* b b))) b) |
| 2.0ms | (sqrt (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 116.0ms | (/ (- (sqrt (fma (* a -4) c (* b b))) b) (* 2 a)) |
| 75.0ms | (- (sqrt (fma (* a -4) c (* b b))) b) |
| 46.0ms | (sqrt (fma (* a -4) c (* b b))) |
| 32.0ms | (fma (* a -4) c (* b b)) |
34 calls. Slowest were:
| 768.0ms | (- (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b)))) (* b b)) |
| 102.0ms | (- (pow (sqrt (fma (* a -4) c (* b b))) 3) (pow b 3)) |
| 75.0ms | (- (log (- (sqrt (fma (* a -4) c (* b b))) b)) (log (* 2 a))) |
6 alts after pruning (6 fresh and 0 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 | (fma (* a -4) c (* b b)) |
| 0.2b | (+ (sqrt (fma (* a -4) c (* b b))) b) |
| 0.2b | (* (* 2 a) (+ (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))) |
| 8.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:
| 131.0ms | (/ (* (* a -4) c) (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b))) |
| 100.0ms | (* (* 2 a) (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 55.0ms | (+ (sqrt (fma (* a -4) c (* b b))) b) |
| 35.0ms | (fma (* a -4) c (* b b)) |
37 calls. Slowest were:
| 675.0ms | (- (log (* (* a -4) c)) (+ (log (* 2 a)) (log (+ (sqrt (fma (* a -4) c (* b b))) b)))) |
| 624.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)))) |
| 544.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)))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.2b | (fma (* a -4) c (* b b)) |
| 0.2b | (+ (sqrt (fma (* a -4) c (* b b))) b) |
| 0.2b | (/ c (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 0.1b | (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:
| 124.0ms | (/ c (+ (sqrt (fma (* a -4) c (* b b))) b)) |
| 70.0ms | (+ (sqrt (fma (* a -4) c (* b b))) b) |
| 47.0ms | (sqrt (fma (* a -4) c (* b b))) |
| 35.0ms | (fma (* a -4) c (* b b)) |
26 calls. Slowest were:
| 242.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))) |
| 40.0ms | (- (pow b 2) (* 4 (* a c))) |
| 37.0ms | (- (pow b 2) (* 4 (* a c))) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.5b | (cbrt (* (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) (fma (* a -4) c (* b b)))) |
| 0.4b | (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) |
| 0.4b | (* (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) (fma (* a -4) c (* b b))) |
| 0.2b | (fma (* a -4) c (* b b)) |
4 calls. Slowest were:
| 12.0ms | (* (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) (fma (* a -4) c (* b b))) |
| 6.0ms | (cbrt (* (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) (fma (* a -4) c (* b b)))) |
| 4.0ms | (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) |
4 calls. Slowest were:
| 98.0ms | (* (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) (fma (* a -4) c (* b b))) |
| 51.0ms | (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) |
| 37.0ms | (fma (* a -4) c (* b b)) |
| 37.0ms | (cbrt (* (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) (fma (* a -4) c (* b b)))) |
49 calls. Slowest were:
| 1.0s | (* (* (* (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) (fma (* a -4) c (* b b))) (* (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) (fma (* a -4) c (* b b)))) (* (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) (fma (* a -4) c (* b b)))) |
| 784.0ms | (* (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) (* b b)) |
| 446.0ms | (* (* (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) (fma (* a -4) c (* b b))) (* (* (fma (* a -4) c (* b b)) (fma (* a -4) c (* b b))) (fma (* a -4) c (* b b)))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 0b
0% (0.2b remaining)
Error of 0.3b against oracle of 0.0b and baseline of 0.3b