



Bits error versus a




Bits error versus b




Bits error versus c
| Original | 33.3 |
|---|---|
| Target | 20.4 |
| Herbie | 9.7 |
if b < -1.7715196556134555e+46Initial program 56.2
Taylor expanded around inf 56.2
Simplified56.2
rmApplied add-sqr-sqrt56.2
Applied sqrt-prod58.3
Applied add-sqr-sqrt56.2
Applied prod-diff56.3
Simplified60.6
Simplified56.2
Taylor expanded around -inf 4.3
Simplified4.3
if -1.7715196556134555e+46 < b < -1.1090157738250754e-76Initial program 42.2
Taylor expanded around inf 42.2
Simplified42.2
rmApplied flip--42.3
Applied associate-/l/45.0
Simplified18.1
if -1.1090157738250754e-76 < b < 5.973421650485517e+105Initial program 13.0
Taylor expanded around inf 13.0
Simplified13.0
rmApplied add-sqr-sqrt13.0
Applied sqrt-prod13.2
Applied add-sqr-sqrt48.6
Applied prod-diff48.6
Simplified13.1
Simplified13.0
if 5.973421650485517e+105 < b Initial program 46.5
Taylor expanded around inf 46.5
Simplified46.4
rmApplied add-sqr-sqrt46.4
Applied sqrt-prod46.5
Applied add-sqr-sqrt62.3
Applied prod-diff62.3
Simplified47.6
Simplified46.5
Taylor expanded around inf 3.8
Final simplification9.7
herbie shell --seed 2019002 +o rules:numerics
(FPCore (a b c)
:name "The quadratic formula (r2)"
:herbie-target
(if (< b 0) (/ c (* 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))))) (* 2 a)))
Time bar (total: 29.9s)Debug log
1 calls. Slowest were:
| 76.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
2 alts after pruning (2 fresh and 0 done)
Merged error: 32.3b
Found 4 expressions with local error:
| 22.2b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 1.9b | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 0.6b | (- (* b b) (* 4 (* a c))) |
| 0.1b | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
4 calls. Slowest were:
| 36.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 14.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 6.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
4 calls. Slowest were:
| 79.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 46.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 41.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 31.0ms | (- (* b b) (* 4 (* a c))) |
98 calls. Slowest were:
| 405.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c)))))) |
| 382.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
| 346.0ms | (- (pow (- b) 3) (pow (sqrt (- (* b b) (* 4 (* a c)))) 3)) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 5.0b
Found 3 expressions with local error:
| 22.2b | (sqrt (fma (* a -4) c (* b b))) |
| 1.9b | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 0.1b | (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (* 2 a)) |
3 calls. Slowest were:
| 24.0ms | (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (* 2 a)) |
| 10.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 2.0ms | (sqrt (fma (* a -4) c (* b b))) |
3 calls. Slowest were:
| 70.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 69.0ms | (/ (- (- b) (sqrt (fma (* a -4) c (* b b)))) (* 2 a)) |
| 68.0ms | (sqrt (fma (* a -4) c (* b b))) |
90 calls. Slowest were:
| 256.0ms | (- (* (- b) (- b)) (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b))))) |
| 221.0ms | (- (pow (- b) 3) (pow (sqrt (fma (* a -4) c (* b b))) 3)) |
| 106.0ms | (/ (* (* (- (- b) (sqrt (fma (* a -4) c (* b b)))) (- (- b) (sqrt (fma (* a -4) c (* b b))))) (- (- b) (sqrt (fma (* a -4) c (* b b))))) (* (* (* 2 a) (* 2 a)) (* 2 a))) |
6 alts after pruning (6 fresh and 0 done)
Merged error: 5.0b
Found 4 expressions with local error:
| 22.2b | (sqrt (fma -4 (* a c) (* b b))) |
| 1.9b | (- (- b) (sqrt (fma -4 (* a c) (* b b)))) |
| 0.6b | (fma -4 (* a c) (* b b)) |
| 0.1b | (/ (+ (- (- b) (sqrt (fma -4 (* a c) (* b b)))) 0) (* 2 a)) |
4 calls. Slowest were:
| 87.0ms | (/ (+ (- (- b) (sqrt (fma -4 (* a c) (* b b)))) 0) (* 2 a)) |
| 5.0ms | (- (- b) (sqrt (fma -4 (* a c) (* b b)))) |
| 1.0ms | (sqrt (fma -4 (* a c) (* b b))) |
4 calls. Slowest were:
| 62.0ms | (/ (+ (- (- b) (sqrt (fma -4 (* a c) (* b b)))) 0) (* 2 a)) |
| 62.0ms | (- (- b) (sqrt (fma -4 (* a c) (* b b)))) |
| 44.0ms | (sqrt (fma -4 (* a c) (* b b))) |
| 30.0ms | (fma -4 (* a c) (* b b)) |
85 calls. Slowest were:
| 148.0ms | (- (* (- (- b) (sqrt (fma -4 (* a c) (* b b)))) (- (- b) (sqrt (fma -4 (* a c) (* b b))))) (* 0 0)) |
| 127.0ms | (/ (+ (- (- b) (sqrt (fma -4 (* a c) (* b b)))) 0) a) |
| 125.0ms | (- (log (+ (- (- b) (sqrt (fma -4 (* a c) (* b b)))) 0)) (log (* 2 a))) |
6 alts after pruning (5 fresh and 1 done)
Merged error: 5.0b
Found 4 expressions with local error:
| 22.2b | (sqrt (fma (* a -4) c (* b b))) |
| 1.9b | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 0.2b | (/ 1 (/ (* 2 a) (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 0.1b | (/ (* 2 a) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
4 calls. Slowest were:
| 8.0ms | (/ (* 2 a) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
| 5.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 4.0ms | (/ 1 (/ (* 2 a) (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
4 calls. Slowest were:
| 110.0ms | (/ (* 2 a) (- (- b) (sqrt (fma (* a -4) c (* b b))))) |
| 81.0ms | (/ 1 (/ (* 2 a) (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 72.0ms | (- (- b) (sqrt (fma (* a -4) c (* b b)))) |
| 46.0ms | (sqrt (fma (* a -4) c (* b b))) |
108 calls. Slowest were:
| 532.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b)))) (* (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 406.0ms | (cbrt (/ (* 2 a) (- (- b) (sqrt (fma (* a -4) c (* b b)))))) |
| 404.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (fma (* a -4) c (* b b))) (sqrt (fma (* a -4) c (* b b)))) (* (- b) (sqrt (fma (* a -4) c (* b b)))))) |
6 alts after pruning (4 fresh and 2 done)
Merged error: 5.0b
85.3% (4.1b remaining)
Error of 9.7b against oracle of 5.7b and baseline of 33.3b