



Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 34.1 |
|---|---|
| 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 "quadp (p42, positive)"
: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: 28.3s)Debug log
1 calls. Slowest were:
| 91.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) (* (* 4 a) c))) |
| 1.5b | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 1.2b | (- (* b b) (* (* 4 a) c)) |
| 0.1b | (/ (- (sqrt (- (* b b) (* (* 4 a) c))) b) (* 2 a)) |
4 calls. Slowest were:
| 40.0ms | (/ (- (sqrt (- (* b b) (* (* 4 a) c))) b) (* 2 a)) |
| 21.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 4.0ms | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 97.0ms | (/ (- (sqrt (- (* b b) (* (* 4 a) c))) b) (* 2 a)) |
| 70.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 43.0ms | (- (* b b) (* (* 4 a) c)) |
| 43.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
32 calls. Slowest were:
| 378.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 352.0ms | (- (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c)))) (* b b)) |
| 250.0ms | (- (pow (sqrt (- (* b b) (* (* 4 a) c))) 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) (* (* 4 a) c))) |
| 1.5b | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 1.2b | (- (* b b) (* (* 4 a) c)) |
| 0.1b | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* 4 a) c))) b))) |
4 calls. Slowest were:
| 20.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 7.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 4.0ms | (- (* b b) (* (* 4 a) c)) |
4 calls. Slowest were:
| 94.0ms | (/ 1 (/ (* 2 a) (- (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 76.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 55.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 47.0ms | (- (* b b) (* (* 4 a) c)) |
30 calls. Slowest were:
| 380.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 228.0ms | (+ (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c)))) (+ (* b b) (* (sqrt (- (* b b) (* (* 4 a) c))) b))) |
| 162.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
10 alts after pruning (8 fresh and 2 done)
Merged error: 6.1b
Found 4 expressions with local error:
| 24.1b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 24.1b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 1.5b | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 1.5b | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
4 calls. Slowest were:
| 11.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 11.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 2.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
4 calls. Slowest were:
| 62.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
| 49.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 48.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 46.0ms | (- (sqrt (- (* b b) (* (* 4 a) c))) b) |
24 calls. Slowest were:
| 386.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 366.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 169.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) (* (* 4 a) c))) |
| 24.9b | (cbrt (- (* b b) (* (* 4 a) c))) |
| 24.9b | (cbrt (- (* b b) (* (* 4 a) c))) |
| 24.1b | (sqrt (* (* (cbrt (- (* b b) (* (* 4 a) c))) (cbrt (- (* b b) (* (* 4 a) c)))) (cbrt (- (* b b) (* (* 4 a) c))))) |
4 calls. Slowest were:
| 8.0ms | (sqrt (* (* (cbrt (- (* b b) (* (* 4 a) c))) (cbrt (- (* b b) (* (* 4 a) c)))) (cbrt (- (* b b) (* (* 4 a) c))))) |
| 2.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 1.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
4 calls. Slowest were:
| 122.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 108.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 103.0ms | (cbrt (- (* b b) (* (* 4 a) c))) |
| 50.0ms | (sqrt (* (* (cbrt (- (* b b) (* (* 4 a) c))) (cbrt (- (* b b) (* (* 4 a) c)))) (cbrt (- (* b b) (* (* 4 a) c))))) |
54 calls. Slowest were:
| 575.0ms | (sqrt (* (* (cbrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) (cbrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3)))) (cbrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))))) |
| 527.0ms | (sqrt (* (* (cbrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) (cbrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c))))) (cbrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))))) |
| 523.0ms | (sqrt (* (* (cbrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) (cbrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c))))) (cbrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 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.2b