


Bits error versus a



Bits error versus b_2



Bits error versus c
Results
if b_2 < -3.434444833633684e+56Initial program 37.6
Simplified37.6
Taylor expanded around inf 37.6
Taylor expanded around -inf 5.1
if -3.434444833633684e+56 < b_2 < 5.633475377269129e-146Initial program 11.2
Simplified11.2
Taylor expanded around inf 11.2
if 5.633475377269129e-146 < b_2 Initial program 49.7
Simplified49.7
Taylor expanded around inf 49.7
Taylor expanded around inf 12.5
Final simplification10.6
herbie shell --seed 2019021
(FPCore (a b_2 c)
:name "quad2p (problem 3.2.1, positive)"
(/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a))
Time bar (total: 14.9s)Debug log
| 43.0ms | 76× | body | 80 | nan |
| 20.0ms | 37× | body | 2560 | valid |
| 16.0ms | 189× | body | 80 | valid |
| 9.0ms | 22× | body | 1280 | valid |
| 6.0ms | 18× | body | 640 | valid |
| 3.0ms | 11× | body | 320 | valid |
| 2.0ms | 4× | body | 5120 | valid |
| 0.0ms | 3× | body | 160 | valid |
| 437.0ms | (/ (+ (- b_2) (sqrt (- (* b_2 b_2) (* a c)))) a) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 33.3b
Found 4 expressions with local error:
| 22.9b | (sqrt (- (* b_2 b_2) (* a c))) |
| 1.1b | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 0.6b | (- (* b_2 b_2) (* a c)) |
| 0.1b | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 7× | add-sqr-sqrt |
| 6× | add-log-exp |
| 6× | add-cube-cbrt |
| 6× | *-un-lft-identity |
| 4× | flip-- |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | flip3-- |
| 4× | pow1 |
| 3× | associate-/l* |
| 3× | sqrt-prod |
| 2× | associate-/l/ |
| 2× | sqrt-div |
| 2× | sub-neg |
| 1× | div-inv |
| 1× | div-sub |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 20.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 15.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 6.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 2.0ms | (- (* b_2 b_2) (* a c)) |
| 82.0ms | (/ (- (sqrt (- (* b_2 b_2) (* a c))) b_2) a) |
| 65.0ms | (- (sqrt (- (* b_2 b_2) (* a c))) b_2) |
| 56.0ms | (sqrt (- (* b_2 b_2) (* a c))) |
| 33.0ms | (- (* b_2 b_2) (* a c)) |
| 328.0ms | (sqrt (- (pow (* b_2 b_2) 3) (pow (* a c) 3))) |
| 300.0ms | (- (* (sqrt (- (* b_2 b_2) (* a c))) (sqrt (- (* b_2 b_2) (* a c)))) (* b_2 b_2)) |
| 225.0ms | (sqrt (- (* (* b_2 b_2) (* b_2 b_2)) (* (* a c) (* a c)))) |
| 132.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
| 73.0ms | (- (pow (sqrt (- (* b_2 b_2) (* a c))) 3) (pow b_2 3)) |
11 alts after pruning (11 fresh and 0 done)
Merged error: 5.7b
Found 4 expressions with local error:
| 22.9b | (sqrt (- (pow b_2 2) (* a c))) |
| 1.1b | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 0.6b | (- (pow b_2 2) (* a c)) |
| 0.1b | (/ (- (sqrt (- (pow b_2 2) (* a c))) b_2) a) |
| 7× | add-sqr-sqrt |
| 6× | add-log-exp |
| 6× | add-cube-cbrt |
| 6× | *-un-lft-identity |
| 4× | flip-- |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | flip3-- |
| 4× | pow1 |
| 3× | associate-/l* |
| 3× | sqrt-prod |
| 2× | associate-/l/ |
| 2× | sqrt-div |
| 2× | sub-neg |
| 1× | div-inv |
| 1× | div-sub |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 54.0ms | (/ (- (sqrt (- (pow b_2 2) (* a c))) b_2) a) |
| 14.0ms | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 8.0ms | (- (pow b_2 2) (* a c)) |
| 3.0ms | (sqrt (- (pow b_2 2) (* a c))) |
| 77.0ms | (/ (- (sqrt (- (pow b_2 2) (* a c))) b_2) a) |
| 54.0ms | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 51.0ms | (sqrt (- (pow b_2 2) (* a c))) |
| 35.0ms | (- (pow b_2 2) (* a c)) |
| 365.0ms | (- (* (sqrt (- (pow b_2 2) (* a c))) (sqrt (- (pow b_2 2) (* a c)))) (* b_2 b_2)) |
| 312.0ms | (sqrt (- (pow (pow b_2 2) 3) (pow (* a c) 3))) |
| 267.0ms | (sqrt (- (* (pow b_2 2) (pow b_2 2)) (* (* a c) (* a c)))) |
| 144.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
| 91.0ms | (- (pow (sqrt (- (pow b_2 2) (* a c))) 3) (pow b_2 3)) |
11 alts after pruning (10 fresh and 1 done)
Merged error: 5.7b
Found 4 expressions with local error:
| 22.9b | (sqrt (- (pow b_2 2) (* a c))) |
| 1.1b | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 0.6b | (- (pow b_2 2) (* a c)) |
| 0.2b | (/ 1 (/ a (- (sqrt (- (pow b_2 2) (* a c))) b_2))) |
| 7× | add-sqr-sqrt |
| 6× | add-log-exp |
| 6× | add-cube-cbrt |
| 6× | *-un-lft-identity |
| 6× | associate-/r* |
| 5× | add-exp-log |
| 5× | pow1 |
| 4× | flip-- |
| 4× | add-cbrt-cube |
| 4× | flip3-- |
| 3× | sqrt-prod |
| 3× | associate-/r/ |
| 2× | div-inv |
| 2× | sqrt-div |
| 2× | sub-neg |
| 1× | inv-pow |
| 1× | pow-flip |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | rec-exp |
| 1× | rem-sqrt-square |
| 23.0ms | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 6.0ms | (sqrt (- (pow b_2 2) (* a c))) |
| 3.0ms | (- (pow b_2 2) (* a c)) |
| 2.0ms | (/ 1 (/ a (- (sqrt (- (pow b_2 2) (* a c))) b_2))) |
| 88.0ms | (/ 1 (/ a (- (sqrt (- (pow b_2 2) (* a c))) b_2))) |
| 65.0ms | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 50.0ms | (sqrt (- (pow b_2 2) (* a c))) |
| 32.0ms | (- (pow b_2 2) (* a c)) |
| 302.0ms | (sqrt (- (pow (pow b_2 2) 3) (pow (* a c) 3))) |
| 238.0ms | (sqrt (- (* (pow b_2 2) (pow b_2 2)) (* (* a c) (* a c)))) |
| 145.0ms | (+ (* (sqrt (- (pow b_2 2) (* a c))) (sqrt (- (pow b_2 2) (* a c)))) (+ (* b_2 b_2) (* (sqrt (- (pow b_2 2) (* a c))) b_2))) |
| 107.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
| 70.0ms | (- (* 1/2 (/ (* a c) b_2)) b_2) |
10 alts after pruning (9 fresh and 1 done)
Merged error: 5.7b
Found 4 expressions with local error:
| 22.9b | (sqrt (- (pow b_2 2) (* a c))) |
| 22.9b | (sqrt (- (pow b_2 2) (* a c))) |
| 1.1b | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 1.1b | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 8× | add-sqr-sqrt |
| 6× | sqrt-prod |
| 6× | add-cube-cbrt |
| 6× | *-un-lft-identity |
| 4× | add-log-exp |
| 4× | flip-- |
| 4× | add-exp-log |
| 4× | sqrt-div |
| 4× | add-cbrt-cube |
| 4× | flip3-- |
| 4× | pow1 |
| 2× | pow1/2 |
| 2× | sub-neg |
| 2× | rem-sqrt-square |
| 26.0ms | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 15.0ms | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 5.0ms | (sqrt (- (pow b_2 2) (* a c))) |
| 3.0ms | (sqrt (- (pow b_2 2) (* a c))) |
| 61.0ms | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 48.0ms | (- (sqrt (- (pow b_2 2) (* a c))) b_2) |
| 43.0ms | (sqrt (- (pow b_2 2) (* a c))) |
| 41.0ms | (sqrt (- (pow b_2 2) (* a c))) |
| 346.0ms | (sqrt (- (pow (pow b_2 2) 3) (pow (* a c) 3))) |
| 280.0ms | (sqrt (- (* (pow b_2 2) (pow b_2 2)) (* (* a c) (* a c)))) |
| 271.0ms | (sqrt (- (pow (pow b_2 2) 3) (pow (* a c) 3))) |
| 204.0ms | (sqrt (- (* (pow b_2 2) (pow b_2 2)) (* (* a c) (* a c)))) |
| 125.0ms | (- b_2 (* 1/2 (/ (* a c) b_2))) |
10 alts after pruning (8 fresh and 2 done)
Merged error: 5.7b
81.3% (5.2b remaining)
Error of 10.6b against oracle of 5.4b and baseline of 33.4b
| 1.1s | 5833× | body | 80 | valid |
| 826.0ms | 1092× | body | 2560 | valid |
| 682.0ms | 912× | body | 1280 | valid |
| 377.0ms | 2857× | body | 80 | nan |
| 204.0ms | 471× | body | 640 | valid |
| 116.0ms | 112× | body | 5120 | valid |
| 105.0ms | 224× | body | 320 | valid |
| 23.0ms | 100× | body | 160 | valid |