



Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.5 |
|---|---|
| Target | 20.5 |
| Herbie | 9.7 |
if b < -2.835340391058544e+148Initial program 58.7
Simplified58.7
Taylor expanded around -inf 2.1
if -2.835340391058544e+148 < b < 3.587545867948427e-80Initial program 12.5
Simplified12.5
if 3.587545867948427e-80 < b Initial program 52.6
Simplified52.6
Taylor expanded around inf 8.5
Simplified8.5
Final simplification9.7
herbie shell --seed 2019004
(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: 15.4s)Debug log
1 calls. Slowest were:
| 97.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 32.6b
Found 4 expressions with local error:
| 22.0b | (sqrt (- (* b b) (* (* c a) 4))) |
| 0.8b | (- (* b b) (* (* c a) 4)) |
| 0.3b | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 0.0b | (* (* c a) 4) |
4 calls. Slowest were:
| 39.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 2.0ms | (* (* c a) 4) |
| 2.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 64.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 56.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 37.0ms | (- (* b b) (* (* c a) 4)) |
| 17.0ms | (* (* c a) 4) |
32 calls. Slowest were:
| 307.0ms | (- (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (* b b)) |
| 254.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 143.0ms | (- (pow (sqrt (- (* b b) (* (* c a) 4))) 3) (pow b 3)) |
8 alts after pruning (7 fresh and 1 done)
Merged error: 4.8b
Found 4 expressions with local error:
| 22.0b | (sqrt (- (* b b) (* (* c a) 4))) |
| 0.8b | (- (* b b) (* (* c a) 4)) |
| 0.4b | (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 0.2b | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
4 calls. Slowest were:
| 16.0ms | (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 7.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 2.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 145.0ms | (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 86.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 49.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 26.0ms | (- (* b b) (* (* c a) 4)) |
43 calls. Slowest were:
| 282.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 163.0ms | (+ (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (+ (* b b) (* (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 149.0ms | (+ (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (+ (* b b) (* (sqrt (- (* b b) (* (* c a) 4))) b))) |
9 alts after pruning (7 fresh and 2 done)
Merged error: 4.8b
Found 4 expressions with local error:
| 22.0b | (sqrt (- (* b b) (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4))))) |
| 16.3b | (cbrt (* (* c a) 4)) |
| 16.3b | (cbrt (* (* c a) 4)) |
| 16.3b | (cbrt (* (* c a) 4)) |
4 calls. Slowest were:
| 2.0ms | (sqrt (- (* b b) (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4))))) |
| 1.0ms | (cbrt (* (* c a) 4)) |
| 1.0ms | (cbrt (* (* c a) 4)) |
4 calls. Slowest were:
| 122.0ms | (cbrt (* (* c a) 4)) |
| 114.0ms | (cbrt (* (* c a) 4)) |
| 111.0ms | (sqrt (- (* b b) (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4))))) |
| 109.0ms | (cbrt (* (* c a) 4)) |
21 calls. Slowest were:
| 428.0ms | (- b (* 2 (/ (* a c) b))) |
| 332.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4))) 3))) |
| 205.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4))) (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4)))))) |
9 alts after pruning (6 fresh and 3 done)
Merged error: 4.8b
Found 4 expressions with local error:
| 22.0b | (sqrt (- (* b b) (* (* c a) 4))) |
| 22.0b | (sqrt (- (* b b) (* (* c a) 4))) |
| 0.8b | (- (* b b) (* (* c a) 4)) |
| 0.8b | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 4.0ms | (- (* b b) (* (* c a) 4)) |
| 4.0ms | (- (* b b) (* (* c a) 4)) |
| 4.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
4 calls. Slowest were:
| 52.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 50.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 26.0ms | (- (* b b) (* (* c a) 4)) |
| 26.0ms | (- (* b b) (* (* c a) 4)) |
26 calls. Slowest were:
| 284.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 275.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 142.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
9 alts after pruning (5 fresh and 4 done)
Merged error: 4.8b
85.1% (4.2b remaining)
Error of 9.7b against oracle of 5.6b and baseline of 33.5b