


Bits error versus a



Bits error versus b



Bits error versus c
Results
if b < -6.341219225938887e+106Initial program 45.4
rmApplied div-inv45.4
Taylor expanded around -inf 3.1
if -6.341219225938887e+106 < b < 3.45159933628946e-306Initial program 8.9
rmApplied div-inv9.0
rmApplied un-div-inv8.9
if 3.45159933628946e-306 < b < 2.408724732699167e+98Initial program 32.5
rmApplied div-inv32.5
rmApplied flip-+32.6
Applied associate-*l/32.7
Simplified16.4
Taylor expanded around 0 9.1
rmApplied *-un-lft-identity9.1
Applied associate-/r*9.1
if 2.408724732699167e+98 < b Initial program 58.7
rmApplied div-inv58.7
Taylor expanded around inf 2.9
Simplified2.9
Final simplification6.7
herbie shell --seed 2019008
(FPCore (a b c)
:name "Quadratic roots, full range"
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
Time bar (total: 21.2s)Debug log
1 calls. Slowest were:
| 75.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 34.5b
Found 4 expressions with local error:
| 23.3b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 1.0b | (- (* b b) (* (* 4 a) c)) |
| 0.9b | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 0.1b | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
4 calls. Slowest were:
| 44.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 9.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 2.0ms | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 93.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 72.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 61.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 42.0ms | (- (* b b) (* (* 4 a) c)) |
36 calls. Slowest were:
| 417.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 260.0ms | (+ (pow (- b) 3) (pow (sqrt (- (* b b) (* (* 4 a) c))) 3)) |
| 233.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c))))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 6.9b
Found 4 expressions with local error:
| 23.3b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 1.0b | (- (* b b) (* (* 4 a) c)) |
| 0.9b | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 0.2b | (* (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (/ 1 (* 2 a))) |
4 calls. Slowest were:
| 15.0ms | (* (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (/ 1 (* 2 a))) |
| 11.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 4.0ms | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 99.0ms | (* (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (/ 1 (* 2 a))) |
| 59.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 52.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 45.0ms | (- (* b b) (* (* 4 a) c)) |
39 calls. Slowest were:
| 384.0ms | (* (- (* (- b) (- b)) (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c))))) (/ 1 (* 2 a))) |
| 374.0ms | (* (- (* (- b) (- b)) (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c))))) 1) |
| 363.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
9 alts after pruning (9 fresh and 0 done)
Merged error: 6.9b
Found 4 expressions with local error:
| 23.3b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 15.2b | (* (* (* c a) 4) (/ 1/2 a)) |
| 1.0b | (- (* b b) (* (* 4 a) c)) |
| 0.8b | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
4 calls. Slowest were:
| 13.0ms | (* (* (* c a) 4) (/ 1/2 a)) |
| 11.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 4.0ms | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 76.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 60.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 44.0ms | (- (* b b) (* (* 4 a) c)) |
| 8.0ms | (* (* (* c a) 4) (/ 1/2 a)) |
33 calls. Slowest were:
| 425.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 386.0ms | (+ (log (* (* c a) 4)) (log (/ 1/2 a))) |
| 352.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 6.9b
Found 4 expressions with local error:
| 23.3b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 1.0b | (- (* b b) (* (* 4 a) c)) |
| 0.8b | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 0.1b | (/ (* 2 c) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) |
4 calls. Slowest were:
| 10.0ms | (/ (* 2 c) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 8.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 2.0ms | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 165.0ms | (/ (* 2 c) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 77.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 50.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 34.0ms | (- (* b b) (* (* 4 a) c)) |
47 calls. Slowest were:
| 424.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c)))) (* (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 399.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 164.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 6.9b
95.3% (1.1b remaining)
Error of 6.7b against oracle of 5.6b and baseline of 28.9b