



Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 33.2 |
|---|---|
| Target | 20.7 |
| Herbie | 8.7 |
if b < -4.8021396801207425e+125Initial program 50.9
Taylor expanded around inf 50.9
Taylor expanded around -inf 3.1
if -4.8021396801207425e+125 < b < 1.13066952980576e-152Initial program 10.9
Taylor expanded around inf 10.8
if 1.13066952980576e-152 < b < 3.0304262084228798e+63Initial program 37.3
rmApplied flip-+37.4
Applied associate-/l/40.6
Simplified20.1
rmApplied associate-/r*15.9
if 3.0304262084228798e+63 < b Initial program 56.5
Taylor expanded around inf 56.5
Taylor expanded around inf 3.5
Simplified3.5
Final simplification8.7
herbie shell --seed 2019016
(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: 24.1s)Debug log
| 378× | (pre true 80) |
| 187× | (body real 80) |
| 93× | (body nan 80) |
| 36× | (body real 2560) |
| 30× | (body real 1280) |
| 16× | (body real 640) |
| 8× | (body real 160) |
| 6× | (body real 320) |
| 2× | (body real 5120) |
| 67.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
1 alts after pruning (1 fresh and 0 done)
Merged error: 35.0b
Found 4 expressions with local error:
| 23.2b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 2.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)) |
| 13× | *-un-lft-identity |
| 8× | add-log-exp |
| 8× | add-sqr-sqrt |
| 7× | add-cube-cbrt |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 4× | associate-/l* |
| 4× | times-frac |
| 4× | pow1 |
| 3× | distribute-lft-out |
| 3× | sqrt-prod |
| 2× | associate-/l/ |
| 2× | flip-- |
| 2× | flip-+ |
| 2× | associate-+l- |
| 2× | sqrt-div |
| 2× | flip3-- |
| 2× | flip3-+ |
| 2× | neg-sub0 |
| 1× | div-inv |
| 1× | div-exp |
| 1× | div-sub |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | sum-log |
| 1× | associate-/r* |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 1× | cbrt-undiv |
| 36.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 12.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 4.0ms | (- (* b b) (* (* 4 a) c)) |
| 4.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 78.0ms | (/ (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (* 2 a)) |
| 72.0ms | (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 45.0ms | (- (* b b) (* (* 4 a) c)) |
| 42.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 452.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 289.0ms | (- (* (- b) (- b)) (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c))))) |
| 244.0ms | (+ (pow (- b) 3) (pow (sqrt (- (* b b) (* (* 4 a) c))) 3)) |
| 204.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
| 173.0ms | (/ (* (* (+ (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (+ (- b) (sqrt (- (* b b) (* (* 4 a) c))))) (+ (- b) (sqrt (- (* b b) (* (* 4 a) c))))) (* (* (* 2 a) (* 2 a)) (* 2 a))) |
8 alts after pruning (8 fresh and 0 done)
Merged error: 5.2b
Found 4 expressions with local error:
| 23.2b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 16.3b | (/ (* 4 (* c a)) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 1.2b | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 1.2b | (- (* b b) (* (* 4 a) c)) |
| 10× | add-sqr-sqrt |
| 9× | add-exp-log |
| 9× | add-cbrt-cube |
| 8× | add-log-exp |
| 7× | *-un-lft-identity |
| 5× | add-cube-cbrt |
| 4× | flip-- |
| 4× | sqrt-prod |
| 4× | flip3-- |
| 4× | pow1 |
| 2× | difference-of-squares |
| 2× | associate-/r/ |
| 2× | associate-*r/ |
| 2× | div-exp |
| 2× | diff-log |
| 2× | sqrt-div |
| 2× | sub-neg |
| 2× | cbrt-undiv |
| 1× | distribute-lft-out-- |
| 1× | associate-/l* |
| 1× | div-inv |
| 1× | cbrt-unprod |
| 1× | prod-exp |
| 1× | pow1/2 |
| 1× | times-frac |
| 1× | frac-2neg |
| 1× | associate-/r* |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | neg-sub0 |
| 1× | associate--l- |
| 25.0ms | (/ (* 4 (* c a)) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 10.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 4.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 4.0ms | (- (* b b) (* (* 4 a) c)) |
| 157.0ms | (/ (* 4 (* c a)) (* (* 2 a) (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 66.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 59.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 27.0ms | (- (* b b) (* (* 4 a) c)) |
| 1.2s | (/ (* (* (* 4 (* c a)) (* 4 (* c a))) (* 4 (* c a))) (* (* (* (* 2 a) (* 2 a)) (* 2 a)) (* (* (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))) (- (- b) (sqrt (- (* b b) (* (* 4 a) c))))))) |
| 784.0ms | (/ (* (* (* 4 (* c a)) (* 4 (* c a))) (* 4 (* c a))) (* (* (* (* 2 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))))))) |
| 506.0ms | (+ (* (- b) (- b)) (+ (* (sqrt (- (* b b) (* (* 4 a) c))) (sqrt (- (* b b) (* (* 4 a) c)))) (* (- b) (sqrt (- (* b b) (* (* 4 a) c)))))) |
| 402.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 177.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
7 alts after pruning (7 fresh and 0 done)
Merged error: 5.2b
Found 4 expressions with local error:
| 23.2b | (sqrt (- (* b b) (* (* 4 a) c))) |
| 13.9b | (/ (* 4 (* c a)) (* 2 a)) |
| 1.2b | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 1.2b | (- (* b b) (* (* 4 a) c)) |
| 10× | add-sqr-sqrt |
| 8× | add-log-exp |
| 7× | *-un-lft-identity |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 5× | add-cube-cbrt |
| 4× | sqrt-prod |
| 4× | pow1 |
| 3× | flip-- |
| 3× | flip3-- |
| 2× | difference-of-squares |
| 2× | diff-log |
| 2× | sqrt-div |
| 2× | sub-neg |
| 1× | distribute-lft-out-- |
| 1× | associate-/l* |
| 1× | div-inv |
| 1× | div-exp |
| 1× | pow1/2 |
| 1× | times-frac |
| 1× | frac-2neg |
| 1× | associate-/r* |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | cbrt-undiv |
| 1× | neg-sub0 |
| 1× | associate--l- |
| 16.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 7.0ms | (/ (* 4 (* c a)) (* 2 a)) |
| 4.0ms | (- (* b b) (* (* 4 a) c)) |
| 4.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 69.0ms | (sqrt (- (* b b) (* (* 4 a) c))) |
| 61.0ms | (- (- b) (sqrt (- (* b b) (* (* 4 a) c)))) |
| 41.0ms | (- (* b b) (* (* 4 a) c)) |
| 8.0ms | (/ (* 4 (* c a)) (* 2 a)) |
| 828.0ms | (/ (* (* (* 4 (* c a)) (* 4 (* c a))) (* 4 (* c a))) (* (* (* 2 a) (* 2 a)) (* 2 a))) |
| 450.0ms | (- (log (* 4 (* c a))) (log (* 2 a))) |
| 423.0ms | (sqrt (- (pow (* b b) 3) (pow (* (* 4 a) c) 3))) |
| 189.0ms | (sqrt (- (* (* b b) (* b b)) (* (* (* 4 a) c) (* (* 4 a) c)))) |
| 137.0ms | (- b (* 2 (/ (* a c) b))) |
7 alts after pruning (6 fresh and 1 done)
Merged error: 5.2b
Found 4 expressions with local error:
| 23.2b | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 2.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)) |
| 13× | *-un-lft-identity |
| 8× | add-log-exp |
| 8× | add-sqr-sqrt |
| 7× | add-cube-cbrt |
| 6× | add-exp-log |
| 6× | add-cbrt-cube |
| 4× | associate-/l* |
| 4× | times-frac |
| 4× | pow1 |
| 3× | distribute-lft-out |
| 3× | sqrt-prod |
| 2× | associate-/l/ |
| 2× | flip-- |
| 2× | flip-+ |
| 2× | associate-+l- |
| 2× | sqrt-div |
| 2× | flip3-- |
| 2× | flip3-+ |
| 2× | neg-sub0 |
| 1× | div-inv |
| 1× | div-exp |
| 1× | div-sub |
| 1× | diff-log |
| 1× | pow1/2 |
| 1× | frac-2neg |
| 1× | sub-neg |
| 1× | sum-log |
| 1× | associate-/r* |
| 1× | clear-num |
| 1× | rem-sqrt-square |
| 1× | +-commutative |
| 1× | cbrt-undiv |
| 23.0ms | (/ (+ (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) (* 2 a)) |
| 19.0ms | (+ (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 4.0ms | (- (pow b 2) (* 4 (* a c))) |
| 4.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 80.0ms | (/ (+ (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) (* 2 a)) |
| 68.0ms | (+ (- b) (sqrt (- (pow b 2) (* 4 (* a c))))) |
| 41.0ms | (sqrt (- (pow b 2) (* 4 (* a c)))) |
| 26.0ms | (- (pow b 2) (* 4 (* a c))) |
| 1.1s | (- (* (- b) (- b)) (* (sqrt (- (pow b 2) (* 4 (* a c)))) (sqrt (- (pow b 2) (* 4 (* a c)))))) |
| 674.0ms | (+ (pow (- b) 3) (pow (sqrt (- (pow b 2) (* 4 (* a c)))) 3)) |
| 439.0ms | (sqrt (- (pow (pow b 2) 3) (pow (* 4 (* a c)) 3))) |
| 174.0ms | (sqrt (- (* (pow b 2) (pow b 2)) (* (* 4 (* a c)) (* 4 (* a c))))) |
| 128.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))) |
7 alts after pruning (5 fresh and 2 done)
Merged error: 5.2b
88.4% (3.0b remaining)
Error of 8.7b against oracle of 5.7b and baseline of 31.9b
| 11500× | (pre true 80) |
| 5932× | (body real 80) |
| 2746× | (body nan 80) |
| 1013× | (body real 2560) |
| 888× | (body real 1280) |
| 455× | (body real 640) |
| 223× | (body real 320) |
| 126× | (body real 160) |
| 117× | (body real 5120) |