



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: 24.1s)Debug log
1 calls. Slowest were:
| 71.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:
| 39.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 25.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 4.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 61.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 48.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 45.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 41.0ms | (- (* b b) (* (* c a) 4)) |
32 calls. Slowest were:
| 303.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 279.0ms | (- (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (* b b)) |
| 233.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:
| 11.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 7.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 2.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 105.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 64.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 44.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 38.0ms | (- (* b b) (* (* c a) 4)) |
30 calls. Slowest were:
| 284.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 172.0ms | (+ (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (+ (* b b) (* (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 130.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:
| 20.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 11.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 2.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 67.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 53.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
| 47.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 46.0ms | (- (sqrt (- (* b b) (* (* c a) 4))) b) |
24 calls. Slowest were:
| 350.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 275.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 127.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:
| 8.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:
| 96.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
| 92.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
| 77.0ms | (cbrt (- (* b b) (* (* c a) 4))) |
| 45.0ms | (sqrt (* (* (cbrt (- (* b b) (* (* c a) 4))) (cbrt (- (* b b) (* (* c a) 4)))) (cbrt (- (* b b) (* (* c a) 4))))) |
54 calls. Slowest were:
| 577.0ms | (sqrt (* (* (cbrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) (cbrt (- (* b b) (* (* c a) 4)))) (cbrt (- (* b b) (* (* c a) 4))))) |
| 501.0ms | (sqrt (* (* (cbrt (- (* b b) (* (* c a) 4))) (cbrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4))))) (cbrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))))) |
| 491.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))))) |
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