



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 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 5.0
Final simplification10.4
herbie shell --seed 2019016
(FPCore (a b c)
:name "quadm (p42, negative)"
: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: 21.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) |
| 113.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- |
| 21.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 6.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 2.0ms | (- (* b b) (* 4 (* a c))) |
| 2.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 76.0ms | (/ (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) (* 2 a)) |
| 69.0ms | (- (- b) (sqrt (- (* b b) (* 4 (* a c))))) |
| 65.0ms | (sqrt (- (* b b) (* 4 (* a c)))) |
| 27.0ms | (- (* b b) (* 4 (* a c))) |
| 319.0ms | (sqrt (- (pow (* b b) 3) (pow (* 4 (* a c)) 3))) |
| 274.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* 4 (* a c)))) (sqrt (- (* b b) (* 4 (* a c)))))) |
| 199.0ms | (- (pow (- b) 3) (pow (sqrt (- (* b b) (* 4 (* a c)))) 3)) |
| 172.0ms | (sqrt (- (* (* b b) (* b b)) (* (* 4 (* a c)) (* 4 (* a c))))) |
| 129.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- |
| 39.0ms | (/ (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) (* 2 a)) |
| 18.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 7.0ms | (- (pow b 2) (* 4 (* a c))) |
| 4.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 58.0ms | (/ (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) (* 2 a)) |
| 54.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 46.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 31.0ms | (- (pow b 2) (* 4 (* a c))) |
| 973.0ms | (- (* (- b) (- b)) (* (sqrt (- (pow b 2) (* 4 (* a c)))) (sqrt (- (pow b 2) (* 4 (* a c)))))) |
| 294.0ms | (sqrt (- (pow (pow b 2) 3) (pow (* 4 (* a c)) 3))) |
| 280.0ms | (- (pow (- b) 3) (pow (sqrt (- (pow b 2) (* 4 (* a c)))) 3)) |
| 165.0ms | (sqrt (- (* (pow b 2) (pow b 2)) (* (* 4 (* a c)) (* 4 (* a c))))) |
| 98.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- |
| 7.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 7.0ms | (/ 1 (/ (* 2 a) (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))))) |
| 3.0ms | (- (pow b 2) (* 4 (* a c))) |
| 3.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 64.0ms | (/ 1 (/ (* 2 a) (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))))) |
| 52.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 40.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 25.0ms | (- (pow b 2) (* 4 (* a c))) |
| 281.0ms | (sqrt (- (pow (pow b 2) 3) (pow (* 4 (* a c)) 3))) |
| 223.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))))))) |
| 151.0ms | (sqrt (- (* (pow b 2) (pow b 2)) (* (* 4 (* a c)) (* 4 (* a c))))) |
| 95.0ms | (/ (* 2 a) (- (- b) (sqrt (- (pow b 2) (* 4 (* a c)))))) |
| 93.0ms | (- b (* 2 (/ (* a c) b))) |
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- |
| 7.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 3.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 3.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 3.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 72.0ms | (- (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 54.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 53.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 45.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 367.0ms | (sqrt (- (pow (pow b 2) 3) (pow (* 4 (* a c)) 3))) |
| 342.0ms | (sqrt (- (pow (pow b 2) 3) (pow (* 4 (* a c)) 3))) |
| 341.0ms | (sqrt (- (pow (pow b 2) 3) (pow (* 4 (* a c)) 3))) |
| 168.0ms | (sqrt (- (* (pow b 2) (pow b 2)) (* (* 4 (* a c)) (* 4 (* a c))))) |
| 168.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) |