



Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 34.2 |
|---|---|
| Target | 21.5 |
| Herbie | 10.3 |
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 "The quadratic formula (r1)"
:herbie-target
(if (< b 0) (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) (/ c (* a (/ (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))))
(/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)))
Time bar (total: 30.1s)Debug log
1 calls. Slowest were:
| 55.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:
| 41.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 19.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 4.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 62.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 44.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 40.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 26.0ms | (- (* b b) (* (* c a) 4)) |
32 calls. Slowest were:
| 417.0ms | (- (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (* b b)) |
| 303.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 239.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:
| 16.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 12.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 4.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 106.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 66.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 56.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 42.0ms | (- (* b b) (* (* c a) 4)) |
30 calls. Slowest were:
| 308.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 195.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
| 146.0ms | (+ (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (+ (* b b) (* (sqrt (- (* b b) (* (* c a) 4))) b))) |
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:
| 24.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 21.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 4.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 62.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 60.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 48.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 47.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
24 calls. Slowest were:
| 376.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 273.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 144.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
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:
| 13.0ms | (sqrt (* (* (cbrt (- (* b b) (* (* c a) 4))) (cbrt (- (* b b) (* (* c a) 4)))) (cbrt (- (* b b) (* (* c a) 4))))) |
| 1.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
| 1.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 104.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
| 103.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
| 96.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
| 68.0ms | (sqrt (* (* (cbrt (- (* b b) (* (* c a) 4))) (cbrt (- (* b b) (* (* c a) 4)))) (cbrt (- (* b b) (* (* c a) 4))))) |
54 calls. Slowest were:
| 645.0ms | (sqrt (* (* (cbrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) (cbrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3)))) (cbrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))))) |
| 607.0ms | (sqrt (* (* (cbrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) (cbrt (- (* b b) (* (* c a) 4)))) (cbrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))))) |
| 600.0ms | (sqrt (* (* (cbrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) (cbrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3)))) (cbrt (- (* b b) (* (* 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