



Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.5 |
|---|---|
| Target | 20.6 |
| Herbie | 10.4 |
if b < -2.3604286795685533e-29Initial program 54.6
Taylor expanded around -inf 6.8
Simplified6.8
if -2.3604286795685533e-29 < b < 3.328175344273382e+63Initial program 15.3
Taylor expanded around inf 15.3
if 3.328175344273382e+63 < b Initial program 37.3
Taylor expanded around inf 37.3
Taylor expanded around inf 5.0
Final simplification10.4
herbie shell --seed 2019016
(FPCore (a b c)
:name "The quadratic formula (r2)"
:herbie-target
(if (< b 0) (/ c (* a (/ (+ (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))) (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))
(/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)))
Time bar (total: 23.4s)Debug log
| 368× | (pre true 80) |
| 182× | (body real 80) |
| 91× | (body nan 80) |
| 35× | (body real 1280) |
| 32× | (body real 2560) |
| 13× | (body real 640) |
| 6× | (body real 320) |
| 5× | (body real 160) |
| 4× | (body real 5120) |
| 96.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 35.2b
Found 4 expressions with local error:
| 24.0b | (sqrt (- (* b b) (* 4 (* a c)))) |
| 1.2b | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 1.2b | (- (* b b) (* 4 (* a c))) |
| 0.1b | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 20× | add-sqr-sqrt |
| 13× | *-un-lft-identity |
| 8× | add-log-exp |
| 7× | add-cube-cbrt |
| 6× | difference-of-squares |
| 6× | associate-/l* |
| 6× | sqrt-prod |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 6× | times-frac |
| 4× | flip-- |
| 4× | flip3-- |
| 4× | pow1 |
| 3× | distribute-lft-out-- |
| 2× | associate-/l/ |
| 2× | diff-log |
| 2× | sqrt-div |
| 2× | sub-neg |
| 1× | div-inv |
| 1× | div-exp |
| 1× | div-sub |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | associate-/r* |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | cbrt-undiv |
| 1× | neg-sub0 |
| 1× | associate--l- |
| 37.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 13.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 4.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 4.0ms | (- (* b b) (* 4 (* a c))) |
| 86.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 64.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 47.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 44.0ms | (- (* b b) (* 4 (* a c))) |
| 361.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
| 283.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c)))))) |
| 215.0ms | (- (pow (- b) 3) (pow (sqrt (- (* b b) (* 4 (* a c)))) 3)) |
| 162.0ms | (sqrt (- (* (* b b) (* b b)) (* (* 4 (* a c)) (* 4 (* a c))))) |
| 106.0ms | (- b (* 2 (/ (* a c) b))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 4.9b
Found 4 expressions with local error:
| 24.0b | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 1.2b | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 1.2b | (- (pow b 2) (* 4 (* a c))) |
| 0.1b | (/ (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) (* 2 a)) |
| 20× | add-sqr-sqrt |
| 13× | *-un-lft-identity |
| 8× | add-log-exp |
| 7× | add-cube-cbrt |
| 6× | difference-of-squares |
| 6× | associate-/l* |
| 6× | sqrt-prod |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 6× | times-frac |
| 4× | flip-- |
| 4× | flip3-- |
| 4× | pow1 |
| 3× | distribute-lft-out-- |
| 2× | associate-/l/ |
| 2× | diff-log |
| 2× | sqrt-div |
| 2× | sub-neg |
| 1× | div-inv |
| 1× | div-exp |
| 1× | div-sub |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | associate-/r* |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | cbrt-undiv |
| 1× | neg-sub0 |
| 1× | associate--l- |
| 24.0ms | (/ (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) (* 2 a)) |
| 12.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 6.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 4.0ms | (- (pow b 2) (* 4 (* a c))) |
| 68.0ms | (/ (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) (* 2 a)) |
| 67.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 60.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 40.0ms | (- (pow b 2) (* 4 (* a c))) |
| 1.0s | (- (* (- b) (- b)) (* (sqrt (- (pow b 2) (* 4 (* a c)))) (sqrt (- (pow b 2) (* 4 (* a c)))))) |
| 315.0ms | (sqrt (- (pow (pow b 2) 3) (pow (* 4 (* a c)) 3))) |
| 236.0ms | (- (pow (- b) 3) (pow (sqrt (- (pow b 2) (* 4 (* a c)))) 3)) |
| 156.0ms | (sqrt (- (* (pow b 2) (pow b 2)) (* (* 4 (* a c)) (* 4 (* a c))))) |
| 104.0ms | (/ (* (* (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) (- (- b) (sqrt (- (pow b 2) (* 4 (* a c)))))) (- (- b) (sqrt (- (pow b 2) (* 4 (* a c)))))) (* (* (* 2 a) (* 2 a)) (* 2 a))) |
9 alts after pruning (8 fresh and 1 done)
Merged error: 4.9b
Found 4 expressions with local error:
| 24.0b | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 1.2b | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 1.2b | (- (pow b 2) (* 4 (* a c))) |
| 0.2b | (/ 1 (/ (* 2 a) (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))))) |
| 16× | add-sqr-sqrt |
| 12× | associate-/r* |
| 11× | *-un-lft-identity |
| 8× | add-log-exp |
| 7× | add-cube-cbrt |
| 7× | add-exp-log |
| 6× | times-frac |
| 5× | sqrt-prod |
| 5× | pow1 |
| 4× | difference-of-squares |
| 4× | flip-- |
| 4× | add-cbrt-cube |
| 4× | flip3-- |
| 3× | associate-/r/ |
| 2× | distribute-lft-out-- |
| 2× | div-inv |
| 2× | diff-log |
| 2× | sqrt-div |
| 2× | sub-neg |
| 2× | rec-exp |
| 1× | inv-pow |
| 1× | pow-flip |
| 1× | div-exp |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | neg-sub0 |
| 1× | associate--l- |
| 14.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 14.0ms | (/ 1 (/ (* 2 a) (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))))) |
| 7.0ms | (- (pow b 2) (* 4 (* a c))) |
| 6.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 92.0ms | (/ 1 (/ (* 2 a) (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))))) |
| 64.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 61.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 38.0ms | (- (pow b 2) (* 4 (* a c))) |
| 396.0ms | (sqrt (- (pow (pow b 2) 3) (pow (* 4 (* a c)) 3))) |
| 266.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (- (pow b 2) (* 4 (* a c)))) (sqrt (- (pow b 2) (* 4 (* a c))))) (* (- b) (sqrt (- (pow b 2) (* 4 (* a c))))))) |
| 166.0ms | (sqrt (- (* (pow b 2) (pow b 2)) (* (* 4 (* a c)) (* 4 (* a c))))) |
| 125.0ms | (- b (* 2 (/ (* a c) b))) |
| 90.0ms | (cbrt (/ (* 2 a) (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))))) |
10 alts after pruning (8 fresh and 2 done)
Merged error: 4.9b
Found 4 expressions with local error:
| 24.0b | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 24.0b | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 24.0b | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 1.2b | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 14× | add-sqr-sqrt |
| 10× | sqrt-prod |
| 9× | *-un-lft-identity |
| 7× | add-cube-cbrt |
| 6× | add-log-exp |
| 6× | sqrt-div |
| 4× | flip-- |
| 4× | add-exp-log |
| 4× | add-cbrt-cube |
| 4× | flip3-- |
| 4× | pow1 |
| 3× | pow1/2 |
| 3× | rem-sqrt-square |
| 2× | difference-of-squares |
| 1× | distribute-lft-out-- |
| 1× | diff-log |
| 1× | sub-neg |
| 1× | neg-sub0 |
| 1× | associate--l- |
| 14.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 6.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 6.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 6.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 59.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 59.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 55.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 40.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 436.0ms | (sqrt (- (pow (pow b 2) 3) (pow (* 4 (* a c)) 3))) |
| 371.0ms | (sqrt (- (pow (pow b 2) 3) (pow (* 4 (* a c)) 3))) |
| 358.0ms | (sqrt (- (pow (pow b 2) 3) (pow (* 4 (* a c)) 3))) |
| 222.0ms | (sqrt (- (* (pow b 2) (pow b 2)) (* (* 4 (* a c)) (* 4 (* a c))))) |
| 181.0ms | (sqrt (- (* (pow b 2) (pow b 2)) (* (* 4 (* a c)) (* 4 (* a c))))) |
10 alts after pruning (7 fresh and 3 done)
Merged error: 4.9b
83.2% (4.7b remaining)
Error of 10.4b against oracle of 5.7b and baseline of 33.5b
| 11455× | (pre true 80) |
| 5803× | (body real 80) |
| 2738× | (body nan 80) |
| 1107× | (body real 2560) |
| 947× | (body real 1280) |
| 433× | (body real 640) |
| 217× | (body real 320) |
| 109× | (body real 5120) |
| 101× | (body real 160) |