



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: 19.6s)Debug log
1 calls. Slowest were:
| 105.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:
| 27.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 4.0ms | (* (* c a) 4) |
| 2.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 67.0ms | (/ (- (sqrt (- (* b b) (* (* c a) 4))) b) (* 2 a)) |
| 54.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 25.0ms | (- (* b b) (* (* c a) 4)) |
| 16.0ms | (* (* c a) 4) |
32 calls. Slowest were:
| 382.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 242.0ms | (- (* (sqrt (- (* b b) (* (* c a) 4))) (sqrt (- (* b b) (* (* c a) 4)))) (* b b)) |
| 187.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:
| 11.0ms | (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 6.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 6.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 130.0ms | (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b)) |
| 82.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* c a) 4))) b))) |
| 46.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 33.0ms | (- (* b b) (* (* c a) 4)) |
43 calls. Slowest were:
| 334.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 180.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* c a) 4) (* (* c a) 4)))) |
| 169.0ms | (- b (* 2 (/ (* a c) 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:
| 133.0ms | (cbrt (* (* c a) 4)) |
| 119.0ms | (sqrt (- (* b b) (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4))))) |
| 115.0ms | (cbrt (* (* c a) 4)) |
| 97.0ms | (cbrt (* (* c a) 4)) |
21 calls. Slowest were:
| 505.0ms | (* (exp (* -1/3 (+ (log (/ 1 a)) (log (/ 1 c))))) (cbrt 4)) |
| 467.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* (cbrt (* (* c a) 4)) (cbrt (* (* c a) 4))) (cbrt (* (* c a) 4))) 3))) |
| 229.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:
| 2.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 2.0ms | (- (* b b) (* (* c a) 4)) |
| 2.0ms | (- (* b b) (* (* c a) 4)) |
4 calls. Slowest were:
| 65.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 56.0ms | (sqrt (- (* b b) (* (* c a) 4))) |
| 27.0ms | (- (* b b) (* (* c a) 4)) |
| 26.0ms | (- (* b b) (* (* c a) 4)) |
26 calls. Slowest were:
| 377.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 330.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* c a) 4) 3))) |
| 169.0ms | (- b (* 2 (/ (* a c) b))) |
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