


Bits error versus a



Bits error versus b



Bits error versus c
Initial program 28.7
rmApplied flip-+28.8
Applied associate-/l/28.8
Simplified0.5
rmApplied associate-/r*0.3
Simplified0.3
rmApplied times-frac0.3
Applied associate-/l*0.4
Simplified0.4
rmApplied associate-/r/0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2018360 +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: 1.4m)Debug log
1 alts after pruning (1 fresh and 0 done)
Merged error: 29.1b
Found 4 expressions with local error:
| 28.9b | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 0.2b | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 0.2b | (- (* b b) (* (* 4 a) c)) |
| 0.1b | (sqrt (- (* b b) (* (* 4 a) c))) |
4 calls. Slowest were:
| 17.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 7.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 6.0ms | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 77.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 45.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 43.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 30.0ms | (- (* b b) (* (* 4 a) c)) |
40 calls. Slowest were:
| 379.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 311.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 308.0ms | (+ (pow (- b) 3) (pow (sqrt (- (* b b) (* (* 4 a) c))) 3)) |
4 alts after pruning (4 fresh and 0 done)
Merged error: 0.4b
Found 4 expressions with local error:
| 0.4b | (/ (* 4 (* c a)) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 0.2b | (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 0.2b | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 0.2b | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 14.0ms | (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 9.0ms | (/ (* 4 (* c a)) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 7.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
4 calls. Slowest were:
| 168.0ms | (/ (* 4 (* c a)) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 96.0ms | (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 67.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 42.0ms | (- (* b b) (* (* 4 a) c)) |
143 calls. Slowest were:
| 1.1s | (/ (* (* (* 4 (* c a)) (* 4 (* c a))) (* 4 (* c a))) (* (* (* (* 2 a) (* 2 a)) (* 2 a)) (* (* (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))))) |
| 758.0ms | (/ (* (* (* 4 (* c a)) (* 4 (* c a))) (* 4 (* c a))) (* (* (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))))) |
| 555.0ms | (* (* 2 a) (fma (- (sqrt (sqrt (- (* b b) (* (* 4 a) c))))) (sqrt (sqrt (- (* b b) (* (* 4 a) c)))) (* (sqrt (sqrt (- (* b b) (* (* 4 a) c)))) (sqrt (sqrt (- (* b b) (* (* 4 a) c))))))) |
10 alts after pruning (10 fresh and 0 done)
Merged error: 0.0b
Found 4 expressions with local error:
| 0.2b | (/ (/ (* 4 (* c a)) (* 2 a)) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
| 0.2b | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 0.2b | (fma (* a -4) c (* b b)) |
| 0.1b | (/ (* 4 (* c a)) (* 2 a)) |
4 calls. Slowest were:
| 9.0ms | (/ (/ (* 4 (* c a)) (* 2 a)) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
| 6.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 5.0ms | (/ (* 4 (* c a)) (* 2 a)) |
4 calls. Slowest were:
| 122.0ms | (/ (/ (* 4 (* c a)) (* 2 a)) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
| 42.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 38.0ms | (fma (* a -4) c (* b b)) |
| 7.0ms | (/ (* 4 (* c a)) (* 2 a)) |
149 calls. Slowest were:
| 1.1s | (sqrt (/ (* 4 (* c a)) (* 2 a))) |
| 1.1s | (/ (sqrt (/ (* 4 (* c a)) (* 2 a))) (sqrt (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 1.1s | (/ (/ (* (* (* 4 (* c a)) (* 4 (* c a))) (* 4 (* c a))) (* (* (* 2 a) (* 2 a)) (* 2 a))) (* (* (- (- b) (sqrt (fma (* a -4) c (* b b)))) (- (- b) (sqrt (fma (* a -4) c (* b b))))) (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 0b
Found 4 expressions with local error:
| 0.3b | (/ 2 (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (/ (* c a) a))) |
| 0.2b | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 0.2b | (fma (* a -4) c (* b b)) |
| 0.2b | (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (/ (* c a) a)) |
4 calls. Slowest were:
| 20.0ms | (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (/ (* c a) a)) |
| 9.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 5.0ms | (/ 2 (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (/ (* c a) a))) |
4 calls. Slowest were:
| 69.0ms | (/ 2 (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (/ (* c a) a))) |
| 51.0ms | (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (/ (* c a) a)) |
| 50.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 31.0ms | (fma (* a -4) c (* b b)) |
158 calls. Slowest were:
| 879.0ms | (- (* (- b) (- b)) (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b))))) |
| 732.0ms | (/ (- (sqrt (- b)) (sqrt (sqrt (fma (* a -4) c (* b b))))) (/ (* c a) a)) |
| 545.0ms | (sqrt (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (/ (* c a) a))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 0b
0% (0.4b remaining)
Error of 0.4b against oracle of 0.0b and baseline of 0.4b