


Bits error versus a



Bits error versus b



Bits error versus c
Results
if b < -4.21731702933355e+148Initial program 58.9
Simplified58.9
Taylor expanded around -inf 2.3
if -4.21731702933355e+148 < b < 1.2612246360442721e-22Initial program 14.8
Simplified14.8
if 1.2612246360442721e-22 < b Initial program 54.5
Simplified54.5
Taylor expanded around inf 6.4
Simplified6.4
Final simplification10.3
herbie shell --seed 2019007
(FPCore (a b c)
:name "Quadratic roots, full range"
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
Time bar (total: 29.2s)Debug log
1 calls. Slowest were:
| 64.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 34.1b
Found 4 expressions with local error:
| 24.1b | (sqrt (- (* b b) (* (* c a) 4))) |
| 1.5b | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 1.2b | (- (* b b) (* (* c a) 4)) |
| 0.1b | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
4 calls. Slowest were:
| 20.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 10.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 2.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 61.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 47.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 44.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 34.0ms | (- (* b b) (* (* c a) 4)) |
32 calls. Slowest were:
| 346.0ms | (- (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (* b b)) |
| 319.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 238.0ms | (- (pow (sqrt (- (* b b) (* (* c a) 4))) 3) (pow b 3)) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 6.1b
Found 4 expressions with local error:
| 24.1b | (sqrt (- (* b b) (* (* c a) 4))) |
| 1.5b | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 1.2b | (- (* b b) (* (* c a) 4)) |
| 0.1b | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
4 calls. Slowest were:
| 12.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 12.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 3.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 124.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 55.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 47.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 26.0ms | (- (* b b) (* (* c a) 4)) |
30 calls. Slowest were:
| 362.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 198.0ms | (+ (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (+ (* b b) (* (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 165.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
10 alts after pruning (8 fresh and 2 done)
Merged error: 6.1b
Found 4 expressions with local error:
| 24.1b | (sqrt (- (* b b) (* (* c a) 4))) |
| 24.1b | (sqrt (- (* b b) (* (* c a) 4))) |
| 1.5b | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 1.5b | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
4 calls. Slowest were:
| 11.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 10.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 2.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 66.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 62.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 45.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 43.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
24 calls. Slowest were:
| 386.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 375.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 162.0ms | (- b (* 2 (/ (* a c) b))) |
10 alts after pruning (7 fresh and 3 done)
Merged error: 6.1b
Found 4 expressions with local error:
| 24.9b | (cbrt (- (* b b) (* (* c a) 4))) |
| 24.9b | (cbrt (- (* b b) (* (* c a) 4))) |
| 24.9b | (cbrt (- (* b b) (* (* c a) 4))) |
| 24.1b | (sqrt (* (* (cbrt (- (* b b) (* (* c a) 4))) (cbrt (- (* b b) (* (* c a) 4)))) (cbrt (- (* b b) (* (* c a) 4))))) |
4 calls. Slowest were:
| 18.0ms | (sqrt (* (* (cbrt (- (* b b) (* (* c a) 4))) (cbrt (- (* b b) (* (* c a) 4)))) (cbrt (- (* b b) (* (* c a) 4))))) |
| 3.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
| 3.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 110.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
| 101.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
| 93.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
| 52.0ms | (sqrt (* (* (cbrt (- (* b b) (* (* c a) 4))) (cbrt (- (* b b) (* (* c a) 4)))) (cbrt (- (* b b) (* (* c a) 4))))) |
54 calls. Slowest were:
| 795.0ms | (sqrt (* (* (cbrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) (cbrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3)))) (cbrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))))) |
| 593.0ms | (sqrt (* (* (cbrt (- (* b b) (* (* c a) 4))) (cbrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3)))) (cbrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))))) |
| 577.0ms | (sqrt (* (* (cbrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) (cbrt (- (* b b) (* (* c a) 4)))) (cbrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))))) |
10 alts after pruning (6 fresh and 4 done)
Merged error: 6.1b
83.3% (4.8b remaining)
Error of 10.3b against oracle of 5.5b and baseline of 34.1b